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If I can't picture it, I can't understand it. 

You know, it would be sufficient to really understand the electron. 



A. Einstein 



• If a spinning particle is not quite a point particle, nor a solid three dimensional top, what 
can it be? What is the structure which can appear under probing with electromagnetic 
fields as a point charge, yet as far as spin and wave properties are concerned exhibits a 
size of the order of the Compton wavelength? 



A.O. Barut 
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Preface 



The present notes contain some basic materials, physical and mathematical, of the general 
formalism for analyzing elementary particles, which under the general name of Kinemati- 
cal Formalism of Elementary Spinning Particles, I have been working during the last years. 
The name kinematical makes reference to its close relationship with the kinematical group of 
space-time transformations associated to the Restricted Relativity Principle which a theoretical 
framework must necessarily satisfy. 

In a certain sense it is a revision of the basic fundamentals of the Lagrangian formalism 
which leads to Euler-Lagrange equations, Noether's theorem, etc., but looking for solutions 
which go through the postulated initial and final states of the variational formalism. This 
produces a classical formalism which is going to be expressed in terms of the end point variables 
of the dynamical evolution. This formalism is, therefore, closer to the quantum mechanical 
dynamical theory and it is through Feynman's path integral approach that we can find the 
bridge between them. 

These end point variables, which I propose to call them kinematical variables, in the case 
of elementary particles will necessarilly span a homogeneous space of the kinematical group. 
In this way, the kinematical group not only refiects the symmetries of the system. It also 
supplies the necessary variables to describe elementary matter. It is crucial for the description 
of matter to improve in our knowledge of this kinematical group. In the present notes we shall 
deal mainly with the Galilei and Poincare groups, but the formalism is so general that it can 
accomodate any further group we consider as the basic symmetry group of matter. 

The notes pretend to be selfcontained and in this way we have included at the end of 
the first three chapters some mathematical appendices which contain not very well spread 
materials. The lecture notes are organised as follows. We begin with an enumeration of some 
pros and cons the formalism has. Probably the right place of this section will be at the end 
of the notes as a kind of conclusions and comments. I like to put it first because it suggests 
many of the features we can find when analyzing the classical models and may create some 
provocation in the reader. The remaining materials are collected into five chapters which more 
or less will cover one lecture each. 

We begin with the general Lagrangian formalism just to enhance the role of the kinematical 
variables in defining a concept of elementary particle. Lecture two will be devoted to the 
analysis of several relativistic and nonrelativistic models, to show how the standard methods 
of analyzing symmetries leads to the definition of the relevant observables. In particular, we 
shall pay attention to the definition of the spin. The spin, as any other observable, will be 
defined in the classical case in terms of the degrees of freedom and their derivatives, and we 
shall analyze its mathematical structure. 

The next two lectures will cover the quantization of the formalism and the analysis of some 
relativistic and nonrelativistic examples. The separate fourth chapter is devoted to the model 
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which satisfies Dirac's equation. Special attention is paid to the analysis of Dirac's algebra 
and its relationship with the classical observables and to show a geometrical interpretation of 
the difi^erence in chirality between matter and antimatter. 

Finally, some physical features which are related to the spin of the elementary particles, 
are described. In some places, the lectures will be complemented with numerical simulations 
whenever the theoretical solution is not available or very difficult to interpret because of the 
mathematical complexity. 

I am very grateful to professor Oleg Teryaev for his kind invitation to the Joint Institute 
for Nuclear Research at Dubna. 



Martin Rivas 



Bilbao-Dubna, September 2005. 



Pros and cons of the kinematical 
formalism 



I will present in what follows a kind of general introduction to the kinematical formalism, 
in the form of some pros and cons the formalism has, from my point of view, with some 
additional comments. There will be probably many more cons than the ones quoted, and the 
pros I consider could be nonsense for other readers. You will probably put some of the pros in 
the cons part. Do it, please. Science is a collective job and probably the greatest endeavour of 
mankind. Any new proposal needs a thorough analysis to find its contradictions. It is the only 
way to improve knowledge. I am strongly convinced that the kinematical formalism I propose 
has more advantages than previous approaches and it is simpler. It opens new perspectives to 
establish a deeper formalism to deal with spinning particles from the very begining. 

Things like the prediction of formation of electromagnetic bound pairs of electrons, the 
chiral difference between particles and antiparticles or the justification of gi = 2 by pure 
kinematical arguments are sufficient, from my point of view, to deserve some time to it, even 
to analyze its contradictions. One feature of the formalism is that it is not finished yet. Another 
is that it shows that the space-time symmetry group of the described elementary particles is 
larger than the group we start with. We have to accomplish, and finish, this task. 

We accept a variational formalism for describing the dynamics of elementary particles. 
To be consistent with the variational statements we have to look for solutions of the corre- 
sponding differential equations, passing through the fixed end points. This implies that the 
variables which define the end points of the formalism, which will be called from now on kine- 
matical variables, will play in the classical case the same role as the wave function in the 
quantum formalism, to characterising the states of the system. Our first task will be to rewrite 
the variational formalism in terms of these kinematical variables. It is for this reason we have 
included a first lecture about a generalized Lagrangian formalism. 

PROS 

• An elementary particle is by definition the simplest mechanical system. An elementary 
particle can be annihilated but it can never be deformed. Its intrinsic attributes are 
not modified by any interaction. It has no excited states and therefore all its possible 
states are just kinematical modifications of any one of them. If the state of the particle 
changes it is always possible to find a new inertial observer who describes the particle 
in the same state as before. In the quantum case this leads to the conclusion that the 
Hilbert space which describes its pure states carries an irreducible representation 
of the kinematical group. In the classical case the kinematical Lagrangian space is a 
homogeneous space of the kinematical group. 
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It is the Restricted Relativity Principle which characterizes not only the space-time 
symmetry group of the theory, usually called the kinematical group, but also which 
supplies the classical variables to describe an elementary spinning particle. 

The formalism is independent of the kinematical group. The classical variables which 
characterize the kinematical state of the system are not postulated. They are related to 
the variables which characterize the parameterization of the kinematical group, and the 
manifold they span is necessarily a homogeneous space of the kinematical group. 

There is no need to use Grassmann or spinor variables to characterize the states of a 
classical elementary spinning particle. If we restrict ourselves to the Galilei or Poincare 
groups, the only and most general variables to characterize the states of a classical 
elementary particle are the 10 variables t, r, v, cx, interpreted respectively as the time, 
position of a point where the charge is located, velocity of this point and finally the 
orientation of the system around this point. These variables are in fact, the variables 
which define a parametrization of any of the mentioned groups. 

The point particle, either relativistic or non-relativistic, is an elementary particle accord- 
ing to this formalism. It is not postulated, it is a consequence of the formalism. Its 
initial and final states are just characterized by the kinematical variables t and r. It 
is the simplest localized system the formalism allows to describe but it represents an 
spinless object. It seems that there are no spinless elementary particles in nature and 
therefore the use of spinless point particles for the description of physical phenomena 
is an approximate one. If the whole formalism of particle physics starts by dealing at 
first with spinless point particles and afterwards this system is endowed with spin and 
other internal properties in some ad hoc manner, this will produce a different framework 
than starting with spinning elementary particles from the very beginning. There is a 
quotation by Albert Einstein that: "Things should be made as simple as possible, but 
not simpler." 

The spin is an angular momentum and therefore it is always defined with respect to some 
fixed point. This point has to be clearly identified. 

If we understand as the spin of the electron as the angular momentum of the electron 
with respect to its center of mass, then the spin of the electron is twofold. One part is 
related to the rotation of the electron (a variables and its derivative the angular velocity 
uj) and the other to the separation between the position of the charge and the center of 
mass, which is a different point, and its relative orbital motion. This comes from the use 
of the V variables and its derivative, the acceleration a. The rotation does not produce 
magnetic moment. The magnetic moment is the result of the motion of the charge around 
the center of mass, so that when the magnetic moment is expressed in terms of the total 
spin this produces a clear interpretation of the gyromagnetic ratio. 

The size and shape of any elementary particle are approximate observables. They are 
geometrical aspects which can be defined for macroscopic objects but they loose their 
geometrical meaning when talking about elementary particles. 

The center of mass of the spinning particle, which is a very well defined point from the 
classical point of view, can be related to the Newton- Wigner position operator or to the 
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position vector defined by the Foldy-Wouthuysen transformation, in the quantum case. 
Dirac spin operator is not the angular momentum of the electron with respect to its 
center of mass. It is the angular momentum with respect to the center of charge. 

• The photon and the charge of the electron move at the speed of light. This motion is not 
altered by any external interaction. Only the center of mass motion is affected. This is 
not contradictory with special relativity because the center of mass or center of energy 
moves always at a velocity v < c. 

• For any inertial observer the charge of the electron is never at rest and therefore it always 
measures magnetic moment and electric dipole moment with respect to the center of mass. 

• In a certain sense electromagnetism means speed of light for the carriers of the electro- 
magnetic interaction and also for the sources of the field. 

• The charge of the electron is located at a single point and therefore we have no problems 
associated to the charge distribution of the extended models. But at the same time it 
moves and oscillates in a region of radius half Compton's wavelength, where at this scale 
quantum phenomena appear. 

• In the relativistic case the most invariant way to define the internal structure of the 
electron is precisely to assume that the charge is moving at the speed of light. This 
velocity cannot be altered by any interaction, and therefore this corroborates the idea 
that an elementary particle cannot be deformed. 

• Classical particles which move in straight lines at the speed of light are massless particles 
which rotate around the direction of motion with the spin S either parallel or antiparallel 
to the velocity, but in the same direction as the angular velocity a;. Spin is invariant 
under Poincare transformations but the angular velocity transforms to produce Doppler 
effect. The spin, when quantized, is not restricted. It can take all greater than zero 
values 1/2, 1, 3/2, .... The spinless case is not predicted. The energy of the spin 1 
object is just S ■ u) = hv, being u the frequency of the rotation of the particle. 

• Particles whose charge moves in circles at the speed of light, in a plane orthogonal to 
the spin, although the center of mass is moving below c, satisfy when quantized, Dirac's 
equation. All of them represent massive charged spin 1/2 particles. Mass and charge are 
unrestricted in this formalism. They are the only systems, predicted by this formalism, 
which satisfy Dirac equation when quantized. 

• Particles whose charge moves faster than light, although the center of mass is moving 
below c, all of them represent, when quantized, massive charged spin 1 particles. Mass 
and charge are also unrestricted. Is it a plausible description of classical massive 
bosons? 

• Elementary spinning particles show a clear chirality. Once the spin direction is fixed, 
the position of the charge of the particle has a definite direction of motion, while the 
antiparticle has the opposite one. They are mirror images of each other. 

• The electron is just a moving charge which can also rotate. The electric and magnetic 
moment are not intrinsic properties. They are derived observables produced by the sep- 
aration between the center of mass and the position of the charge and the motion of the 
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charge, respectively. This justifies that quantum electrodynamics is completely deter- 
mined by the minimal coupling between the electron current and the external potentials, 
and without any anomalous electric or magnetic coupling. 

Dirac algebra of 4 x 4 matrices is completely generated by complex linear combination 
of 16 linearly independent hermitian matrices which represent very precise observables. 
Their real linear combinations will produce all translation invariant observables of the 
electron. One of them is the unit matrix and the other are 15 traceless hermitian matrices 
which represent the following observables: 3 are the spin components in the laboratory 
frame, other 3 are the spin components in the body frame and the remaining 9 are 
the 9 components of the three unit vectors of the body frame. We must remember 
that the elementary particle has as degrees of freedom the position of a point and its 
orientation in space which can be characterized by the description of the body frame. 
Once the translational degrees of freedom are supressed, by analyzing for instance the 
electron in the center of mass frame, the orientation observables completely characterize 
its translation invariant structure. From the algebraic point of view Dirac algebra can 
be generated by products and linear combinations of any 4 of these 9 components of its 
body frame. 

When analized the interaction between two spinning electrons different processes can 
be described according to the energy involved and the separation between particles. In 
low energy we can describe elastic scattering and also the possibility of formation of a 
spin 1 bound state of two electrons, provided the energy is below a certain value, the 
spins are parallel and their respective center of masses are separated by a distance below 
Compton's wavelength. In high energy, for relative velocities of the center of masses of 
the particles greater than 0.1c, deep inelastic scattering processes appear, where by deep 
we mean that the two charges approach each other below Compton's wavelength. 

The bound state of two equal charged particles is metastable. The mass of the system is 
greater than the sum of the two masses of the elementary particles. 

The formalism is complete in the sense that all Lagrangian systems whose kinematical 
space is a homogeneous space of the kinematical group can be quantized and the cor- 
responding Hilbert space carries a projective unitary irreducible representation of the 
kinematical group. All known one-particle wave equations are produced by quantization 
of the corresponding Lagrangian models. 

To quantize these systems it is not necessary to use any constrained Hamiltonian formal- 
ism like the one proposed by Dirac, because for every system we have a very well defined 
non-singular Lagrangian, as suggested by Feynman, and the quantization is performed 
through Feynman's path integral approach. 

The formalism predicts, if PCT invariance is assumed, that the magnetic moment and 
spin of an elementary particle and its antiparticle must necessarily have the same relative 
orientation, either parallel or antiparallel. It is postulated that electrons and positrons 
have opposite relative orientation. To my knowledge no clear experimental evidence of 
this relative measurement for free electrons and positrons have ever been performed. 
The same thing happens for /x"*" and /i' . All very accurate measurements of g — 2 are 
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precession experiments which do not discriminate whether spin and magnetic moment 
are parallel or antiparallel. 

• The description of the electron can be done in terms of dimensionless variables, so that 
the symmetry group of the model is larger than the Poincare group. Space-time dila- 
tions are among the new symmetries. Because the description of the orientation has to 
be independent of how we choose the local inertial body frame, we have at least, in the 
quantum formulation an additional SU{2) x U{1) kinematical group of space-time trans- 
formations. What is the physical meaning of the generators of these groups? Although 
we started the formalism by assuming Poincare invariance we find that the elementary 
spinning objects it describes have a larger kinematical group of symmetries. We have 
to start again, so that the new group will give us new classical variables to describe 
more internal structure. But this additional 'internal' structure is related to space-time 
symmetries. Poincare invariance only describes mass and spin as intrinsic (observer in- 
dependent) properties. It is through the analysis of the Casimir operators of the new 
enlarged group that new intrinsic properties could be defined. This new group does not 
commute with the Poincare group so that the enlargement is not simply a direct product 
of the two groups. This has not been completed yet. 



CONS 

• We consider massive particles whose charge is moving at the speed of light. The charge 
of a massive elementary spinning particle is located at a point. The classical system 
that when quantized satisfies Dirac equation is that one whose charge moves at the 
speed of light. This is not contradictory with special relativity because no information 
is traveling faster than light for distances greater than Compton's wavelength, which is 
the spatial domain for the quantum phenomena. Energy moves with velocity below c. 
This possibility of motion at the speed of light is also contained in Dirac's analysis of the 
electron. 

• The charge of a free electron moves in circles at the speed of light in the center of 
mass frame. It is therefore accelerated and classical mechanics predicts radiation. But 
classical mechanics predicts that this internal motion of the charge is stationary, and 
therefore quantum mechanics allows nonradiating stationary states. The classical theory 
of radiation has to be revisited and radiation has to be associated with the acceleration 
of the center of mass. 

• The electric field created by the electron is neither static nor Coulomb like. It behaves 
like 1/r from the retarded charge position. Nevertheless, the time average value of the 
electric field over one complete turn of the charge, is static and Coulomb like from around 
a distance of five Compton's wavelength from the center of mass up to infinity and in 
any direction. It does not diverge at the origin where it goes to zero. The instantaneous 
electromagnetic energy density in the surroundings of the charge goes like and is 
greater than that for the strict Coulomb field which goes like 1/r''. 
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The only divergences of the electromagnetic field appear in the zitterbewegung plane but 
all they go like 1/r when r — > and therefore the electromagnetic energy density has 
no divergence there. Although they are promissing features I have not been able for the 
moment, to renormalize these classical models. 

The position of the charge of the electron satisfies a fourth order differential equation. 
This fourth order differential equation is most difficult to analyze. But it can be separated 
in a system of coupled second order differential equations for the center of mass and center 
of charge. One possibility is that we have Newton-like differential equations for the center 
of mass in terms of the external force and a harmonic oscillator like equations for the 
motion of the charge around the center of mass, and which is independent of the external 
force, thus confirming the idea that the internal structure is unaffected by the interaction. 
Nevertheless a fourth order differential equation is the most general differential equation 
a point can satisfy in three-dimensional space, as derived from Frenet-Serret equations. 
External interactions on a spinning particle will produce forces and torques. It will 
therefore change the curvature and torsion of the trajectory, thus justifying the need of 
a fourth order differential equation. 

It is usually said that a point cannot rotate. It is clear that all matter that surround 
us moves and rotates. In the usual classical approach, it seems that the only exception 
are the elementary particles. It is out of logic that the basic constituents of matter are 
excluded from rotation. Why? This implies that the consideration of the point particle 
(and therefore spinless) as the starting object to build all material systems as is done 
in Newtonian mechanics is only an approximate formalism. In our formalism we have 
to attach orientation to the point to describe its possible rotation. And this rotation 
contributes to the spin of the system, although there is another contribution coming 
from the orbital motion of the charge around the center of mass. 
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Lagrangian formalism 



1.1 Generalized Lagrangian formalism 

The Lagrangian formalism of generalized systems depending on higher order derivatives 
was already worked out by Ostrogradsky. ^ We shall outline it briefly here, mainly to analyze 
the generalized Lagrangians not only in terms of the independent degrees of freedom but also 
as functions of what we shall call the kinematical variables of the system, i.e., of the end 
point variables of the variational formulation. 

Let us consider a mechanical system of n degrees of freedom, characterized by a Lagrangian 
that depends on time t and on the n essential coordinates qiit), that represent the n indepen- 
dent degrees of freedom, and their derivatives up to a finite order k. Because we can have 
time derivatives of arbitrary order we use a superindex enclosed in brackets to represent the 
corresponding fc-th derivative, i.e., q"^\t) = d!' qi{t) / dt^ . The action functional is defined by: 

A[q] = r Lit,q,it),qi'\t),...,qf'\t))dt, (1.1) 

where i = 1, . . . ,n. Using a more compact notation we define q^^^ = g^, and therefore we shall 
write 

L(t,g,(t),g«(t),...,gf (t))=L(t,gf'(^)), 

for s = 0, . . . , A;. 

The trajectory followed by the mechanical system is that path which passing through 
the fixed end-points at initial and final times ti and t2, qi^\ti) and ql^\t2), ^ = 1, ...,n, s = 
0,l,...,/c — 1, makes extremal the action functional (jl.lll . Note that we need to fix as boundary 
values of the variational principle some particular values of time t, the n degrees of freedom 
and their derivatives up to order k — 1, i.e., one order less than the highest derivative of each 
variable in the Lagrangian, at both end points of the problem. In other words we can say 
that the Lagrangian of any arbitrary generalized system is in general an explicit function of 
the variables we keep fixed as end points of the variational formulation and also of their next 
order time derivative. 

Once the action functional is defined for some particular path qi {t) , to analyze its variation let 
us produce an infinitesimal modification of the functions qi (t) , qi (t) qi (t) + Sqi (t) while leaving 

^ M. Ostrogradsky, Memoire sur les equations difFerentielles relatives au probleme des isoperimetres, Mem. 
Acad. St. Petersburg, 6(4), 385-517 (1850). 
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fixed the end-points of the variational problem, i.e., such that at ti and t2 the modification of 
the generalized coordinates and their derivatives up to order fc — 1 vanish, and thus Sq^'^\tl) = 
^Qi^\t2) — 0, for i — 1, . . . ,n and s = 0, 1, . . . , fc — 1. Then, the variation of the derivatives of 
the qi(t) is given by q'l'^^t) qi^\t) + Sq'l^\t) — (7-''''(t) + d'^ Sqi{t) / df^ , since the modification of 
the s-th derivative function is just the s-th derivative of the modification of the corresponding 
function. This produces a variation in the action functional 5 A = A[q + Sq] — A[q], given by: 



SA^ 



L{t, ql'\t) + Sq^'\t))dt - / L{t, <lt\t))dt 



dtj: 



^^X (1) . , 9L (fc) 



after expanding to lowest order the first integral. The term 
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dL 



dq. 
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dq. 
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dL 



dt 



.dq.. 
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Sqi, 



and by partial integration of this expression between ti and t2, it gives: 



ti dq. 



. dt 



dL 



dq. 



(1) 



Sqidt 



dL 



ti 



dt 



(1) 



Sqi dt, 



because the variations dqi{ti) and Sqi{t2), vanish. Similarly for the next term: 
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6qi dt, 



because Sqi and Sq^^^ vanish at ti and t2, and finally for the last term 



dL 



= (-1)'/ jp: 



t2 



dL 



dq. 



(k) 



Sqi dt, 



(1.2) 



so that each term of H1.2|l is written only in terms of the variations of the degrees of freedom 6qi 
and not of their higher order derivatives. Remark that to reach these final expressions, it has 
been necessary to assume the vanishing of all Sq'^'^\ for s = 0, . . . , fc — 1, at times ti and t2. By 
collecting all terms we get 



6A^ 



dtj: 



dL_d_ 
dqi dt 



dL 



Jk 



dL 



6q^. 



If the action functional is extremal along the path qi{t), its variation must vanish, SA = 0. The 
variations Sqi are arbitrary and therefore all terms between squared brackets cancel out. We 
obtain a system of n differential equations. 



dL 



dq^ dt \gq 



dL 



(1) 



+ ••• + (-1)* 



dL 



dt^ 



(fc) 



0, 1^1,. 



(1.3) 



the Euler-Lagrange equations, which can be written in condensed form as: 

k 



dL 



s=0 



0, 



i = 1, . . . ,n. 



(1.4) 
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1.2 Kinematical variables 

In general, the system (II .41) is a system of n ordinary differential equations of order 2k, 
and thus existence and uniqueness theorems guarantee only the existence of a solution of this 
system for the 2kn boundary conditions ql^\ti), i = 1, . . . ,n and s = 0, 1, . . . , 2A; — 1, at the 
initial instant ti. However the variational problem has been stated by the requirement that 
the solution goes through the two fixed endpoints, a condition that does not guarantee neither 
the existence nor the uniqueness of the solution. Nevertheless, let us assume that with the 
fixed endpoint conditions of the variational problem, and g-'^^(t2), i = l,...,n and 

s = 0, 1, . . . , /c — 1, at times ti and ^2, respectively, there exists a solution of ()1.4p perhaps non- 
unique. This implies that the 2kn boundary conditions at time ti required by the existence 
and uniqueness theorems, can be expressed perhaps in a non-uniform way, as functions of the 
kn conditions at each of the two endpoints. From now on, we shall consider systems in which 
this condition is satisfied. It turns out that a particular solution passing through these points 
will be expressed as a function of time with some explicit dependence of the end point values 



j, / = 1, . . . , n, r = 0, 1, ... A; — 1, in terms of these boundary end point conditions. 

Definition: The Action Function ^ of the system along a classical path is the 
value of the action functional ()1.1|1 when we introduce in the integrand a particular 
solution (II. 5|) passing through those endpoints: 



Once the time integration is performed, we see that it will be an explicit function of the 
kn + 1 variables at the initial instant, qj'\ti), r = 0, . . . , k — 1 including the time ti, and of 
the corresponding kn + 1 variables at final time ^2- We write it as 



We thus arrive at the following 

Definition: The kinematical variables of the system are the time t and the n 
degrees of freedom and their time derivatives up to order k — 1. The manifold 
X they span is the kinematical space of the system. 

The kinematical space for ordinary Lagrangians is just the configuration space spanned by 
variables qi enlarged with the time variable t. It is usually called the enlarged configuration 

^Please remark that we use the same letter A{ ) for the action function, followed by normal brackets 
containing the variables of which it depends, and for the action functional A[ ] which is followed by squared 
brackets to enhance that it is not a function but rather a functional. 



m^q^{t■,q^ih),ql^\t,)), 



(1.5) 




(1.6) 
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space. But for generalized Lagrangians it also includes higher order derivatives up to one 
order less than the highest derivative. Thus, the action function of a system becomes 
a function of the values the kinematical variables take at the end points of the 
trajectory, xi and X2. From now on we shall consider systems for which the action function 
is defined and is a continuous and differentiable function of the kinematical variables at the 
end points of its possible evolution. This function clearly has the property A{x,x) = 0. 

1.2.1 Replacement of time as evolution parameter 

The constancy of speed of light in special relativity brings space and time variables on the 
same footing. So, the next step is to remove the time observable as the evolution parameter of 
the variational formalism and express the evolution as a function of some arbitrary parameter 
to be chosen properly. Then, let us assume that the trajectory of the system can be expressed 
in parametric form, in terms of some arbitrary evolution parameter r, {t(r), gj(r)}. The 
functional ()1.H1 can be rewritten in terms of the kinematical variables and their derivatives 
and becomes: 

A[t, q] = £ L (^t{T), g,(r), . . . , iir)dr 

'L{x{r),x{T))dr, (1.7) 

1 

where the dot means derivative with respect to the evolution variable r that without loss of 
generality can be taken dimensionless. Therefore L = L{t{T),qf'' /i{T))i{T) has dimensions of 
action. 

It seems that ()1.7p represents the variational problem of a Lagrangian system depending 
only on first order derivatives and of kn + 1 degrees of freedom. However the kinematical 
variables, considered as generalized coordinates, are not all independent. There exist among 
them the following {k — l)n differential constraints 

qf\T)=i'-'\r)/i{r), z = l,...,n, . = 1, . . . , A; - 1. (1.8) 

We can also see that the integrand L is a homogeneous function of first degree as a function 
of the derivatives of the kinematical variables. In fact, each time derivative function q\^\t) 
has been replaced by the quotient qf^^\T)/i{T) of two derivatives with respect to r. Even the 
highest order /c-th derivative function qf''' = qf''^'' /i is expressed in terms of the derivatives of 
the kinematical variables qf' and t. Thus the original function L is a homogeneous function 
of zero degree in the derivatives of the kinematical variables. Finally, the last term i(r), gives 
to the new defined L the character of a homogeneous function of first degree. Then, Euler's 
theorem on homogeneous functions gives rise to the additional relation: 

L(a;(r),i(r)) = ^i^' = Y.F.j{x,x)xK (1.9) 

i i 

With the above {k — l)n differentiable constraints among the kinematical variables (|1.8I1 and 
condition (II. 9|) . it reduces to n the number of essential degrees of freedom of the system (|1.7jl . 
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This possibility of expressing the Lagrangian as a homogeneous function of first degree of 
the derivatives was already considered in 1933 by Dirac on aesthetical grounds. It is this 
homogeneity of first degree in terms of the derivatives which will allow us later to transform 
the variational formalism into a geodesic problem on the kinematical space X, but where the 
metric gij{x,x) will be direction dependent, and thus the particle trajectory is a geodesic, not 
in a Riemannian manifold but rather in a Finsler space. ^ 

Function L is not an explicit function of the evolution parameter r and thus we can see that 
the variational problem ()1.7p . is invariant with respect to any arbitrary change of evolution 
parameter r. ^ 

In fact, if we change the evolution parameter t = r(6'), then the derivative t(r) = {dt/ d9){d9 / dr) 
and (t) = {dqf\e) / de){de / dr) such that the quotients 

c^f\T) _ {d<^'\e)/de)9{T) ^4'\e) 
i{T) {dt{e)/de)0{T) i{e) ' 

where once again this last dot means derivation with respect to 9. It turns out that lll.7|l can be 
written as: 

= [ ' L[x{e),x{e))de. (i.io) 



1.2.2 Recovering the Lagrangian from the Action function 

The formalism thus stated has the advantage that it is independent of the evolution pa- 
rameter, and if we want to come back to a time evolution description, we just use our time 
as the evolution parameter and make the replacement t = t, and therefore t = 1. From now 
on we shall consider those systems for which the evolution can be described in a parametric 
form, and we shall delete the symbol ^ over the Lagrangian, which is understood as written 
in terms of the kinematical variables and their first order derivatives. 

If what we know is the action function of any system A{xi,X2), as a function of the kine- 
matical variables at the end points we can proceed conversely and recover the Lagrangian 
L{x, x) by the limiting process: 

dA{x,y) ■ 

L(x,x) = hm — — X"*, 1-11 
^ ^ y^x dy^ ^ ^ 

where the usual addition convention on repeated or dummy index j, extended to the whole set 
of kinematical variables, has been assumed. 

If in lll.7|l we consider two very close points xi = x and X2 = x + dx, we have that the action 
function A{x,x-\-dx) = A{x,x + xdT) — L{x,x)dT and making a Taylor expansion of the function 
A with the condition A{x, x) — we get 

^ P. A.M. Dirac, Proc. Cam. Phil. Soc. 29, 389 (1933): "a greater elegance is obtained", "a symmetrical 
treatment suitable for relativity." 

^G.S. Asanov, Finsler geometry, Relativity and Gauge theories, Reidel Pub. Co, Dordrecht (1985). 

^ R. Courant, D. Hilbert, Methods of Mathematical Physics, Vol. 1, Interscience, N.Y. (1970); I.M. Gelfand, 
S.V. Fomin, Calculus of Variations Prentice Hall, Englewood Cliffs, N.J. (1963). 
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1.3 Generalized Noether's theorem 

Noether's analysis for generalized Lagrangian systems also states the following 

Theorem: To every one-parameter group of continuous transformations that 
transform the action function of the system, leaving dynamical equations invariant, 
in the form 

A{5gxi, Sgx2) = Xs) + B{x2)5g - B{xi)5g, 

and where B{x) is a function defined on the kinematical space, there is associated 
a classical observable A^, which is a constant of the motion. 



Let us assume the existence of a r-parameter continuous group of transformations G, of 
the enlarged configuration space (t,qi), that can be extended as a transformation group to 
the whole kinematical space X. Let 6g be an infinitesimal element of G with coordinates 
5(7", a = 1, . . . , r and its action on these variables be given by: 

t~^t' = t + 5t = t + Ma{t,q)5g", (1.12) 
q^it)^q'At') = q,it) + Sq,{t)=q,it) + M^^\t,q)Sg", (1.13) 

and its extension on the remaining kinematical variables by 

qf\t') = q^\t) + 5q^\t) = q'^\t) + M^l\t,q,q^'^)5g^, (1.14) 

and in general 

qf\t') = )(t) + 5qi'\t) = qi'\t) + {t, g, • • • , g^^V?", . = 0, 1, . . . , A; - 1, (1.15) 

where M« and M^^ are functions only of qi and t while the functions m'i^ with s > 1, 
obtained in terms of the derivatives of the previous ones, will be functions of the time t and 
of the variables qi and their time derivatives up to order s. 

For instance, 

^ dq'S') ^ d{q,{t) + M^l^g^) dt_ 
'^^ ^^>- dt' dt df 

but up to first order in 5g 

^ = l-M^{t,q)5g'^, 

and thus 

^(f) = ^iT) + {^-^^^ - .r^M.(.,,)) 5g-. 
and comparing with II1.14|I we get 

M^:!{t^q,q^^y)^'-^^^^~c^pM^[t^ql 

where the total time derivative 

dMl'J {t, q) ^ 9mL"^ {t, q) ^ dMl^'J {t, q) (i) 
dt dt ^ dq, ■ 

(s) (s — \) 

The remaining M^J for s > 1, are obtained in the same way from the previous ' . 
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Under 5g the change of the action functional of the system is: 



5A[q] = f'L{t',qf{t'))dt'- rL{t,q^\t))dt 
Jt[ Jti 

L(t + 5t, q^\t) + 5q^\t))dt' - r L(t, q^\t))dt. 

Ju 



By replacing in the first integral the integration range (t[,t2) by (^1,^2) having in mind the 
Jacobian of t' in terms of t, this implies that the differential dt' = (1 + d{dt)/dt)dt, and thus: 

6A[q] = £ Lit + 5t, g« + 5gf )) {l + ^) dt - Lit, <^^)dt 
f'^ ( ,di5t) dL^ dL ^ A , 

keeping only for the Lagrangian L(t + g*^"*^ + 6q^^^), first order terms in its Taylor expansion. 

Now, in the total variation of 6qf\t) = q''f\t') — qi''\t) is contained a variation in the 
form of the function qi'^\t) and a variation in its argument t, that is also affected by the 
transformation of the group, i.e., 

Sq^'^ = qf {t + 6t) - ) it) = qf (t) - g« (t) + (rfgf ^ it)/dt)6t 

= Sq^\t) + qt'\t)St, 

where 5q^^\t) is the variation in form of the function qi^\t) at the instant of time t. Taking 
into account that for the variation in form 

5q'f\t) = d\5qi{t))/dt^ = d{5qt'\t))/dt, 

it follows that 



*2 f d{L6t) , dL ^js) 



+ ^S€'it) dt. (1.16) 



dt dqf^ 



Making the replacements 



dL- dL- 

oqi oqi 

dL ^ (I) dL diSqi) d ( dL ^ \ d ( dL . 
Mi = . m J. = 31 -^7T^^<1^ - 31 TTTTT 



dqf^ dq^^ dt dt ^Qq^) ) dt \Qqf) 



dqP dt \dqP ) dt ^agf) 

d ( dL - d ( d ( dL\- \ I dL , ^ 

= s ( ) - S I S ( ^ j ) + ( ^ I 
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dL 



d dL 



and collecting terms we get 
6A[q]- 



(fc-i) 
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dL 



did 

dt\dt\ ^qf) 



Sq, 



(fe-2) 



r 9L 

+oqi 



ft2 

dt < 
, 1 


'd{L8t) 
dt 




d 




+ ••■ 


~ Jt 


[dq?) 


' dL 


d f 


dL \ 




~ ~dt\ 





d^ 



dL 



[k) 



4t I 



dL d dL 



+ 



dqP dt {dq^f^ 
dL 



+ 



+ (-1 



.,d^ 
dt''-' \dql 
c/'=-2 / dL 



dL 



ife-2 



df"^ 



Mi 



(k) 



dq 



(k) 



The terms between squared brackets are precisely the conjugate momenta of order s, p^^^-^, 
except the first one, which is the left-hand side of and vanishes identically if the functions 
qi satisfy the dynamical equations. 

In ordinary Lagrangian systems that depend only on first order derivatives of the independent de- 
grees of freedom, the canonical approach associates to every generaHzed coordinate qi a dynamical 
variable pi, called its canonical conjugate momentum and defined by 

dL 



As a generalization of this, in Lagrangian systems with higher order derivatives, a generalized 
canonical formalism can be obtained by defining various canonical conjugate momenta (up to a 
total of k of them) associated to each of the independent degrees of freedom qi : ^ 



k — s 



dq. 



(r+s) 



S — 1 , . . . ^ , Z — 1 , . . . , ?Tj 



(1.17) 



which are precisely the above terms between the squared barckets. It is said that p^^,-) is the 
conjugate momentum of order s of the variable qi. 



Now if we introduce in the integrand the variables qi that satisfy Euler-Lagrange equations, 
the variation of the action functional ()1.16|l is transformed into the variation of the action 
function along the classical trajectory, and therefore, the variation of the action function can 
be written as, 

6A{xu X2) = I {L6t + [dqipl^ + 5qf^p\^^ + ■■■ + 5qt'\\k)) } dt. (1.18) 
^ E.T.Whittaker, Analytical Dynamics, Cambridge University Press, Cambridge (1927), p. 265. 
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with given in p.l7|) . If we replace in (I1.18|l the form variation Sq^'^^ = Sq^'^^ — qi^'^^\ then 
6A{xu X,) = j'^ I {L5t + 5gf V(.+i) - g?^'V(.+i)5t} dt (1.19) 



with the usual addition convention. By substitution of the variations 5t and Sqf'^ in terms of 
the infinitesimal element of the group 5g", ()1. 13111. 15| . we get: 

M(xi, X2) = I |(l - pl^qf^) M„ + p;„+,)M,S:) } 5g^dt, (1.20) 

with the following range for repeated indexes for the addition convention, i = 1, . . . ,n, s = 
1, . . . , /c, M = 0, 1, . . . , /c — 1 and a = 1, . . . , r. 

In the above integral we are using the solution of the dynamical equations, and therefore 
the variation of the action function is 

5A{xi, X2) = A{6gxi, Sgx2) — A{xi,X2). 

If it happens to be of first order in the group parameters in the form 

6A{x,, X2) = B^{x2)6g'' - B^{x,)6g'', (1.21) 

then equating to (ll.2()|l we can perform the trivial time integral on the right hand side. By 
considering that the group parameters 6g" are arbitrary, rearranging terms depending on ti 
and t2 on the left- and right-hand side, respectively, we get several observables that take the 
same values at the two arbitrary times ti and ^2- They are thus constants of the motion and 
represent the time conserved physical quantities, 

iV„ = B^ix) -{l- pl^^qf^ M„ - pJ,+i)M,^:\ « = 1, . . . , r, (1.22) 

(s) 

where the term within brackets H = p^^^^^ql — L is the generalized Hamiltonian. 

These are the r Noether constants of the motion related to the infinitesimal transformations 
()1.21|) of the action function under the corresponding r-parameter Lie group. 

To express the different magnitudes in terms of the kinematical variables, let us define the 
variables x^ according to the rule: = t, = qi, x"^"*"* = g-^\ . . . ^ a;*^'^"^)""'"* = Since 
L = and q^^^ = qf~^^ /x^ , the derivatives in the definition of the canonical momenta can 

be written as: 

dL_ ^ d{L/x^) ^ dL ^ 

in terms of the functions of the expansion (|1.9jl of the Lagrangian. The different conjugate 
momenta appear in the form: 

r=0 

in terms of the functions Fj and their time derivatives. Therefore the Noether constants of the 
motion are written as 

Na = B4x) - (^F,- - pl^ -—^) Mo. - p[s+DM^i- (1-25) 
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We see that the Noether constants of the motion Na are finally expressed in terms of 
the functions Fi and their time derivatives, of the functions Mj;^ which represent the way 
the different kinematical variables transform under infinitesimal transformations, and of the 
functions Ba which, as we shall see below, are related to the exponents of the group G. 
Functions F^ and their time derivatives are homogeneous functions of zero degree in terms of 
the derivatives of the kinematical variables x\ Functions Ba{x) and M^^{x) depend only on the 
kinematical variables. Consequently, Noether constants of the motion are also homogeneous 
functions of zero degree in terms of the derivatives of kinematical variables and thus invariant 
under arbitrary changes of evolution parameter. 

1.4 Lagrangian gauge functions 

In the variational formulation of classical mechanics 

A[q] = [ L{t,q'f\t))dt= [ L{x,x)dT, (1.26) 

Jti Jti 

A[q] is a path functional, i.e., it takes in general different values for the different paths joining 
the fixed end points xi and X2- Then it is necessary that Ldr be a non-exact differential. 
Otherwise, if Ldt = d\, then ^[g] = A2— Ai and the functional does not distinguish between the 
different paths and the action function of the system from Xi to X2, A{xi, X2) = A(x2) — A(xi), 
is expressed in terms of the potential function A(x), and is thus, path independent. 

If A(x) is a real function defined on the kinematical space X of a Lagrangian system with 
action function A{xi,X2), then the function A'{xi,X2) = A{xi,X2) + \{x2) — A(xi) is another 
action function equivalent to A(xi, X2). In fact it gives rise by (jl.ll|l to the Lagrangian L' that 
differs from L in a total r-derivative. ^ 

Using (jl.llll . we have 

L'{x,x)=L{x,x) + ^, (1.27) 
dr 

and therefore L and L' produce the same dynamical equations and A{xi,X2) and A'{xi,X2) 
are termed as equivalent action functions. 

Let G be a transformation group of the enlarged configuration space it,qi), that can be 
extended to a transformation group of the kinematical space X. Let g & G he an arbitrary 
element of G and x' = gx, the transform of x. Consider a mechanical system characterized by 
the action function A{xi,X2) that under the transformation g is changed into A{x[,X2). If G 
is a symmetry group of the system, i.e., the dynamical equations in terms of the variables x' 
are the same as those in terms of the variables x, this implies that A{x[,x'2) and A{xi,X2) are 
necessarily equivalent action functions, and thus they will be related by: 

Ma^i, 9x2) = A{xi, X2) + a{g; X2) - a{g; Xi). (1.28) 

The function a will be in general a continuous function of g and x. This real function a{g; x) 
defined on G x X is called a gauge function of the group G for the kinematical space X. 
Because of the continuity of the group it satisfies a{e;x) = 0, e being the neutral element of 
G. If the transformation g is infinitesimal, let us represent it by the coordinates Sg'' , then 



J.M. Levy-Leblond, Comm. Math. Phys. 12, 64 (1969). 
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a{6g;x) = 5g"B„{x) to first order in the group parameters. The transformation of the action 
function takes the form 

A{5gxi, 6gx2) = X2) + 5g''B„{x2) - 6g'^Ba{xi), 

i.e., in the form required by Noether's theorem to obtain the corresponding conserved quanti- 
ties. In general, B„ functions for gauge-variant Lagrangians are obtained by 



j,^^^^^daig;x) 



dg- 



(1.29) 

9=0 



Because of the associative property of the group law, any gauge function satisfies the 
identity 

«(fl''; gx) + a{g; x) - a{g'g; x) = ^{g, g), (1.30) 

where the function ^, defined on G x G, is independent of x and is an exponent of the group 
G. 

This can be seen by the mentioned associative property of the group law. From lll.28|l we get: 
A{9'9Xi,g'gx2) = A(xi,X2) + a{g'g;x2) - a{g'g;xi), (1.31) 

and also 

Mg' gxi, g' gx2) = A{gxi,gx2) + a{g'; 3x3) - a{g'] gxi) 

= A{xi,X2) + a{g; X2) ~ a(g; xi) + a{g'; gx2) - a{g'; gxi), 

and therefore by identification of this with the above ljl.31|l , when collecting terms with the same 
X argument we get 

"(5'; 32^2) + oi{g; X2) - a{g'g; X2) = gxi) + a{g\ x{) - a{g' g; xi), 

and since xi and X2 are two arbitrary points of X, this expression is l|1.3fl|l and defines a function 
£,{9'ig)y independent of a;. 

If we substitute this function C,{g',g) into (ll.45jl we see that it is satisfied identically. For 
g' = g = e, it reduces to ^(e, e) = a(e; x) = 0, and thus ^ is an exponent of G. 

It is shown by Levy-Leblond in the previous reference that if X is a homogeneous space 
of G, i.e., if there exists a subgroup H of G such that X = G/H, then, the exponent ^ is 
equivalent to zero on the subgroup H, and gauge functions for homogeneous spaces become: 

ai9;x) = ^{g,K), (1.32) 

where hx is any group element of the coset space represented by x G G/H. 

For the Poincare group V all its exponents are equivalent to zero and thus the gauge 
functions when X is a homogeneous space of V are identically zero. Lagrangians of relativistic 
systems whose kinematical spaces are homogeneous spaces of P can be taken strictly invariant. 

However, the Galilei group Q has nontrivial exponents, that are characterized by a param- 
eter m that is interpreted as the total mass of the system, and thus Galilei Lagrangians for 
massive systems are not in general invariant under Q. In the quantum formalism, the Hilbert 
space of states of a massive nonrelativistic system carries a projective unitary representation 
of the Galilei group instead of a true unitary representation. ^ 



* see ref.7 and also J.M. Levy-Leblond, Galilei Group and Galilean Invariance, in E.M. Loebl, Group Theory 
and its applications, Acad. Press, NY (1971), vol. 2, p. 221. 
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1.5 Elementary systems 

In Newtonian mechanics the simplest geometrical object is a point of mass m. Starting 
with massive points we can construct arbitrary systems of any mass and shape, and thus any 
distribution of matter. The massive point can be considered as the elementary particle of 
Newtonian mechanics. In the modern view of particle physics it corresponds to a spinless 
particle. We know that there exist spinning objects like electrons, muons, photons, neutrinos, 
quarks and perhaps many others, that can be considered as elementary particles in the sense 
that they cannot be considered as compound systems of other objects. Even more, we do not 
find in Nature any spinless elementary particles. It is clear that the Newtonian point does 
not give account of the spin structure of particles and the existence of spin is a fundamental 
intrinsic attribute of an elementary particle, which is lacking in Newtonian mechanics, but it 
has to be accounted for. 

In quantum mechanics, Wigner's work ^ on the representations of the inhomogeneous 
Lorentz group provides a very precise mathematical definition of the concept of elementary 
particle. An elementary particle is a quantum mechanical system whose Hilbert space of 
pure states is the representation space of a projective unitary irreducible representation of the 
Poincare group. Irreducible representations of the Poincare group are characterized by two 
invariant parameters m and S, the mass and the spin of the system, respectively. By finding 
the different irreducible representations, we can obtain the quantum description of massless 
and massive particles of any spin. 

The very important expression of the above mathematical definition, with physical conse- 
quences, lies in the term irreducible. Mathematically it means that the Hilbert space is an 
invariant vector space under the group action and that it has no other invariant subspaces. 
But it also means that there are no other states for a single particle than those that can be 
obtained by just taking any arbitrary vector state, form all its possible images in the different 
inertial frames and finally produce the closure of all finite linear combinations of these vectors. 

We see that starting from a single state and by a simple change of inertial observer, we 
obtain the state of the particle described in this new frame. Take the orthogonal part of this 
vector to the previous one and normalize it. Repeat this operation with another kinematical 
transformation acting on the same first state, followed by the corresponding orthonormaliza- 
tion procedure, as many times as necessary to finally obtain a complete orthonormal basis of 
the whole Hilbert space of states. All states in this basis are characterized by the physical 
parameters that define the first state and a countable collection of group transformations of 
the kinematical group G. And this can be done starting from any arbitrary state. 

This idea allows us to define a concept of physical equivalence among states of any arbi- 
trary quantum mechanical system in the following way: Two states are said to be physically 
equivalent if they can produce by the above method an orthonormal basis of the same Hilbert 
subspace, or in an equivalent way, if they belong to the same invariant subspace under the 
group action. It is easy to see that this is an equivalence relation. But if the representation is 
irreducible, all states are equivalent as basic pieces of physical information for describing the 
elementary system. There is one and only one single piece of basic physical information to 
describe an elementary object. That is what the term elementary might mean. 

But this definition of elementary particle is a pure group theoretical one. The only quantum 
mechanical ingredient is that the group operates on a Hilbert space. Then one question arises. 



^ see ref.l. 
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Can we translate this quantum mechanical definition into the classical domain and obtain an 
equivalent group theoretical definition for a classical elementary particle? 

Following with the above idea, in classical mechanics we have no vector space structure 
to describe the states of a system. What we have are manifolds of points where each point 
represents either the configuration state, the kinematical state or the phase state of the system 
depending on which manifold we work. But the idea that any point that represents the state 
of an elementary particle is physically equivalent to any other, is in fact the very mathematical 
concept of homogeneity of the manifold under the corresponding group action. In this way, 
the irreducibility assumption of the quantum mechanical definition is translated into the realm 
of classical mechanics in the concept of homogeneity of the corresponding manifold under the 
Poincare group or any other kinematical group we consider as the symmetry group of the 
theory. But, what manifold? Configuration space? Phase space? The answer as has been 
shown in previous works, is that the appropriate manifold is the kinematical space. 

In the Lagrangian approach of classical mechanics, the kinematical space X is the manifold 
where the dynamics is developed as an input-output formalism. When quantizing the system we 
will obtain the natural link between the classical and quantum formalisms through Feynman's 
path integral approach, as will be shown later. This manifold is the natural space on which 
to define the Hilbert space structure of the quantized system. In a formal way we can say 
that each point x ^ X that represents the kinematical state of a system is spread out and 
is transformed through Feynman's quantization into the particle wave function i'^'(x) defined 
around x. This wave function is a squared integrable complex function defined on X. 

We can also analyze the elementarity condition from a different point of view. Let us 
consider an inertial observer O that is measuring a certain observable A{t) of an arbitrary 
system at an instant r. This observable takes the value A'{t) for a different inertial observer 
O'. It can be expressed in terms of A{t) in the form A'{t) = f{A{T),g), where g is the 
kinematical transformation between both observers. At instant r + dr, the corresponding 
measured values of that observable will have changed but A'{t + dr) = f{A{T + dT),g) with 
the same g as before, and assuming that the evolution parameter r is group invariant. 

But if the system is elementary, we take as an assumption that the modifications of the 
observables produced by the dynamics can always be compensated by a change of inertial 
reference frame. Then, given an observer O, it is always possible to find at instant t + dr 
another inertial observer O' who measures the value of an essential observable A'{t + dr) with 
the same value as O does at instant r, i.e., A'{T + dT) = A{t). If the system is not elementary, 
this will not be possible in general because the external interaction might change its internal 
structure, and thus it will not be possible to compensate the modification of the observable 
by a simple change of inertial observer. Think about a non-relativistic description of an atom 
that goes into some excited state. The new internal energy, which is Galilei invariant, cannot 
be transformed into the old one by a simple change of reference frame. 

But the essential observables are the kinematical variables. From the dynamical point of 
view we can take as initial and final points any Xi and X2 E X , compatible with the causality 
requirements. This means that any x can be considered as the initial point of the variational 
formalism. In this way, at any instant r if the system is elementary, we can find an infinitesimal 
kinematical transformation 6g{T) such that 

x'{t + dr) = f{x{T + dr), 5g{T)) = x(t). 



M. Rivas, J. Phys. A 18, 1971 (1985); J. Math. Phys. 30, 318 (1989); J. Math. Phys. 35, 3380 (1994). 
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or by taking the inverse of this transformation, 

x{T + dT)^f-\xiT),Sg{T)). 

This equation represents the dynamical evolution equation in X space. Knowledge of the 
initial state xi and the function Sg{T) completely determines the evolution of the system. In 
general, Sgi^r) will depend on the instant r, because the change of the observables depends on 
the external interaction. But if the system is elementary and the motion is free, all Sg{T) have 
necessarily to be the same, and thus r independent. We cannot distinguish in a free motion 
one instant from any other. Then, starting from Xi we shall arrive at X2 by the continuous 
action of the same infinitesimal group element Sg, and the free particle motion is the action 
of the one-parameter group generated by Sg on the initial state. Therefore, there should exist 
a finite group element g E G such that X2 = gxi. If the evolution is not free, the composition 
of all infinitesimal group elements Sg{T) for all intermediate values of r e [Ti,r2], will also 
produce a finite group element g, and thus, X2 — gxi. We thus arrive at the: 

Definition: A classical elementary particle is a Lagrangian system whose 
kinematical space X is a homogeneous space of the kinematical group G. 

Usually the Lagrangian of any classical Newtonian system is restricted to depend only on 
the first order derivative of each of the coordinates qi that represent the independent degrees 
of freedom, or equivalently, that the Qi satisfy second order differential equations. But at this 
stage, if we do not know what are the basic variables we need to describe our elementary 
system, how can we state that they necessarily satisfy second order differential equations? If 
some of the degrees of freedom, say qi, q2 and q^, represent the center of mass position of the 
system, Newtonian mechanics implies that in this particular case L will depend on the first 
order derivatives of these three variables. But what about other degrees of freedom? It is 
this condition on the kinematical space to be considered as a homogeneous space of G, as the 
mathematical statement of elementarity, that will restrict the dependence of the Lagrangian 
on these higher order derivatives. It is this definition of elementary particle with the proper 
election of the kinematical group, which will supply information about the structure of the 
Lagrangian. 

The Galilei and Poincare groups are ten-parameter Lie groups and therefore the largest 
homogeneous space we can find for these groups is a ten-dimensional manifold. The variables 
that define the different homogeneous spaces will share the same domains and dimensions 
as the corresponding variables we use to parameterize the group. Both groups, as we shall 
see later, are parameterized in terms of the following variables (6, a, v, oc) with domains and 
dimensions respectively like 6 e M that represents the time parameter of the time translation 
and a e R^, the three spatial coordinates for the space translation. Parameter G are the 
three components of the relative velocity between the inertial observers, restricted to v < c 
in the Poincare case. Finally a G 5*0(3) are three dimensionless variables which characterize 
the relative orientation of the corresponding Cartesian frames and whose compact domain is 
expressed in terms of a suitable parametrization of the rotation group. 

In this way the maximum number of kinematical variables, for a classical elementary parti- 
cle, is also ten. We represent them by a; = (t, r, u, tx) with the same domains and dimensions 
as above and interpret them respectively as the time, position, velocity and orientation 
of the particle. 



1.5. ELEMENTARY SYSTEMS 
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Because the Lagrangian must also depend on the next order derivatives of the kinematical 
variables, we arrive at the conclusion that L must also depend on the acceleration and angular 
velocity of the particle. The particle is a system of six degrees of freedom, three r, represent 
the position of a point and other three a, its orientation in space. We can visualize this by 
assuming a system of three orthogonal unit vectors linked to point r as a body frame. But the 
Lagrangian will depend up to the second time derivative of r, or acceleration of that point, 
and on the first derivative of a, i.e., on the angular velocity. The Galilei and Poincare groups 
lead to generalized Lagrangians depending up to second order derivatives of the position. 

By this definition it is the kinematical group G that implements the special Relativity Prin- 
ciple that completely determines the structure of the kinematical space where the Lagrangians 
that represent classical elementary particles have to be defined. Point particles are particu- 
lar cases of the above definition and their kinematical spaces are just the quotient structures 
between the group G and subgroup of rotations and boosts, and thus their kinematical vari- 
ables reduce only to time and position (t,r). Therefore, the larger the kinematical group of 
space-time transformations, the greater the number of allowed classical variables to describe 
elementary objects with a more detailed and complex structure. In this way, the proposed 
formalism can be accommodated to any symmetry group. It is the proper definition of this 
group which contains the physical information of the elementary particles. 



1.5.1 Elementary Lagrangian systems 

An elementary Lagrangian system will be characterized by the Lagrangian function L[x, x) 
where the variables x e X lie in a homogeneous space X of G. L is a homogeneous function 
of first degree of the derivatives of the kinematical variables, and this allows us to write 



L{x, x) = Fi{x, x) x^ 



(1.33) 



Functions Fi{x,x) are therefore homogeneous functions of zero degree in the variables and 
summation convention on repeated indexes as usual is assumed. 

Under G, x transforms as x' = gx or more explicitly its coordinates by x'^ = f^{g,x), and 
their derivative variables 

dx 



X 



dx^ 



x\ 



transform like the components of a contravariant vector. 
The Lagrangian transforms under G, 

L{x'{x),x'{x,x)) — L{x,x) + 



da{g; x) 
dr 



(1.34) 



(1.35) 



I.e., 



Fi{x',x')x' = Fj{x,x)x^ + ^^^fi-i'- 
Taking into account the way the different variables transform, we thus arrive at: 



(1.36) 



Fi{x',x') 



dx'-> 
dx"" 



Fj{x, x) + 



da{g; x) 
dx^ 



see ref.l4. 



(1.37) 
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In the case when a{g;x) = 0, they transform like the components of a covariant vector over 
the kinematical space X. But in general this will not be the case and a{g;x) contains basic 
physical information about the system. 

We thus find that for a fixed kinematical space X, the knowledge of the group action of G 
on X, and the gauge function a{g;x), will give us information about the possible structure of 
the functions Fi{x,x), and therefore about the structure of the Lagrangian. 

In practice, if we restrict ourselves to the Galilei Q and Poincare V groups, we see that V 
has gauge functions equivalent to zero and thus Poincare Lagrangians that describe elementary 
particles can be taken strictly invariant. In the case of the Galilei group, it has only one class 
of gauge functions that define the mass of the system, and thus nonrelativistic Lagrangians 
will be in general not invariant. In the particular case of Galilei invariant Lagrangians, they 
will describe massless systems. 



1.6 Appendix: Lie groups of transformations 

Let us introduce the notation and general features of the action of Lie groups on continuous 
manifolds to analyze the transformation properties of the different magnitudes we can work 
with in either classical or quantum mechanics. We shall use these features all throughout this 
book. 

Let us consider the transformation of an n-dimensional manifold X , x' — gx given by n 
continuous and differentiable functions depending on a set e G of r continuous parameters 
of the form 

x"^f(x^;g''), yxeX, Vg E G, i,j^l,...,n, a = l,...,r. 

This transformation is said to be the action of a Lie group of transformations if it fulfils the 
two conditions: 

(i) G is a Lie group, i.e., there exists a group composition law c = (p{ci,b) e G, \/a,b e G, in 
terms of r continuous and differentiable functions (j)'^. 

(ii) The transformation equations satisfy 

x" = fix'; b) = fifix; a); 6) = f{x; c) = f{x; 4>{a, b)). 

The group parametrization can be chosen such that the coordinates that characterize the 
neutral element e of the group are e = (0, . . . , 0), so that an infinitesimal element of the group 
is the one with infinitesimal coordinates 8g" , a = 1, . . . , r. 

Under the action of an infinitesimal element bg of the group G, the change in the coordinates 
x^ of a point a; G X is given by 



x' + dx' = f{x; 5g) = + ■ 



dg^ 



5g% 



g=e 



after a Taylor expansion up to first order in the group parameters and with = /'(x;0). 
There are nr auxiliary functions of the group that are defined as 



. df{x;g) 
u^{x) - 



dg'^ 



(1.38) 

g=e 



and therefore to first order in the group parameters, dx^ = u^^{x)5g" . 
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The group action on the manifold X can be extended to the action on the set J^{X) of 
continuous and differentiable functions defined on X by means of: 

g : h{x) h'{x) = h{gx). (1.39) 

If the group element is infinitesimal, then 

h\x) = h{x' + dx') = h{x' + ui{x)6g'') = h{x) + ^(s)^^", 

after a Taylor expansion to first order in the infinitesimal group parameters. The infinitesimal 
transformation on J-'{X) can be represented by the action of a differential operator in the form 

h'{x) = (l + 5g''ul{x)4-] h{x) = iI + 5g''X,)h{x) = U{5g)h{x), 



dx ^ 

where I is the identity operator and the linear differential operators 

X. = <(x)A. (1.40) 

In particular, when acting with the operator U {6g) = (I + 5g"X^) on the coordinate x^ we get 
x^ + dx^ = x^ + u^^{x)6g'^ . 

The operators X^ are called the generators of the infinitesimal transformations. They 
are r linearly independent operators that span an r-dimensional real vector space such that its 
commutator also belongs to the same vector space, i.e., 

[X., Xa] = C X„, a, a, A = 1, ... , r. (1.41) 

The coefficients d^^ ^ ^^t of real constant numbers, called the structure constants of 
the group, and the vector space spanned by the generators is named the Lie algebra C{G), 
associated to the Lie group G. The structure constants are antisymmetric in their lower indexes 
c'^^ = — c"^, and satisfy Jacobi's indentitites: 

Cxia + cl,c^.a + clA. = 0, Va, A, /i, /? = 1, . . . , r. 

Equations (11.4111 are the commutation relations that characterize the structure of the Lie 
algebra of the group. 

If a finite group transformation of parameters g^ can be done in n smaller steps of param- 
eters g" /n, with n sufficiently large, then a finite transformation U{g)h{x) can be obtained 

as ^ 

U{g)h{x) = lim (I + —X^] h{x) = exp(c/"X,) h{x). 

n-»oo \ n J 

This defines the exponential mapping and in this case the group parameters g" are called 
normal or canonical parameters. In the normal parameterization the composition law of one- 
parameter subgroups reduces to the addition of the corresponding parameters of the involved 
group elements. 

Consider J-'{X) a Hilbert space of states of a quantum system; ()1.39|1 can be interpreted 
as the transformed wave function under the group element g. Then if the operator U{g) is 
unitary it is usually written in the explicit form 

U{g) = exp g'^X^ 
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in terms of the imaginary unit i and Planck's constant h, such that in this case the new X„ 
above are self-adjoint operators and therefore represent certain observables of the system. The 
physical dimensions of these observables depend on the dimensions of the group parameters g'^, 
since the argument of the exponential function is dimensionless and because of the introduction 
of Planck's constant this implies that g'^Xa- has dimensions of action. These observables, taking 
into account ()1.4n|) . are represented in a unitary representation by the differential operators 



However, (|1.39p is not the most general form of transformation of the wave function of a 
quantum system, as we shall see in Chapter 3, but once we know the way it transforms we 
shall be able to obtain the explicit expression of the group generators by a similar procedure 
as the one developed so far. In general the wave function transforms under continuous groups 
with what is called a projective unitary representation of the group, which involves in general 
some additional phase factors. 

1.6.1 Casimir operators 

When we have a representation of a Lie group either by linear operators or by matrices 
acting on a linear space, we can define there what are called the Casimir operators. They are 
operators C that can be expressed as functions of the generators of the Lie algebra with the 
property that they commute with all of them, i.e., they satisfy [C, Xg.] = 0, Vo" = 1, . . . ,r. 
In general they are not expressed as real linear combinations of the X^ and therefore they do 
not belong to the Lie algebra of the group. They belong to what is called the group algebra, 
i.e., the associative, but in general non-commutative algebra, spanned by the real or complex 
linear combinations of products of the X^-, in the corresponding group representation. 

In those representations where the X^- are represented by self-adjoint operators as in a 
quantum formalism, the Casimir operators may be also self-adjoint and will represent those 
observables that remain invariant under the group transformations. In particular, when we 
consider later the kinematical groups that relate the space-time measurements between in- 
ertial observers, the Casimir operators of these groups will represent the intrinsic properties 
of the system. They are those properties of the physical system whose measured values are 
independent of the inertial observers. 

For semisimple groups, i.e., for groups that do not have Abelian invariant subgroups like the 
rotation group SO (3), the unitary groups SU{n) and many others, it is shown that the Casimir 
operators are real homogeneous polynomials of the generators X^-, but this is no longer the 
case for general Lie groups. Nevertheless, for most of the interesting Lie groups in physics, like 
Galilei, Poincare, De Sitter, S'L(4, M), the inhomogeneous JS'L(4, M) and Conformal SU{2,2) 
groups, the Casimir operators can be taken as real polynomial functions of the generators. 

1.6.2 Exponents of a group 

The concept of exponent of a continuous group G was developed by Bargmann in his work 
on the projective unitary representations of continuous groups. 



X^ = -ul{x) 



_d_ 



(1.42) 



12 V. Bargmann, Ann. Math. 59, 1 (1954). 
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Wigner's theorem about the symmetries of a physical system is well known in Quantum 
Mechanics. 

It states that if 7i is a Hilbert space that characterizes the pure quantum states of a system, 
and the system has a symmetry S, then there exists a unitary or antiunitary operator U{S), 
defined up to a phase, that implements that symmetry on 7i, i.e., if and ip E Ti. are two 
possible vector states of the system and | < > p is the transition probability between 
them and U{S)(f) and U{S)ip represent the transformed states under the operation S, then 

I < U{S)<p\U{S)ij > |2 = I < 01^ > |2. 

If the system has a whole group of symmetry operations G, then to each element g E G 
there is associated an operator U{g) unitary or antiunitary, but if G is a continuous group, in 
that case U{g) is necessarily unitary. This can be seen by the fact that the product of two 
antiunitary operators is a unitary one. 

Because there is an ambiguity in the election of the phase of the unitary operator U{g), 
it implies that in general U{gi)U{g2) 7^ U{gig2) and therefore the transformation of the wave 
function is not given by an expression of the form (ll.39jl , but it also involves in general a phase 
factor. However in the case of continuous groups we can properly choose the corresponding 
phases of all elements in such a way that 

U{gi)U{g2) = uj{gug2)U{gig2), (1.43) 

where uj{gi, g2) = exp{i^{gi, g2)} is a phase that is a continuous function of its arguments. The 
real continuous function on G x G, ^{gi, (72) is called an exponent of G. The operators U{g) 
do not reproduce the composition law of the group G and (ll.43jl represents what Bargmann 
calls a projective representation of the group. 

If we use the associative property of the group law, we get 

iU{g,)Uig2))U{gs) = cuig,, g2)U{g^g2)U{gs) 

= ^{.gi,g2)uj{,gig2,gz)U{gig2gz), 

and also 

U{gi){U{g2)U{g,)) = U{g,)uj{g2.g,)U{g2g,) 

= ^{.gi^g2g-i)^{.g2,gz)U{gig2gz)- 

Therefore 

^{gi^ g2)^{gig2, gs) = ^{gi, g2gz)^{g2, gs), (i-44) 

which in terms of the exponents becomes: 

^(^1, ^2) + ^{gm, gs) = ^{gu ^2^3) + ^(^2, gs)- (1-45) 

Because of the continuity of the exponents, 

^{g,e) = ae,g) = 0, Vg E G, (1.46) 

where e is the neutral element of the group. 

E.P. Wigner, Group theory and its application to the quantum mechanics of atomic spectra, Acad. Press, 
NY (1959); V. Bargmann, J. Math. Phys. 5, 862 (1964). 
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Any continuous function on G, (j){g), with the condition 0(e) = 0, can generate a trivial 
exponent by 

that satisfies (jl.45|) and (II .461) . All trivial exponents are equivalent to zero exponents, and 
in a unitary representation ()1.43|) can be compensated into the phases of the factors, thus 
transforming the projective representation ()1.43p into a true unitary one. 

Given a continuous group, the existence or not of non-trivial exponents is an intrinsic group 
property related to the existence or not of central extensions of the group. 



1.6.3 Homogeneous space of a group 

A manifold X is called a homogeneous space of a group G, if Vxi, X2 G X there exists at 
least one element g & G such that X2 = gxi. In that case it is said that G acts on X in a 
transitive way. The term homogeneous reminds us that the local properties of the manifold at 
a point X are translated to any other point of the manifold by means of the group action, and 
therefore all points of X share the same local properties. 

The orbit of a point x is the set of points of the form gx, Wg G G, such that if X is a 
homogeneous space of G, then the whole X is the orbit of any of its points. 

Given a point xq G X, the stabilizer group (little group) of xq is the subgroup H^cg of G, 
that leaves invariant the point xq, i.e., Wh G H^q, hxo = xq- 

If if is a subgroup of G, then every element g E G can be written as 5^ = g'h, where h E H, 
and g' is an element of G/H, the set of left cosets generated by the subgroup H. If X is a 
homogeneous space of G, it can be generated by the action of G on an arbitrary point xq G X. 
Then Vx G X, x = gxo = g'hxo = g'xo, and thus the homogeneous space X is isomorphic to 
the manifold G/H^q- 

The homogeneous spaces of a group can be constructed as quotient manifolds of the group 
by all its possible continuous subgroups. Conversely, it can also be shown that if X a homo- 
geneous space of a group G, then there exists a subgroup H of G such that X is isomorphic 
to G/H. Therefore, the largest homogeneous space of a group is the group itself. 



see ref.7 and J.M. Levy-Leblond, Comm. Math. Phys., 12, 64 (1969); A. A. Kirillov, Elements dela theorie 
des representations, Mir, Moscow (1974). 
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Soluble examples of spinning particles 



Nonrelativistic particles 

2.1 Nonrelativistic point particle 

See the Appendix about the Galilei group at the end of this chapter for the notation used 
through this chapter. 

Let us consider a mechanical system whose kinematical space is the manifold X = Q/TC, 
where H is the six-dimensional subgroup of the homogeneous Galilei transformations of ele- 
ments of the form {0,0,v, fi). See the Appendix at the end of this chapter for the notation 
related to the Galilei group. Then X is a four-dimensional manifold spanned by the variables 
{t, r) = X, with domains t G M, r G M^, similar to the group parameters b and a respectively. 
We assume that they are functions of some evolution parameter r and at any instant r of the 
evolution two different inertial observers relate their measurements by: 

t'ir) = t(r)+6, (2.1) 
r'(r) = R{n)r{T) + vt{T) + a. (2.2) 

Because of the way they transform, we can interpret them respectively as the time and position 
of the system. If we assume that the evolution parameter r is group invariant, by taking the 
r— derivative of both sides of the above expressions, it turns out that the derivatives of the 
kinematical variables at any instant r transform as: 

i'{r) = i{r\ (2.3) 
r'(r) = R{^l)r{T) +vi{T). (2.4) 

There are no constraints among these variables. It is only the homogeneity of the Lagrangian 
in terms of their derivatives (jl.9j) which reduces to three the number of independent degrees 
of freedom. This homogeneity leads to the general form: 

L = Ti + R-r, (2.5) 

where T = dL/di and Ri = dL/drt are still some unknown functions of the kinematical vari- 
ables and their derivatives, which are homogeneous of zero degree in terms of the derivatives. 
Associated to this manifold X, the gauge function for this system is 

"(^?; = ^9, x) = m {vH/2 + V ■ R{fi)r) , (2.6) 
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where the parameter m is interpreted as the mass of the system and ^{g,g') is the exponent 
of Q, so that the transformation of the Lagrangian under the Galilei group is 

L{x', x') = L{x, x)+m {vH/2 + v ■ R{^i)r) . (2.7) 

Then 

dL' fdL 1 ^\ di f dL ^ \ <9ri , , 

dt' \dt 2 J dV \dri ^ ^^'^ Jdt'' ^ ' 

but from (j2.3jl and ()2.4I1 we get di/di' = 1 and dri/di' = —R^^{iJ,)ikVk, respectively, and thus 

T' = T -^mv^ - V ■ R{n)R. (2.9) 

Similarly 

E! = R{ti)R + mv. (2.10) 

The conjugate momenta of the independent degrees of freedom = r,, are pi = dL/drt, and 
consequently Noether's theorem leads to the following constants of the motion: 

a) Under time translations the gauge function p.6|) vanishes, 5t = 56, M = 1, while Sr^ = 
and the constant reduces to the following expression R- dr/dt — L/t = —T. 

b) Under space translations also a{g;x) = 0, (5t = 0, M = 0, while 6ri = 6ai, Mij = 6ij and 
the conserved observable is R. 

c) Under pure Galilei transformations 6t = 6b and M = 0, while 5rj = t6vi and Mij = t6ij, 
but now the gauge function to first order in the velocity parameters is a{6v; x) = mr ■ 6v, and 
we get mr — Rt. 

d) Under rotations a{g; x) = 0, 6t = and M = 0, while 6ri = —eij^rjirikda and Mj^ = —Sijkfj 
the conserved quantity is r x i?. 

Collecting all terms we can give them the following names: 

Energy H = -T, (2.11) 

linear momentum P = R = p, (2-12) 

kinematical momentum K = mr — Pt, (2.13) 

angular momentum J = r x P. (2-14) 



We reserve for these observables the same symbols as the corresponding group generators 
which produce the space-time transformations that leave dynamical equations invariant. Even 
their names make reference to the corresponding group transformation parameter, except the 
energy which in this context should be called the 'temporal momentum'. For the kinemat- 
ical momentum we can find in the literature alternative names like 'Galilei momentum' or 
'static momentum'. Being consistent with this notation, we should call it 'Poincare or Lorentz 
momentum' in a relativistic approach. Nevertheless we shall use the name of kinematical mo- 
mentum for this observable K in either relativistic or non-relativistic formalism. None of these 
conserved quantities associated to the symmetry under one-parameter subgroups are definite 
positive so that the observable H can take any sign. 

The linear momentum takes the general expression P = mr /t = mu because taking the 
r-derivative in (I2.13|l of the kinematical momentum, K = 0, implies P = mu, where u is the 
time derivative of the position of the system, i.e., the velocity of the particle. 



2.2. GALILEI FREE SPINNING PARTICLE 
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The six conditions P = and K = 0, imply u = and r = 0, so that the system is 
at rest and placed at the origin of the observer's reference frame. There is still an arbitrary 
rotation and a time translation to fix a unique inertial observer. Nevertheless we call this class 
of observers, for which P = and K = 0, the center of mass observer. These six conditions 
will be also used as the definition of the center of mass observer for any other system even in 
a relativistic approach. 

From ()2.9p and ()2.in|l we see that the energy and linear momentum transform as: 

H' = H + v ■ R{^i)P +^mv\ (2.15) 
P' = R{^i)P + mv. (2.16) 

Then, if Hq and P = are the energy and linear momentum measured by the center of mass 
observer, for any arbitrary observer who sees the particle moving with velocity n, it follows 
from (pn and (l2l6l) that 

H = Ho + ]^mu^ = Ho + P'^ /2m, P = mu. 
The Lagrangian for the point particle is thus 

• 2 

Tfi V 

L = Tt + R-r = -Ht + P-r = -Hot + —^, (2.17) 

^ t 

with Hq an arbitrary constant which plays no role in the dynamics and can be taken Hq = 0. 
It will be related to the mc^ term of the relativistic point particle. 
If we define the spin of the system, as in (I2.171|) . by 

S = J - —K X P = J -r X P = 0, (2.18) 

m 

it represents the angular momentum of the system with respect to the center of mass r. It 
vanishes, so that the point particle is a spinless system. 



2.2 Galilei free spinning particle 

The most general nonrelativistic particle ^ is the system whose kinematical space X is 
the whole Galilei group Q. Then the kinematical variables are the ten real variables x(r) = 
{t{T),r{T),u{T), p{t)) with domains t G M, r G M^, -u G and p G similarly as the 
corresponding group parameters. The relationship between the values x'{t) and x(r) they 
take at any instant r for two arbitrary inertial observers, is given by: 

t'ir) = t(r)+6, (2.19) 
r'(r) = R{fi)r{T)+vt{T) + a, (2.20) 
u'{t) = R{^i)u{t)+v, (2.21) 



1 M. Rivas, J. Phys. A 18, 1971 (1985). 
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Among these kinematical variables there exist the differential constraints u{t) = r{T)/i{T)^ 
that together with the homogeneity condition of the Lagrangian L in terms of the derivatives 
of the kinematical variables: 

L{x,x) = {dL/dxi)Xi, (2.23) 

reduce from ten to six the essential degrees of freedom of the system. 

These degrees of freedom are the position r{t) and the orientation pit). The Lagrangian 
depends on the second derivative of r{t) and the first derivative of pit). Expression (I2.23jl is 
explicitly given by: 

L = Ti + R-r + U-u + V-p, (2.24) 

where the functions T = dL/di, Ri = dL/dr\ Ui = dL/du\ Vi = dL/dp" will be in general 
functions of the ten kinematical variables (t, r, n, p) and homogeneous functions of zero degree 
in terms of the derivatives {t,r,u, p). By assuming that the evolution parameter r is group 
invariant, these derivatives transform under Q: 

i\r) = i{T), (2.25) 

r'(r) = Rip)rir) +viiT), (2.26) 

ii'ir) = R{p,)u{T), (2.27) 
X (p(^) + /X X p{t)){1 - ■ pjr)) 

= (I-P-Pir))^ + 

/x-p(r)(/x + p(r) + /X X pjr)) 

(I-M-P(r))^ • ^'-''^ 

Instead of the derivative p(r) that transforms in a complicated way, we can define the 
angular velocity of the particle a; as a function of it in the form 

u; = -^{p + pxp). (2.29) 
1 + p"^ 

It is a linear function of p, and transforms as: 

Lj'ir) = R{h)u;{t). (2.30) 

We interpret the rotation matrix R{p) as the rotation that carries the initial frame linked to 
the body at instant r = to the frame at instant r, as in a rigid body. Then, the three 
columns of matrix R{p) represent the Cartesian components of the three unit vectors linked 
to the body when chosen parallel to the laboratory frame at instant r = 0. 

If fc(T) is any internal vector of a rigid body with origin at point r, then its dynamics is contained 
in the expression fc(r) = _R(p(r))fc(0). The velocity of point fe is 

fc(r) = Ripir))m = 7?(p(r))i?-i(p(r))fe(r) = r!(r)fc(r) 

where matrix Q — RR^^ — RR^ is an antisymmetric matrix. At any instant r, R{p{t))R^ [p{t)) = 
I, where superscript T means the transposed matrix and I is the 3x3 unit matrix. Taking the r- 
derivative of this expression, RR^ + RR^ — Q + fl'^ — 0, and thus the three essential components 
of the antisymmetric matrix fl define a three- vector lj 

/ -cj, 

I Lu, 

such that we can also write fc(T) = fl(T)k{T) = u}{t) x A;(t) and a; is interpreted as the instan- 
taneous angular velocity. The different components of uj, expressed as functions of the variables 
p and p are given in II2.29|I . 
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Expression ()2.22|1 corresponds to R{p'{r)) = _R(/x)_R(p(r)). Therefore 

n' = R{p'{T))R^{p'{r)) = i?(M)i?(p(r))i?^(p(r))i?^(A.) 
= R{p)QR-\fi), 



and this leads to the equation ()2.3n|) in terms of the essential components uj of the antisym- 
metric matrix ^l. 

In this way the last part of the Lagrangian {dL/dff)f>'^ can be writen as 

dL ■ dL duj^ 

due to the linearity of uj in terms of p and where Wi = dL/du\ Thus the most general form 



of the Lagrangian of a nonrelativistic particle can also be written instead of (12.2411 as: 

L = Ti + R-r + U- u + W-uj. (2.32) 

Since X is the whole Galilei group Q the most general gauge function is just the group 
exponent: 

a{g; x) = ^{g, K) = m{vh{r)/2 + v ■ R{p)r{r)), (2.33) 

similar to ()2.6|] . and this allows us to interpret the parameter m as the mass of the system. 
Under the action of an arbitrary element of the Galilei group, the Lagrangian L transforms 
according to: 

L{gx{T) , d^gxij)) / (It) = L{x{t),x{t)) + da{g;x{T))/dT. (2.34) 



This leads through some straightforward calculations, similar to the ones performed in (|2T 
(j2.1()j) . to the following form of transformation of the functions: 

T'{t) = T{t) - V R{p)R{t) - mvy2, (2.35) 

R\r) = R{p)R{t) +mv, (2.36) 

U'{t) = R{p)U{t), (2.37) 

W'{r) = R{p)WiT). (2.38) 



2.2.1 Noether constants of the motion 

Using the action of the Galilei group on the kinematical space given by (|2.19|] - (j2.22|) . 
Noether's theorem defines the following constants of the motion: 

a) Under time translation the action function is invariant and as usual we call the corresponding 
conserved quantity, the total energy of the system H. Since 5t = 6b and 5gf ^ = 0, M = 1 
and M^'^^ = 0, by applying ()1.25j) we have: 

H = -{L- p\,)qf^)M = -{L/i - )) = -T - R ■ u - U ■ ii/i - W ■ 

+ {R - dU/dt) - u + U ■ ii/i + V ■ p/i, 
and since W ■ iv = V ■ p, it turns out that 



(2.39) 
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b) Under spatial translation, A{xi,X2) is invariant and this defines the total linear momen- 
tum of the system. We have now: 

6t = 0, M = 0, 6ri = 6ai, M^f = 6ij, 6ui = 0, M^}^ = 0, 
5p, = 0, Mif = 0, 

and then 

P = R-^. (2.40) 

dt ^ ' 

c) Under a pure Galilei transformation of velocity 6v, A{xi, X2) is no longer invariant but taking 
into account (jl.28j) and the gauge function (I2.33jl . it transforms as 6 A = mr2 ■ Sv — mvi ■ 6v 
and this defines the total kinematical momentum K, in the following way: 

St = 0, M = 0, 6ri = 6vit, M^f = Syt, 6ui = 5vi, M^f = Sij, 



and thus 



6p, = 0, M^f 



K = mr- Pt-U. (2.411 



From 1^ = 0, this leads to P = mu — dU/dt, and thus by identification with (j2.4njl . the 
function R = mu irrespective of the particular Lagrangian. The total linear momentum does 
not lie along the velocity of point r. 

d) Finally, under rotations A{xi,X2) remains invariant and the corresponding constant of 
the motion, the total angular momentum of the system, comes from the infinitesimal 
transformation of value = Sai/2, i.e., half of the rotated infinitesimal angle, and then 

St = 0, Mi = 0, Svi = eikjSajVk, M^f = eikjr^, 

Sui = eikjSajUk, M^j^ = eikjUk, 
5 Pi = Sa^iSij + eikjp^ + APi)/2, M'^f = (Sij + ei^jp^ + PiPj)/2, 

which leads to 

' doj^ dpi 

and therefore 

J = rxP + uxU + W = rxP + Z. (2.42) 

We are tempted to consider Z as the spin of the system. Since J = 0, this function Z satisfies 
dZ /dt = Pxu and is not a constant of the motion for a free particle. It is the classical angular 
momentum equivalent to Dirac's spin operator in the quantum case. Because J is the angular 
momentum of the particle with respect to the origin of the observer frame, Z represents the 
angular momentum of the particle with respect to the point r. We shall define the spin for a 
free particle as the angular momentum with respect to its center of mass, once we accurately 
identify the center of mass of the particle. 

The center of mass observer is defined as that inertial observer for whom P = and K = 0. 
These six conditions do not define uniquely an inertial observer but rather a class of them up to 
a rotation and an arbitrary time translation. In fact, the condition P = establishes the class 
of observers for which the center of mass is at rest, and 1^ = is the additional condition to 
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locate it at the origin of coordinates. This comes from the analysis of ()2.4H] . where k = U/m 
is an observable with dimensions of length, and taking the derivative with respect to r of both 
sides, taking into account that P = 0, we have: 

d(r — k) 

K = = mr- Pi-mk, i.e., P = m-— -. (2.43) 

dt ^ ' 

Then the point qr = r — fc is moving at constant speed and we say that it represents the 
position of the center of mass of the system. Thus, the observable = r — q is just the relative 
position of point r with respect to the center of mass. Therefore P = and = give rise 
to dq/dt = 0, and r = k, i.e., g = 0, as we pointed out. With this definition, the kinematical 
momentum can be written as K = mq — Pt, in terms of the center of mass position q and the 
total linear momentum P. 

The spin of the system is defined as the difference between the total angular momentum J 
and the orbital angular momentum of the center of mass motion q x P^ and thus 

1 dk 
S=J~qxP = J K xP=Z + kxP = -mk x — + W. (2.44) 

m dt 

The spin expressed in terms of the constants of the motion J, K and P, is also a constant 
of the motion. 

It is the sum of two terms, one Z = uxU + coming from the new degrees of freedom 
and another k x P, which is the angular momentum of the linear momentum located at point 
r with respect to the center of mass. Alternatively we can describe the spin according to the 
last expression in which the term —k x mdk/dt suggests a contribution of (anti)orbital type 
coming from the motion around the center of mass. It is related to the zitterbewegung or 
more precisely to the function U = mk which reflects the dependence of the Lagrangian on 
the acceleration. The other term W comes from the dependence on the other three degrees of 
freedom pi, and thus on the angular velocity. This zitterbewegung is the motion of the center 
of charge around the center of mass. Point r, as representing the position of the center of 
charge, has been also suggested in previous works for the relativistic electron. ^ 

Because J = 0, and that dW/dr = uj x W and the expression of P, (j2.4()|] . this implies 
the general relation for a free particle 

rxR + uxU + uxW = 0, (2.45) 

which reflects the fact that velocity, acceleration and angular velocity are not independent 
magnitudes, and taking into account that R and r have the same direction, it reduces to 

uxU + ujxW = 0. (2.46) 



2.2.2 Galilei spinning particle of (anti) orbital spin 

To analyze the spin structure of the particle, and therefore the different contributions to 
spin coming from these functions U and W, let us consider the following simpler example. 

Consider a Galilei particle whose kinematical space is X = Q/S0{3), so that any point 
X e X can be characterized by the seven variables x = {t,r,u), u = dr/dt, which are 
interpreted as time, position and velocity of the particle respectively. In this example we have 

2 A.O. Barut and A.J. Bracken, Phys. Rev. D 23, 2454 (1981). 
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no orientation variables. The Lagrangian will also depend on the next order derivatives, i.e., 
on the velocity which is already considered as a kinematical variable and on the acceleration 
of the particle. Rotation and translation invariance implies that L will be a function of only 
n^, {du/dtY and u ■ du/dt = d{u^ /2) / dt, but this last term is a total time derivative and it 
will not be considered here. 

Since from condition (I2.46jl U" ~ -it, let us assume that our elementary system is represented 
by the following Lagrangian, which when written in terms of the three degrees of freedom and 
their derivatives is expressed as 

Parameter m is the mass of the particle because the first term is gauge variant in terms of the 
gauge function (I2.33|l defined by this constant m, while parameter u of dimensions of time~^ 
represents an internal frequency. It is the frequency of the internal zitterbewegung. 

In terms of the kinematical variables and their derivatives, and in terms of some group 
invariant evolution parameter r, the Lagrangian can also be written as 

mir-'^ m V? , , 

L = : 2.48 

2 t 2^2 t ' ^ ' 

where the dot means r-derivative. If we consider that the evolution parameter is dimensionless, 
all terms in the Lagrangian have dimensions of action. Because the Lagrangian is a homoge- 
neous function of first degree in terms of the derivatives of the kinematical variables, L can 
also be written as 

L = Ti + R-r + U-u, (2.49) 

where the functions accompanying the derivatives of the kinematical variables are defined and 
explicitly given by 



dL 


m f dr 


'dt 


~ ~ 2" 


dL 
dr 


dr 


dL 


m d'^r 


dii 


dt^ • 



m 



(<PrV 



= - = (2.60) 
V - l^^-^^S. (2.51) 



Dynamical equations obtained from Lagrangian ()2.47|1 are: 



1 d'^r d'^r 

+ ^ = 0, (2.52) 



u;2 df^ dt^ 
whose general solution is: 

r{t) = A + Bt + Ccosujt + Dsinujt, (2.53) 

in terms of the 12 integration constants A, C and D. 

When applying Noether's theorem to the invariance of dynamical equations under the 
Galilei group, the corresponding constants of the motion can be written in terms of the above 



2.2. GALILEI FREE SPINNING PARTICLE 



37 



functions in the form: 

Energy H = -T-u-^, (2.54) 

^ ^ dU dU , , 

linear momentum P = R — = mu — , (2.55) 

dt dt ^ ' 

kinematical momentum K = mr — Pt — U, (2.56) 
angular momentum J = r x P + u x U. (2-57) 

It is the presence of the U function that distinguishes the features of this system with respect 
to the point particle case. We find that the total linear momentum is not lying along the 
direction of the velocity u, and the spin structure is directly related to the dependence of the 
Lagrangian on the acceleration. 

If we substitute the general solution ()2.53|) in (|2. 54112. 57| we see in fact that the integration 
constants are related to the above conserved quantities 

H = -B'-^{C' + D'), (2.58) 

P = mB, (2.59) 
K = mA, (2.60) 
J = Ax mB -muC X D. (2.61) 



We see that the kinematical momentum K in ()2.56p differs from the point particle case 
()2.13|) in the term —U, such that if we define the vector k = U/m, with dimensions of length, 
then K = leads from (I2.56jl to the equation: 

d(r - k) 
P = m^- -, 

dt ' 

and q = r — k, defines the position of the center of mass of the particle that is a different point 
than r and using (12.5111 is given by 

1 I d^r , , 

q = r--U = r + — —. 2.62 

m 00^ dV 

In terms of it, dynamical equations ()2.52jl can be separated into the form: 

d^q 



= 0, (2.63) 
+ uj\r-q) = 0, (2.64) 



c/t2 

d'^r , 

df^ 

where (12:631) is just eq. (ICTl) after twice differentiating (ITHll . and Equation (IIIM is (ITOll 

after collecting all terms on the left hand side. 

From (|2.63|) we see that point q moves in a straight trajectory at constant velocity while 
the motion of point r, given in ()2.64|1 . is an isotropic harmonic motion of angular frequency oo 
around point q. 

The spin of the system S is defined as 

S = J -qx P = J - —K X P, (2.65) 

m 
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and since it is written in terms of constants of the motion it is clearly a constant of the motion, 
and its magnitude S"^ is also a Galilei invariant quantity that characterizes the system. In 
terms of the integration constants it is expressed as 

S = -mu}CxD. (2.66) 

From its definition we get 

d dk 
S = uxU + kxP = —mir — q) x — (r — q) = —k x m— , (2.67) 

^ dt^ ^ dt ^ ' 

which appears as the (anti) orbital angular momentum of the relative motion of point r around 
the center of mass position q at rest, so that the total angular momentum can be written as 

J = qxP + S = L + S. (2.68) 

It is the sum of the orbital angular momentum L associated to the motion of the center of 
mass and the spin part S. For a free particle both L and S are separately constants of the 
motion. We use the term (anti)orbital to suggest that if vector k represents the position of 
a point mass m, the angular momentum of this motion is in the opposite direction as the 
obtained spin observable. But as we shall see in a moment, vector k does not represent the 
position of the mass m but rather the position of the charge e of the particle. 

2.2.3 Interacting with an external electromagnetic field 

But if q represents the center of mass position, then what position does point r represent? 
Point r represents the position of the charge of the particle. This can be seen by considering 
some interaction with an external field. The homogeneity condition of the Lagrangian in terms 
of the derivatives of the kinematical variables leads us to consider an interaction term of the 
form 

Li = -e(l){t,r)i + eA{t,r) ■ r, (2.69) 

which is linear in the derivatives of the kinematical variables t and r and where the external 
potentials are only functions of t and r. We can also consider more general interaction terms of 
the form N{t, r, u)-u^ and also more general terms in which functions and A also depend on 
u and ii. If the interaction Lagrangian depends on ii this implies that the interaction modifies 
the definition of the observable U = mk which defines the spin of the free system. But if the 
system is elementary the spin definition cannot be changed, so that (I2.69jl is the most general 
interaction term. 

Dynamical equations obtained from L + Lj are 

+ — = —{E{t,r) + ux B{t,r)), (2.70) 



dt^ dt^ m 



where the electric field E and magnetic field B are expressed in terms of the potentials in 
the usual form, E = — V0 — dA/dt, B = V x A. Dynamical equations ()2.7flj) can again be 
separated into the form 

d'^q e 



dt"^ m 
d'^r r, 



[E{t,r) + ux B{t,r)), (2.71) 



+ u'{r-q)=0. (2.72) 
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The center of mass q satisfies Newton's equations under the action of the total external Lorentz 
force, while point r still satisfies the isotropic harmonic motion of angular frequency u around 
point q. But the external force and the fields are defined at point r and not at point q. It is the 
velocity u of point r that appears in the magnetic term of the Lorentz force. Point r clearly 
represents the position of the charge. In fact, this minimal coupling we have considered is the 
coupling of the electromagnetic potentials with the particle current, that in the relativistic case 
can be written as j^A^, but the current is associated to the motion of a charge e at point r. 




Figure 2.1: Charge motion in the CM. frame. 

This charge has an oscillatory motion of very high frequency u that, in the case of the 
relativistic electron is = 2mc^ /h ~ 1.55 x lO^^s^^. The average position of the charge is 
the center of mass, but it is this internal orbital motion, usually known as the zitterbewegung, 
that gives rise to the spin structure for this model and also to the magnetic properties of the 
particle, as we shall see later. 

When analyzed in the center of mass frame (see Fig. 12. HI . q = 0, r = fc, the system reduces 
to a point charge whose motion is in general an ellipse, but if we choose C = D, and C ■ D = 0, 
it reduces to a circle of radius a = C = D, orthogonal to the spin. Then if the particle has 
charge e, it has a magnetic moment that according to the usual classical definition is: ^ 

^ = -y.x,dV = -fcx- = --S, (2.73) 

where j = e6^{r — k)dk/dt is the current associated to the motion of a charge e located at 
point k. The magnetic moment is orthogonal to the zitterbewegung plane and opposite to the 
spin if e > 0. It also has a non- vanishing oscillating electric dipole d = ek, orthogonal to n 
and therefore to S in the center of mass frame, such that its time average value vanishes for 
times larger than the natural period of this internal motion. Although this is a nonrelativistic 
example it is interesting to point out and compare with Dirac's relativistic analysis of the 

^ J.D. Jackson, Classical Electrodynamics, John Wiley & Sons, NY (1998), 3rd. ed. p. 186. 
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electron, ^ in which both momenta /x and d appear, giving rise to two possible interacting 
terms in Dirac's Hamiltonian. We shall come back to this analysis later when we study the 
elementary relativistic particles. 



2.2.4 Spinning Galilei particle with orientation 

Another simple example of spinning particles is the one in which the spin is related only 
to the angular variables that describe orientation. 

Let us assume now a dynamical system whose kinematical space is X = Q/M.^, where 
M.I = {M^, +} is the 3-parameter Abelian subgroup of pure Galilei transformations. Then, 
the kinematical variables are x = {t,r,p), which are interpreted as the time, position and 
orientation respectively. 

The Lagrangian for this model takes the general form 

L = Ti + R-r + W -u:. 

Because of the structure of the exponent ()2.167jl . the gauge function for this system can be 
taken the same as before. The general relationship ()2.46|) leads to x a; = 0, because the 
Lagrangian is independent of it, and therefore W and u) must be collinear. According to the 
transformation properties of the Lagrangian, the third term W ■ uj is Galilei invariant and 
since W and uj are collinear, we can take ~ and one possible Lagrangian that describes 
this model is of the form: 

^ mr I UJ , , 

The different Noether's constants are 



m 



fdrV I 



K = mr - Pt, J = r X P + W, 

where u = dr/dt is the velocity of point r, and Q, = uj/i is the time evolution angular 
velocity. Point r is moving at a constant speed and it also represents the position of the 
center of mass. The spin is just the observable S = W that satisfies the dynamical equation 
dS/dt = UJ X S = 0, and thus the frame linked to the body rotates with a constant angular 
velocity fl. 

The spin takes the constant value S = Jfi, whose absolute value is independent of the 
inertial observer and also the angular velocity Q, = uj/i is constant. The parameter / plays the 
role of a principal moment of inertia, suggesting a linear relationship between the spin and the 
angular velocity, which corresponds to a particle with spherical symmetry. The particle can 
also be considered as an extended object of gyration radius Rq, related to the other particle 
parameters by / = niRl. 

This system corresponds classically to a rigid body with spherical symmetry where the 
orientation variables p can describe for instance, the orientation of its principal axes of inertia 
in a suitable parametrization of the rotation group. This is a system of six degrees of freedom. 
Three represent the position of the center of charge r and the other three p, represent the 
orientation of a Cartesian frame linked to that point r. Since for this system there is no 

^ P. A.M. Dirac, The Principles of Quantum mechanics, Oxford Univ. Press, 4th ed. (1967). 
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dependence on the acceleration, the centers of mass and charge will be represented by the 
same point. 

In the center of mass frame there is no current associated to this particle and therefore it 
has neither magnetic nor electric dipole structure. As seen in previous examples, all magnetic 
properties seem therefore to be related to the zitterbewegung part of spin and are absent in 
this rigid body-like model. 

Relativistic particles 

2.3 Relativistic point particle 

See the Appendix about the Poincare group at the end of this chapter for the group notation 
used throughout this section. 

The kinematical space is the quotient structure X = V/C, where V is the Poincare group 
and the subgroup £ is the Lorentz group. Then every point x G X is characterized by the 
variables x = (t('r), t(t)), with domains t G M, r G as the corresponding group parameters, 
in such a way that under the action of a group element g = (6, a, v, /j,) of V they transform as: 

t'ir) = -rt{T) +^{v R{^i)r{T))/c^ + b, (2.75) 

r'(r) = i?(^)r(r)+7t;t(r)+ / (-t; ■ i?(/x)r(r))t; + a, (2.76) 

and are interpreted as the time and position of the system. If, as usual, we assume that the 
evolution parameter r is invariant under the group, taking the r-derivatives of ()2.75j) and ()2.76|1 

we get 

t'ir) = ^i{T)+^{vR{fi)r{T))/c'', (2.77) 
r'(r) = R{^)r{T)+^viiT) + —^—ivRifi)rir))v. (2.78) 

The homogeneity condition of the Lagrangian, in terms of the derivatives of the kinematical 
variables, reduces to three the number of degrees of freedom of the system. This leads to the 
general expression 

L = Ti + R-r, (2.79) 

where T = dL/dt and Ri = dL/dvi, will be functions of t and r and homogeneous functions 
of zero degree of i(r) and T'(r). Because the Lagrangian is invariant under V, the functions T 
and R transform under the group V in the form: 

r = -fT --f{v- R{fi)R), (2.80) 

R' = R{fi)R--fvT/c^ + -^{v ■ R{n)R)v/c^. (2.81) 

1 + 7 

We thus see that T and R are invariant under translations and therefore they must be functions 
independent of t and r. 

The conjugate momenta of the independent degrees of freedom = are pi = dL/dri, 
and consequently Noether's theorem (ll.25jl leads to the following constants of the motion, 
that are calculated similarly as in the Galilei case except for the invariance under pure Lorentz 
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transformations. We have now no gauge function and the variations are 6t = r ■ 6v/c^, 
Mi = Ti/c^ and 5r = t6v, Mij = t6ij and thus we get: 

Energy H = -T, (2.82) 

linear momentum P = R = p, (2.83) 

kinematical momentum K = Hr/c^ — Pt, (2.84) 

angular momentum J = r x P. (2.85) 

The energy and the linear momentum transform as: 

H\r) = jH{T)+j{vR{^j,)P{r)), (2.86) 

2 

P'(r) = R{t^)Pir) + ^HiT) + j^^-^iv.R{fi)Pir))v. (2.87) 

They transform like the contravariant components of a four- vector P^ = {H/c,P). The 
observables cK and J are the essential components of the antisymmetric tensor J^*^ = — J'^'" = 
x'^P" - x'^P^, cKi = and Jk = eniJ'^ /2. 

Taking the r derivative of the kinematical momentum, K = {], we get P = Hr/c^t = 
Hu/c^, where u = r/i is the velocity of the particle and the point r represents both the center 
of mass and center of charge position of the particle. 

The six conditions P = and K = 0, imply u = and r = 0, so that the system is at 
rest and placed at the origin of the reference frame, similarly as in the nonrelativistic case. We 
again call this class of observers the center of mass observer. 

From ()2.86p and ()2.87|1 we see that the magnitude {H/cY — P^ = rn^c^ is a Poincare 
invariant and a constant of the motion. This defines the mass of the particle which we take as 
a positive number m. By using the expression of P = Hu/c?, we get 



and the sign of H, which is another Poincare invariant property, can be either positive or 
negative. The velocity u < c, otherwise H will be imaginary. If u > c the invariant (H/c)^ — 
P^ < and it is not possible to define the rest mass of the system. By substitution of the 
found expressions for T and R in ()2.79|] . there are two possible Lagrangians for a point particle 
of mass m, characterized by the sign of H 

L = '^mc\/ cH"^ — r^. (2. 
Expansion of this Lagrangian to lowest order in u/c, in the case of positive H, we get 

2 i 



o- m r 

L = —mc t + 



where the first term —mc^i that can be withdrawn is just the equivalent to the Galilei in- 
ternal energy term —HQt of ()2.17|1 . The Lagrangian with H < has as nonrelativistic limit 
— {m/2)r'^/i which is not obtained in the Galilei case. 

The spin of this system, defined similarly as in the nonrelativistic case, 



c2 



S = J - qx P = J - —K X P = 0, (2.89) 



vanishes, so that the relativistic point particle is also a spinless system. 
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2.4 Relativistic spinning particles 

There are three maximal homogeneous spaces of V, all of them at first parameterized by 
the variables (t, r, u, p), where the velocity variable u can be either u < c, u = c or u > c. We 
shall call these kinds of particles by the following names: The first one, since the motion of the 
position of the charge r satisfies m < c, we call a Bradyon, from the Greek term I3pa5v<, = 
slow. Bradyons are thus particles for which point r never reaches the speed of light. The 
second class of particles {u = c) will be called Luxons because point r is always moving at 
the speed of light for every observer, and finally those of the third group, because u > c, are 
called Tachyons, from the Greek raxvt; = fast. 

For the second class we use the Latin denomination Luxons in spite of the Greek one of 
photons, because this class of particles will supply the description not only of classical photons 
but also a classical model of the electron. This class of models is very important and it has 
no nonrelativistic limit. Therefore the models this manifold produce have no nonrelativistic 
equivalent. 

The first class corresponds to a kinematical space that is the Poincare group itself and 
produces models equivalent to the ones analyzed in the non-relativistic case. To describe the 
classical electron and the photon we shall consider next the case of luxons. 



2.5 Luxons 

Let us consider those mechanical systems whose kinematical space is the manifold X gen- 
erated by the variables {t, r, u, p) with domains t G M, r G M^, p G as in the previous case, 
and n G but now with u = c. Since m = c we shall call this kind of particles Luxons. This 
manifold is in fact a homogeneous space of the Poincare group P, and therefore, according to 
our definition of elementary particle has to be considered as a possible candidate for describ- 
ing the kinematical space of an elementary system. In fact, if we consider the point in this 
manifold x = (0, 0, u, 0), the little group that leaves x invariant is the one-parameter subgroup 
Vu of pure Lorentz transformations in the direction of the vector u. Then X ~ V/Vu, is a 
nine-dimensional homogeneous space. 

For this kind of systems the variables t, r transform according to (|2.75jl and (|2.76jl . re- 
spectively and the derivatives as in ()2.77p and ()2.78|) . For the velocity u the transformation is 
obtained from ()2.178p and is 

2 

= ViTi W2^ • (2-90) 

7(1 + v ■ R[p) u{t)/c^) 

The general transformation of the orientation variables p are obtained from (j2.179jl but now 
the functions F and G, which involve some 7(m) factors, become infinite and in the limit m — > c 
they take the form 

N _ + P(^) + X P(^) + Fcjv, P; u{t), p{t)) 

^ l-M-p(r) + G,Kp;n(r),p(r)) ' ^ ^ 

where the functions Fc and Gc are given now by: 

F^{v,p;u,p) = /''^^x 2 [uxv + u{v- p) + v{u ■ p) 

(1 + 7(t;jjc^ 
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+ u X (v X n) + (u X p) X V + {u ■ p){v X fl) 

+ (ux p){v ■ p) + {ux p) X (v X p)], (2.92) 

Gc{v,^i]u,p) = — — ——[u-v + u-{vxp) + v-{uxp) 

(1 + 7(t;))c^ 

- {u- p){v- p) + {ux p)-{vx p)]. (2.93) 

Since u' = u = c, the absolute value of the velocity vector is conserved and it means that u' 
can be obtained from u by an orthogonal transformation, so that the transformation equations 
of the velocity under V can be expressed as: 

u' = R{(i))u, (2.94) 

where the kinematical rotation of parameter (f) is 

^ _ p + F^{v,P]u{t),Q) 

1 + G>,^;x.(r),0)- ^'-''^ 

In this case there also exist among the kinematical variables the constraints u = r/t. 
Equation ()2.9H1 also corresponds to 

R{p')=R{cf>)R{p), (2.96) 

with the same in both cases, as in ()2.95|l . 

Since the variable u{t) = c, during the whole evolution, we can distinguish two different 
kinds of systems, because, by taking the derivative with respect to r of this expression we get 
■u(r) ■ u{t) = 0, i.e., systems for which w = or massless systems as we shall see, and systems 
where -it 7^ but always orthogonal to u. These systems will correspond to massive particles 
whose charge internal motion occurs at the constant velocity c, although their center of mass 
moves with velocity below c. 



2.5.1 Massless particles. (The photon) 

If Tx = 0, n is constant and the system follows a straight trajectory with constant velocity, 
and therefore the kinematical variables reduce simply to (t, r, p) with domains and physical 
meaning as usual as, time, position and orientation, respectively. The derivatives t and r 
transform like ()2.77l) and ()2.78l) and instead of the variable p we shall consider the linear 
function a; defined in ()2.29ll that transforms under V: 

iv'ir) = Rict^Mr), (2.97) 

where, again, <p is given by ()2.95|) . 

In fact, from Ij2.96|l . since it = 0, taking the r-derivative. 

Rip') = R{4>)R{p), 

the antisymmetric matrix Vl = R{p)R^{p) has as essential components the angular velocity a;, 

(0 UJy \ 

LOz -w^ . (2.98) 
-LOy UJx / 
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It transforms as 

n' = R{p')R^{p') = R{cf))R(p)R^{p)R^{(j>) = R{(t})nR'^{(j>), 
and this matrix transformation leads for its essential components to II2.97|I . 

For this system there are no constraints among the kinematical variables, and, since it = 0, 
the general form of its Lagrangian is 

L = Ti + R-r + W-u;. (2.99) 

Funtions T = dL/di, = dL/dr\ Wi = dLjduf-, will depend on the variables {t,r,p) 
and are homogeneous functions of zero degree in terms of the derivatives of the kinematical 
variables {i, r, w). Since i they will be expressed in terms of u = r/i and fl = uj/i, which 
are the true velocity and angular velocity of the particle respectively. 

Invariance of the Lagrangian under V leads to the following transformation form of these 



functions under the group V: 

T' = -iT -^{v R{^i)R), (2.100) 

R! = Rin)R - -fvT/c^ + . , J v ■ R{n)R)v, (2.101) 

(1 + 7)c^ 

W = R{(f>)W. (2.102) 



They are translation invariant and therefore independent of t and r. They will be functions 
of only (p, n, ri), with the constraint u = c. Invariance under rotations forbids the explicit 
dependence on p, so that the dependence of these functions on p and p variables is only 
through the angular velocity uj. 

Noether's theorem gives rise, as before, to the following constants of the motion: 



Energy H = -T, (2.103) 

linear momentum P = R, (2.104) 

kinematical momentum K = Hr /c^ — P t — W x u/c^ , (2.105) 

angular momentum J = rxP + W. (2.106) 



In this case the system has no zitterbewegung because the Lagrangian does not depend on 
n which vanishes. The particle, located at point r, is moving in a straight trajectory at the 
speed of light and therefore it is not possible to find an inertial rest frame observer. Although 
we have no center of mass observer, we define the spin hjS = J — rxP = W. 

If we take in (12.1061) the r-derivative we get dS/dt = P x u. Since P and u are two 
non-vanishing constant vectors, then the spin has a constant time derivative. It represents a 
system with a continuously increasing angular momentum. This is not what we understand by 
an elementary particle except if this constant dS/dt = 0. Therefore for this system the spin is 
a constant of the motion and P and u are collinear vectors. 

Energy and linear momentum are in fact the components of a four- vector and with the spin 
they transform as 

H' = -fH + -f{v R{^J,)P), (2.107) 

P' = R{fi)P + ^vH/c' + 2\ ■ (2-108) 

(1 + 7jc^ 

S' = R{(f>)S. (2.109) 
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The relation between P and u can be obtained from ()2.in5p . taking the r-derivative and the 
condition that the spin W is constant, K = = —Hr/c^ + Pt, i.e., P = Hu/c^. If we take 
the scalar product of this expression with u we also get H = P ■ u. 

Then, from ()2.in7jl and (12.1081) . an invariant and constant of the motion, which vanishes, 
is {H/cY — P^. The mass of this system is zero. It turns out that for this particle both H and 
P are non- vanishing for every inertial observer. Otherwise, if one of them vanishes for a single 
observer they vanish for all of them. By ()2.in9|l . S"^ is another Poincare invariant property of 
the system that is also a constant of the motion. 

The first part of the Lagrangian Tt + R ■ r = —Hi + P ■ r, which can be written as 
— {H — P ■ u)i = 0, also vanishes. Then the Lagrangian is reduced to the third term S ■ uj. 

We see from ()2.94p and 1)2.109^ that the dimensionless magnitude e = S ■ u/ Sc is another 
invariant and constant of the motion, and we thus expect that the Lagrangian will be explicitly 
dependent on both constant parameters S and e. Taking into account the transformation 
properties under Vofu,uj and S, given in ()2.94L (j2.97jl and ()2.1()9jl respectively, it turns out 
that the spin must necessarily be a vector function of u and uj. 

If the spin is not transversal, as it happens for real photons, then S = eSu/c where e = ±1, 
and thus the Lagrangian finally becomes: 



c / t 



(2.110) 



From this Lagrangian we get that the energy is H = —dL/di = S ■ fi, where Q, = uj/i 
is the angular velocity of the particle. The linear momentum is P = dL/dr = eSft/c, and, 
since P and u are parallel vectors, and u must also be parallel, and if the energy is definite 
positive, then O = eflu/c. 

This means that the energy H = Sfl. For photons we know that S = h, and thus H = 
hVt = hu. In this way the frequency of a photon is the frequency of its rotational motion 
around the direction of its trajectory. We thus see that the spin and angular velocity for 
H > particles have the same direction, although they are not analytically related, because 
S is invariant under V while Q is not. 

We say that the Lagrangian (12.11011 represents a photon of spin S and polarization e. A 
set of photons of this kind, all with the same polarization, corresponds to circularly polarized 
light, as has been shown by direct measurement of the angular momentum carried by these 
photons. ^ Left and right polarized photons correspond to e = 1 and e = — 1, respectively. 
Energy is related to the angular frequency H = HQ, and linear momentum to the wave number 
P = hk, that therefore is related to the angular velocity vector by A; = efi/c. If it is possible 
to talk about the 'wave-length' of a single photon this will be the distance run by the particle 
during a complete turn. 

2.5.2 Massive particles. (The electron) 

If we consider now the other possibility, -it 7^ but orthogonal to u, then variables i and r 
transform as in the previous case (I2T7II and (l278ll . but for ii and UJ we have: 

ii' = R{(t))u + R{(t))u, (2.111) 
J = R{(t))uj + uj^, (2.112) 



R. A. Beth, Phys. Rev. 50, 115 (1936). 
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where the rotation of parameter (f) is again given by ()2.95|1 and vector uj^ is: 

7_R'U X D — (7 — X it) + 27^(1; ■ _R(n X + 7)c^ 

= r^— 5-^; ■ (2.113) 

7(c"' + V ■ Ru) 

Expression ()2.11ip is the r-derivative of p.94p and can also be written in the form: 

u' = ^(^)^ (2 114) 

Expression ()2.112jl comes from -R(p') = R{<p)R{p) and taking the r-derivative of this expression 
R{p') = R{<p)R{p) +R{<p)R{p), because parameter (p depends on r through the velocity u{t), 
and therefore 

n' = R{p')R^{p') = R{(j))nR^{cj)) + i?(0)i?^(0). 

R{(f))VtRF {4>) corresponds to R{<j))uj and the antisymmetric matrix = R{<p)R^{<j)) has as 
essential components the uj^f, vector, i.e., equation (12.11311 . 

The homogeneity condition of the Lagrangian leads to the general form 

L = Ti + R-r + U- u + W-LjJ, (2.115) 

where T = dL/di^ Ri = dL/dr\ Ui = dL/dii^ and Wi = dL/duj^, and Noether's theorem 
provides the following constants of the motion: 

Energy H = -T - {dU /dt) ■ u, (2.116) 

linear momentum P = R—{dU/dt), (2.117) 

kinematical momentum K = Hr/c^ — Pt — Z x u/c^, (2.118) 

angular momentum J = rxP + Z. (2.119) 

In this case the function Z is defined as in the Galilei case, by 

Z = uxU + W. (2.120) 

Expressions (I2.116L I2.117jl imply that H/c and P transform like the components of a 
four- vector, similarly as in (12. 86112. 87|l , thus defining the invariant and constant of the motion 
{H/cf — P^ = rri^c^, in terms of the positive parameter m which is interpreted as the mass of 
the particle. 

Observable Z transforms as: 

2 

Z'{r) = iR{p)Z{r) - ■ R{l^)Z{r))v + J(^ x R{p){Z{t) x n)), (2.121) 

an expression that corresponds to the transformation of an antisymmetric tensor Z^^'^ with 
strict components = (Z x n)Yc, and Z""^ = e^-^^Z^. 

By defining the relative position vector k = Z x u/H, the kinematical momentum ()2.118jl 
can be cast into the form 

K = Hq/c^ - Pt, 

where q = r — k, represents the position of the center of mass of the particle. 
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The spin is defined as usual 



c2 



S = J - qx P = J - —K X (2.122) 

H 

and is a constant of the motion. It takes the form 

S = Z + k X P = Z + ^{Z X u) X P. (2.123) 

The helicity S P = Z P = J P, is also a constant of the motion. We can construct the 
constant Pauli-Lubanski four-vector 

wf" = {P ■ S,HS/c), (2.124) 

with —w^w^ = Tn^c^S"^, in terms of the invariant properties m and 5" of the particle. 

If we take in (|2.118|1 the r-derivative and the scalar product with the velocity u we get the 
Poincare invariant relation: 

This will give rise to Dirac's Hamiltonian, H = cP ■ ol + jSmc^ when expressed in the 
quantum case, in terms of the a and (3 Dirac matrices. Since col is usually interpreted as the 
local velocity operator u of the electron, ^ we have H = P ■u + [3mc^ and this relation suggests 
the identification 

(3 = — -Z ■ { — xu 

mc^ \ at 

Here all magnitudes on the right-hand side are measured in the center of mass frame. We shall 
come back to this relation after quantization of this system. 

The center of mass observer is defined by the conditions P = K = For this observer 
Z = S IS constant, H = mc^ and thus from ()2.118p we get 

^ Sxu, (2.126) 



and the internal motion takes place in a plane orthogonal to the constant spin S. The scalar 
product with u leads to r-dr/rft = 0, and thus the zitterbewegung radius is a constant. Taking 
the time derivative of both sides of (12.1261) . we obtain mc^u = {S x du/dt), because the spin 
is constant in this frame, we get that u and S are orthogonal and therefore 

S = muxr. (2.127) 

Since S and u = c are constant, the motion is a circle of radius Rq = S/mc. For the electron 
we take S = h/2, and the radius is h/2meC = 1.93 x 10^^^ m., half the Compton wave length 
of the electron. The frequency of this motion in the CM. frame is z/ = 2me(? /h = 2.47 x 10^'' 
s~\ and u = 27rz/ = 1.55 x 10^^ rad s~^. The ratio of this radius to the so-called classical 
radius Rd = e'^ /SneomeC^ = 1.409 x 10~^^ m, is precisely Rd/Ro = e'^/2eohc = 1/136.97 = a, 
the fine structure constant. 



J.J. Sakurai, Advanced Quantum Mechanics, Addison- Wesley Reading, MA (1967). 
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Motions of this sort, in which the particle is moving at the speed of light, can be found in 
early literature, but the distinction between the motion of center of charge and center of mass 
is not sufficiently clarified. ^' ^ 

Nevertheless, in the model we are analyzing, the idea that the electron has a size of the 
order of the zitterbewegung radius is a plausible macroscopic vision but is not necessary to 
maintain any longer, because the only important point from the dynamical point of view is the 
center of charge position, whose motion completely determines the dynamics of the system. 
In this form, elementary particles, the kind of objects we are describing, look like extended 
objects. Nevertheless, although some kind of related length can be defined, they are dealt 
with as point particles with orientation because the physical attributes are all located at the 
single point r. The dynamics of equation ()2.126|1 can be represented in figur d2.2l where we 
have separated the two contributions to the total spin Z = S = Su + S^^, related respectively 
to the orbital and rotational motion. 




Figure 2.2: Motion of the charge of the electron in the center of mass frame. 

The transformation equation for the function Z, ()2.121jl can also be written as 

Z' = j{l + v ■ R{fj,)u/c^)R{cl))Z, (2.128) 

and therefore Z ii = Z' -u and Z' -u' = 7(1 + v ■ R{ii)u/c^)Z ■ u. Since it is orthogonal to u 
and -it, for the center of mass observer, it is also orthogonal to u and ii for any other inertial 
observer. 

An alternative method of verifying this is to take the time derivative in II2.118|1 and II2.119|I , 
and thus 

Hu - c^P xu - Z X — =0, 

dt dt 

M. Mathisson, Acta. Phys. Pol. 6, 163 (1937); 6, 218 (1937) 

M.H.L. Weyssenhof, Acta Phys. Pol. 9, 46 (1947). M.H.L. Weyssenhof and A. Raabe, Acta Phys. Pol. 
9, 7 (1947); 9, 19 (1947). 
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dZ „ 

i.e., 

Z X ^ = (H-u-P)u. 

dt ^ ' 

and a final scalar product with Z , leads to {H — u ■ P)u ■ Z = 0. The first factor does not vanish 
since the invariant /c^ — = m}(? is positive definite and if H = u - P, then {u ■ P)'^/c^ — 
with u < c is always negative, then Z -u — 0. If we take the time derivative of this last expression, 
with the condition that dZ /dt is orthogonal to m, we obtain Z ■ ii = Q. The observable Z has 
always the direction of the non-vanishing vector tt x u for positive energy particles and the 
opposite direction for particles of negative energy. 

Equation (j2.125jl can be recast into the form 

—ci - P r - \z ■ (iixu) =0, 

c 

where the first two terms give rise to the invariant term P^x^ = mc^icm, and the third to the 
invariant relation 

Z ■ (iixu) = mcHcm- (2.129) 

Here tcm is the time observable measured in the center of mass frame, and the right-hand side, 
which is positive definite for particles, implies that Z has precisely the direction of it x u. In 
the case of antiparticles it has the opposite direction. 

We see that the particle has mass and spin, and the center of charge moves in circles at 
the speed of light in a plane orthogonal to the spin, for the center of mass observer. All these 
features are independent of the particular Lagrangian of the type ()2.115|1 we can consider. All 
that remains is to describe the evolution of the orientation and therefore its angular velocity. 
The analysis developed until now is compatible with many different possibilities for the angular 
velocity. The behaviour of the angular velocity depends on the particular model we work with. 

To end this section and with the above model of the electron in mind, it is convenient to 
remember some of the features that Dirac ^ obtained for the motion of a free electron. Let 
point r be the position vector on which Dirac's spinor ip{t,r) is defined. When computing the 
velocity of point r, Dirac arrives at: 

a) The velocity u = i/h[H, r] = ccx, is expressed in terms of cx matrices and writes, . . a 
measurement of a component of the velocity of a free electron is certain to lead to the result 
±c'. 

b) The linear momentum does not have the direction of this velocity n, but must be related 
to some average value of it: ... 'the xi component of the velocity, cai, consists of two parts, 
a constant part c^piH~^, connected with the momentum by the classical relativistic formula, 
and an oscillatory part, whose frequency is at least 2mc^/h, . . . '. 

c) About the position r: 'The oscillatory part of Xi is small, . . . , which is of order of 
magnitude h/ mc, . . .'. 

And when analyzing, in his original 1928 paper, the interaction of the electron with an 
external electromagnetic field, after performing the square of Dirac's operator, he obtains two 
new interaction terms: 

^ „ ieh „ 

^.B + ^cE, 2.130 

2mc 2mc 



^ P. A.M. Dirac, The Principles of Quantum mechanics, Oxford Univ. Press, 4th ed. Oxford (1967). 
1° P.A.M. Dirac, Proc. Roy. Soc. Lon. A117, 610 (1928). 
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where the electron spin is written as S = and 

^-^0 a- 

in terms of a-Pauli matrices and E and B are the external electric and magnetic fields, 
respectively. He says, 'The electron will therefore behave as though it has a magnetic moment 
{eh/2mc) S and an electric moment {ieh/2mc) ex. The magnetic moment is just that assumed 
in the spinning electron model' (Pauli model). 'The electric moment, being a pure imaginary, 
we should not expect to appear in the model. ' 

However, if we look at our classical model, we see that for the center of mass observer, 
there is a non- vanishing electric and magnetic dipole moment 

e ^ e , dk e , , 

d = ek = — -Sxu, fi=-kx — = -—Y, 2.131 

mc^ 2 at 2m 

where S is the total spin and Y = —mk x dk/dt is the zitterbewegung part of spin. The time 
average value of d is zero, and the average value of n is the constant vector /x. 

This classical model gives rise to the same kinematical prediction as the nonrelativistic 
model described in Sec l2.2.3l If the charge of the particle is negative, the current of Fig l2.2l 
produces a magnetic moment that necessarily has the same direction as the spin. If the 
electron spin and magnetic moments are antiparallel, then we need another contribution to 
the total spin, different from the zitterbewegung. All real experiments to determine very 
accurately the gyromagnetic ratio are based on the determination of precession frequencies, but 
these precession frequencies are independent of the spin orientation. However, the difficulty 
to separate electrons in a Stern-Gerlach type experiment, suggests to perform polarization 
experiments in order to determine in a direct way whether spin and magnetic moment for 
elementary particles are either parallel or antiparallel. We have suggested a couple of plausible 
experiments to determine the relative orientation between the spin and magnetic moment of 
free electrons and also for electrons in the outer shell of atoms^^ 

Another consequence of the classical model is that it enhances the role of the so-called 
minimal coupling interaction jp,A^. The magnetic properties of the electron are produced by 
the current of its internal motion and not by some possible distribution of magnetic dipoles, 
so that the only possible interaction of a point charge at r with the external electromagnetic 
field is that of the current j^, associated to the motion of point r, with the external potentials. 



2.6 The dynamical equation of the spinning electron 

We have seen that for relativistic particles with u = c and u and ii orthogonal vectors, the 
position vector r moves in circles according to the dynamical equation (j2.126jl in the center 
of mass frame, as depicted in figure 12.21 But this solution is independent of the particular 
Lagrangian we choose as an invariant function of the kinematical variables and their derivatives, 
which accomplish with this orthogonality n ■ = 0, requirement. We are going to analyze this 
dynamical equation for any arbitrary inertial observer. 

Let us consider the trajectory r(t), t G [^1,^2] followed by a point of a system for an 
arbitrary inertial observer O. Any other inertial observer O' is related to the previous one by 

^^M. Rivas, Are the electron spin parallel or antiparallel vectors ?, ^ArXiv -.physics /Oil 205 7\ 
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a transformation of a kinematical group such that their relative space-time measurements of 
any space-time event are given by 

t' = T{t, r; 5(1, ... , gr), r' = R{t, r; gi, . . . , gr), 

where the functions T and R define the action of the kinematical group G, of parameters 
{gi, . . . ,gr), on space-time. Then the description of the trajectory of that point for observer 
O' is obtained from 

t'{t) = Tit, rit); gi,..., gr), r'{t) = R{t, r{t); g^, . . . , gr), Vt G [h, ta]. 

If we eliminate t as a function of t' from the first equation and substitute into the second we 
shall get 

r'{t')=r'{t';g^,...,gr). (2.132) 

Since observer O' is arbitrary, equation (|2.132p represents the complete set of trajectories of 
the point for all inertial observers. Elimination of the r group parameters among the function 
r'it') and their time derivatives will give us the differential equation satisfied by the trajectory 
of the point. This differential equation is invariant by construction because it is independent 
of the group parameters and therefore independent of the inertial observer. If G is either the 
Galilei or Poincare group it is a ten-parameter group so that we have to work out in general up 
to the fourth derivative to obtain sufficient equations to eliminate the ten group parameters. 
Therefore the order of the differential equation is dictated by the number of parameters and 
the structure of the kinematical group. 



2.6.1 The relativistic spinning electron 

Let us assume the above electron model. Since the charge is moving at the speed of light 
for the center of mass observer O* it is moving at this speed for every other inertial observer 
O. Now, the relationship of space-time measurements between the center of mass observer and 
any arbitrary inertial observer is given by: 

t{t*;g) = -fit* + v- R{a)r*{t*)) + b, 

rit*;g) = Ria)r* it*) + -fvt* + — — iv ■ Ria)r*it*)) v + a. 

1 + 7 

With the shorthand notation for the following expressions: 

dr*it*) dK dV 
Kit*) = Ria)r*it*), Vit*) = i?(a)^^V^ = 3^ = 'K, 

^ ^ V y V \ J ^ ' dt* dt* dt* 

Bit*)=v-K, Ait*) = v-V = ^, ^ = -B 

^ ^ ' V y dt*' dt* 



we obtain 
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r(3) = ^3^^^^^, (-35(1 + A)K-{1 + A- 3B^)V+ 



1 + 7 



{A{1 + A) + 3B^)v] (2.135) 



r(^) = ^,^^l^y ((1 + A)il-2A- 3A' - 15B')K- 



B{7 + AA- SA^ - 15B'^)V- 

7 



1 + 7 



1-8A-9A^ -15B^)Bv] . (2.136) 



From this we get 



(r«-r«)' = 1, (r«-r(2))=0, (2.137) 

p)..(^)) = -{r^''-r^'') = ::^^^,, (2.138) 

(r(^)..(3)) = -lir^''-r^'')=::^^^,, (2-139) 

(r(3).r(3)) = (1-^^ + 35^), (2.140) 

(r(^).r(^)) = ^^,^^(-1 + 2^ + 3^^ + 95^), (2.141) 

(r(3)..W) = -,^^1^(1 + ^ + 35^) 45. (2.142) 

From equations (12.1 3811 - fl2.1 4011 we can express the magnitudes A, B and 7 in terms of these 
scalar products between the different time derivatives (r*^*) ■ r^^'^). The constraint that the 
velocity is 1 implies that all these and further scalar products for higher derivatives can be 
expressed in terms of only three of them. If the three equations ()2.133|l - p.l35p are solved in 



terms of the unknowns v, V and K and substituded into (I3.24jl . we obtain the differential 
equation satisfied by the charge position, for any arbitrary inertial observer 

(^^(2) . ^(2)) ^ 



^ 2(r(3) ■ r(3)) _ 3(r(^) ■ r^^))^ _ (^^^ .A (2) _ n 

(r(2) . r(2)) 4(r(2) .r(2))2 ' ^ J ~ °- ^^"^^^^ 

It is a fourth order ordinary differential equation which contains as solutions motions at the 
speed of light. In fact, if (r^^^ ■ r^^^) = 1, then by derivation we have (r^^^ ■ r^"^^) = and the 
next derivative leads to (r^"^^ ■ r^"^^) + (r^^^ ■ r'^^^) = 0. If we take this into account and make 
the scalar product of (j2.143|l with r^^^, we get (r^^) ■ r^^'') + 3(r(2) ■ r^^^) = 0, which is another 
relationship between the derivatives as a consequence of \r^^^\ = 1. It corresponds to a helical 
motion since the term in the first derivative r^^^ is lacking. 
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2.6.2 The center of mass 

The center of mass position is defined by 



(^(2) . ^(2))3/2 + (^(3) . ^(3)) _ ^^^^^ "-^^J^ 



(2.144) 



We can check that both q and q^^^ vanish for the center of mass observer. Then, the fourth 
order dynamical equation for the position of the charge can also be rewritten here as a system 
of two second order differential equations for the positions q and r 

^ =0' ^ =— (^r7j^(^-^)' (2.145) 

a free motion for the center of mass and a kind of central motion for the charge around the 
center of mass. 

For the non-relativistic electron we get in the low velocity case g*^^^ and |g — t| = 1, 
the equations of the Galilei case 

q(2)^0, r(^)=g-r. (2.146) 

a free motion for the center of mass and a harmonic motion around q for the position of the 
charge. 

2.6.3 Interaction with some external field 

The free equation for the center of mass motion g*^^^ = 0, represents the conservation of the 
linear momentum dP/dt = 0. But the linear momentum is written in terms of the center of 
mass velocity as P = m'-f{q^^^)q^^\ so that the free dynamical equations ()2.145|1 in the presence 
of an external field should be replaced by 

P^'^ = F, r(^) = (g - r) , (2.147) 

where F is the external force and the second equation is left unchanged because we consider, 
even with interaction, the same definition of the center of mass position. 

— = m7(g«)q(2) +^^(g(i))3(q(i) . q{2))g{i) 
dt 

we get 

,n7(g«)3(q(i).q(2)) = F-g« 

and by leaving the highest derivative q^'^^ on the left hand side we finally get the differential 
equations which describe the evolution of a relativistic spinning electron in the presence of an 
external electromagnetic field: 

^q(2) = _^[^; + ^(i)xS-q«([^; + r«xB]-g«)], (2.148) 

1 _ . ^(1) , , , , 
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2.7 Appendix: Galilei group 

The Galilei group is a group of space-time transformations characterized by ten parameters 
g = {b, a.v, a.). The action of on a space-time point x = {t,r) is given by x' = gx, and is 
considered in the form 

x' — exp{bH) exp(a • P) exp(v • K) exp(a • J) x 

as the action of a rotation followed by a pure Galilei transformation and finally a space and 
time translation. In this way all parameters that define each one-parameter subgroup are 
normal, because the exponential mapping works. Explicitly 

t' = t + b, (2.150) 
r' = R{a)r + vt + a, (2.151) 

and the composition law of the group g" — g'g is: 

b" = b' + b, (2.152) 

a" = R{a')a + v'b + a', (2.153) 

v" = R{a')v + v', (2.154) 

R{a") = R{cx')R{cx). (2.155) 

For rotations we shall alternatively use two different parametrizations. One is the normal 
or canonical parametrization in terms of a three vector a = an, where n is a unit vector along 
the rotation axis, and a e [0, tt] is the clockwise rotation angle in radians, when looking along 
n. Another, in terms of a three vector /* = ntan(Q;/2), which is more suitable to represent 
algebraically the composition of rotations. 

The rotation matrix i?(a) = exp(a • J) is expressed in terms of the normal parameters ctj 
and in terms of the antisymmetric matrix generators Jj which have the usual matrix represen- 
tation 

/O \ / 1\ /O -1 0' 

Ji = -1 , J2 = , J3 = 1 
\0 1 / \-l 0/ \0 

and satisfy the commutation relations [Jj, J^] = ^ikiJh such that if we write the normal param- 
eters a. = an in terms of the rotation angle a and the unit vector n along the rotation axis, 
it is written as 

R{ot)ij — Sijcosa + ninj{l — cos a) — eijkriksma, k — 1,2, 3. (2.156) 
In the parametrization /x = ntan(Q;/2), the rotation matrix is 

= 7-— T ((1 - + - 2eijfc/Xfe) , i,j,k= 1,2,3. (2.157) 

In terms of these variables, R{tJi") = is equivalent to 
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This can be seen in a simple manner by using the homomorphism between the rotation group 
and the group SU{2). The matrix generators of SU{2) are J = — icr/2 in terms of Pauli 
matrices cr. In the normal parametrization the rotation matrix exp(Q; ■ J) = exp{—ict ■ cr/2) 
is written in the form 

R{cx) = cos(a/2)I - t{n ■ cr) sin(a/2). 
By defining /x = ntan(a/2), this rotation matrix is expressed as 

Ri^j,) = -7=L= (I-zfi-a), (2.159) 

where I is the 2x2 unit matrix and in this form we can get the composition law (j2.158jl . 

If the rotation is of value tt, then eqs. 112.15611 or H2.157|l lead to 

R(n,7T)ij = —Sij + 2ninj. 

Even if the two rotations and R{lJ,') involved in (|2. 15811 are of value tt, although tan(7r/2) = 

oo, this expression is defined and gives: 

// n X n' 

n"tan{a"/2) = 



n ■ n' 



The absolute value of this relation leads to tan(a"/2) — tan 6*, i.e., a" = 29, where 6 is the angle 
between the two unit vectors n and n'. We obtain the known result that every rotation of value 
a around an axis n can be obtained as the composition of two rotations of value tt around two 
axes orthogonal to n and separated by an angle a/2. 

For the orientation variables we shall use throughout the book the early Greek variables 
a, /3, . . . whenever we consider the normal parametrization, while for the tan(a/2) parametriza- 
tion we will express rotations in terms of the intermediate Greek variables p, . . . . In this 
last notation, transformation equations (|2. 152112. 155|1 should be replaced by 

b" = b' + b, (2.160) 
a" = R{fi')a + v'b + a', (2.161) 
v" = R{fi')v + v', (2.162) 

„ p,' + fl + ^i' X ^ . . 

M = • 2.163 

The neutral element of the Galilei group is (0, 0, 0, 0) and the inverse of every element is 
(6, a, V, cx)^^ = {—b, —R{—a.){a — bv), —R{—a)v, —a). 
The generators of the group in the realization ()2.15n[ I2.151j) are the differential operators 

H = d/dt, P, = d/dn, Ki = td/dri, Jk = ekHrid/dn (2.164) 
and the commutation rules of the Galilei Lie algebra are 

[J,J] = -J, [J,P] = -P, [J,K] = -K, [J,i/] = 0, (2.165) 
[H, P] = 0, [H, K] = P, [P, P] = 0, [K, K] = 0, [K, P] = 0. (2.166) 



12 



D. Hestenes, Space-time algebra, Gordon and Breach, NY (1966). 
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All throughout this book, except when explicitly stated, we shall use the following shorthand 
notation for commutators of scalar and 3-vector operators, that as usual, are represented by 
bold face characters: 



[A,B] 


= C, = 






[A,B] 


= c, = 




= 6ijC, 


[A,B] 


= c, = 




= Ci, 


[B,A] 


= c, = 


^ [B,Ai] 





where 6ij = 6ji is Kronecker's delta and eijk is the completely antisymmetric symbol, so that 
Latin indexes match on both sides of commutators. 

The group action (j2.15()jl - (l2. 15111 represents the relationship between the coordinates (t, t) 
of a space-time event as measured by the inertial observer O and the corresponding coordinates 
{f, r') of the same space-time event as measured by another inertial observer O' . The ten group 
parameters have the following meaning. If we consider the event (0, 0) measured by O, for 
instance the flashing of a light beam from its origin at time t = 0, it takes the values (6, a) 
in O', where h is the time parameter that represents the time translation and a is the space 
translation. The parameter v of dimensions of velocity represents the velocity of the origin of 
the Cartesian frame of O as measured by O', and finally the parameters a, or R{ol), represent 
the orientation of the Cartesian frame of O as measured by O' . In a certain sense the ten 
parameters (6, a, d, cm.) with dimensions respectively of time, position, velocity and orientation 
describe the relative motion of the Cartesian frame of O by O' . 

The Galilei group has non-trivial exponents given by 

^{g, g') = m (^v%' + v ■ R{a)a'^ . (2.167) 

They are characterized by the non- vanishing parameter m. 

The central extension of the Galilei group is an 11-parameter group with an additional 
generator / which commutes with the other ten, 

[/, H] = [I, P] = [I, K] = [I, J] = 0, (2.168) 

and the remaining commutation relations are the same as above ()2.165| I2.166|l , except the last 
one which appears as 

[Ki,Pj] = -m5ijl, or [K,P] = -mI, (2.169) 

using our shorthand notation, in terms of a non-vanishing parameter m. If we define the 
following polynomial operators on the group algebra 

W = IJ-—KxP, U = IH-—P^, (2.170) 
m 2m 

U commutes with all generators of the extended Galilei group and W satisfies the commutation 
relations: 

[W,W] = -IW, [J,W] = -W, [W,P] = [W,K] = [W,H]=0, 

13 V. Bargmann, Ann. Math. 5, 1 (1954). 

1^ J.M. Levy-Leblond, Galilei Group and Galilean Invariance, in E.M. Loebl, Group Theory and its applica- 
tions, Acad. Press, NY (1971), vol. 2, p. 221. 
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so that W'^ also commutes with all generators. It turns out that the extended Galilei group 
has three functionally independent Casimir operators which, in those representations in which 
the operator / becomes the unit operator, for instance in irreducible representations, are in- 
terpreted as the mass, M = ml, the internal energy Hq = H — P^/2m, and the absolute value 
of spin 

= (j-—K X p] . (2.171) 

\ m J 

The spin operator S in those representations in which / = I, satisfy the commutation relations: 

[S, S] = -S, [J, S] = -5, [S, P] = [5, H] = [S, K] = 0, 

i.e., it is an angular momentum operator, transforms like a vector under rotations and is 
invariant under space and time translations and under Galilei boosts, respectively. It reduces 
to the total angular momentum operator J in those frames in which P = K = 0. 



2.8 Appendix: Poincare group 

The Poincare group is the group of transformations of Minkowski's space-time that leave 
invariant the separation between any two close space-time events ds"^ = rj^ydx'^dx'' . We shall 
consider the contravariant components = (ct, r), and x' = gx is expressed as x"^ = K'^y x" + 
a^, in terms of a constant matrix A and a constant translation four- vector a'^ = {cb,a). We 
take for the covariant components of Minkowski's metric tensor 77^,^ = diag(l, — 1, — 1, — 1). 
Then dx'^ = A^^dx^ and ds"^ = rj^j^ydx'^dx'" = rj^pdx'^dxP implies for the matrix A 

?7^,A^A% = ri^p. (2.172) 

Relations ()2.172p represent ten conditions among the 16 components of the matrix A, so 
that each matrix depends on six essential parameters, which can be chosen in many ways. 
Throughout this book we shall take three of them as the components of the relative velocity 
V between inertial observers and the remaining three as the orientation o: of their Cartesian 
frames, expressed in a suitable parametrization of the rotation group. 

Therefore, every element of the Poincare group V will be represented, as in the previous 
case of the Galilei group, by the ten parameters g = {b, a, v, a.) and the group action on a 
space-time point x = (t, r) will be interpreted in the same way, i.e., x' = gx: 

x' = exp{bH) exp(a ■ P) exp(/3 ■ K) exp{cx ■ J)x, (2.173) 

as the action of a rotation followed by a boost or pure Lorentz transformation and finally a 
space and time translation. It is explicitly given on the space-time variables by 

t' = -ft + j{v-R{n)r)/c^ + b, (2.174) 
r' = R{n)r + -fvt + -f^{v ■ R{fi)r)v/{l + -f)c^ + a. (2.175) 

Parameter (3 in ()2.173|) is the normal parameter for the pure Lorentz transformations, that 
in terms of the relative velocity among observers v is expressed as /3//9tanh/5 = v/c as we 
shall see below. The dimensions and domains of the parameters b, a and /x are the same 
as those of the Galilei group, and the parameter G M^, with the upper bound v < c, has 
also dimensions of velocity. The physical meaning of these ten parameters, that relate any two 
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inertial observers, is the same as in the Galilei case. The parameter v is the velocity of observer 
O, as measured by O', and -R(Ai) represents the orientation of Cartesian frame O relative to 
O', once O' is boosted with velocity v. The factor 7(t>) = (1 — f^/c^)~^/^. 

The composition law of the group is obtained from x" = A'x' + a' = A' (Ax + a) + a' that by 
identification with x" = A"x + a" reduces to A" = A'A and a" = A' a + a', i.e., the composition 
law of the Lorentz transformations, that we will find in the next Section [2. 8. H and a Poincare 
transformation (A', a') of the four-vector a^. In this parametrization g" = g'g, is: 

b" = ih + i{v' ■R{ti')a)/c^ + h', (2.176) 

'2 

a = R{fj,')a + ^'v'b+—^-—^{v'-R{fi')a)v' + a', (2.177) 

(1 + 

R{fi')v + yv' + ,,7' ■ R{li')v)v' 

-i'{l + v' ■R{ti')v/c^) ' ^ ' 

" ^ n' + ti + fx' X n + F{v', fi', V, ^) 

where F{v', /x', v, fi) and G{v', fj,', v, /x) are the real analytic functions: 

F{v',fi',v,fi) = [vxv' + v{v' ■ ti') + v'{v- fi) 

[1 + 7)(1 + 7')c"^ 

+ V X {v X /x') + {v X ^) X v' + {v ■ ^i){v' X fl') 

+ {v X n){v' ■ fx') + {v X fi) X {v' X fi')] , (2.180) 

Giv',^J,',v,^^) = ' [vv' + v-iv'xn') + v'-ivx^) 

(1 + 7)(1 + 7')c"= 

- {v- n){v' ■ n') + {v X n) ■ (v' X n')]. (2.181) 

The unit element of the group is (0, 0, 0, 0) and the inverse of any arbitrary element 
(6, a, V, fi) is 

(-76 + -fv ■ a/c^, -R{-n){a - -fvb + — — —{v ■ a)v), -R{-n)v, -/x). 

(1 + 7)c"= 

The group generators in the realization (j2.174[ I^.175j) . and in terms of the normal param- 
eters (6, a, f3, a), are 

H = d/dt, P, = d/du, Ki = ctd/dn + {n/c)d/dt, Jk = Eki'rid/dn. 

Thus, K and J are dimensionless and the commutation relations become 

[J, J] = -J, [J, P] = -P, [J, K] = -K, [J, H] = 0, [H, P] = 0, (2.182) 

[H, K] = cP, [P, P] = 0, [K, K] = J, [K, P] = -H/c. (2.183) 

If, as usual, we call = ct, Pq = H/c, and Ki = J^i = —Jio and Jk = ^^kirJir, a;^ = rj^^x^ , 
= 0, 1, 2, 3 and dy = d/dx" , then. 



1^ M.Rivas, M.Valle and J.M.Aguirregabiria, Eur. J. Phys. 6, 128 (1986). 
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In covariant notation the commutation relations appear: 

[P,iPu] = 0, 
\Jfj.yi P<j\ ^^.aPu ~l~ ^uaPpi 

\Jhui J P(t\ ^HpJva ^uaJpp ~l~ ^upJpa ~l~ ^fiaJup- 

The Poincare group has two functionally independent Casimir invariants. One is interpreted 
as the squared mass of the system, 

P^^P^ = {H/cf - p2 = m^c^ (2.184) 

and the other is the square of the Pauli-Lubanski four- vector W^. The Pauli-Lubanski four- 
vector is defined as 

W^" = le""''^ PyJ^x = {P J, HJ/c -KxP) = {P-S, HS/c), (2.185) 

which is by construction orthogonal to P^, i.e., W^P^ = 0. 

It is related to the spin of the system S, defined through the relation 

HS/c = HJ/c- K X P, (2.186) 

so that its time component W'^ = P ■ S = P ■ J is the helicity of the particle, and the spatial 
part is the vector (12.18611 . 

The other Casimir operator is thus 

W^'W^ = (P • J)2 - {HJ/c -KxPf = -rn^c^S'^, (2.187) 

where it depends on S"^, the absolute value squared of the spin. We see in the relativistic case 
that the two parameters m and S characterize the two Casimir invariants and therefore they 
are the intrinsic properties of the elementary particle the formalism provides. In the quantum 
case, since the representation must be irreducible S"^ = s{s + for any s = 0, 1/2, 1, . . ., 
depending on the value of the quantized spin of the particle, but in the classical case S'^ can 
take any continuous value. 

These operators satisfy the commutation relations: 

[W", W] = e^'^'^PW^Pp, (2.188) 

where we take e°^^^ = +1, and 

[P^ W] = 0, [M,,, W^] = -v^^W, + 7],^W,. (2.189) 

The Poincare group has no non-trivial exponents, so that gauge functions when restricted to 
homogeneous spaces of V vanish. 

2.8.1 Lorentz group 

The Lorentz group £ is the subgroup of transformations of the form (0, 0, v, /j,), and every 
Lorentz transformation A{v,iJ,) will be interpreted as A{v,fj,) = L{v)R{fj,), as mentioned 
before where L{v) is a boost or pure Lorentz transformation and Ri^fi) a spatial rotation. 
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Expressions (12.1 781 l2T79|l come from A{v",n") = A{v' , n')A{v, fx). Expression (12.17811 is the 
relativistic composition of velocities since 

L{v")R{fj,") = L{v')R{fi')L{v)R{fj,) 

= L{v')R{fi')L{v)R{-fi')R{ti')R{^i), 

but the conjugate of the boost R{^')L{y)R{—^') = L{R{^')v) is another boost and thus 

L{v")R{fi") = L{v')L{R{^i')v)R{tJ,')R{fi). 

The product L{y')L{R{^')v) = L{v")R{w) where v" is the relativistic composition of the 
velocities v' and R{fi')v, and R{w) is the Thomas- Wigner rotation associated to the boosts 
L{v') and L{R{fi')v). 

Therefore, expression (12 .1781] is equivalent to 

L{v") = L{v')L{R{^i')v)R{-w), (2.190) 

and (12.1791) is 

Rifj,") = R{w)R{n')R{fi) = R{ct))R{^). (2.191) 
The Thomas- Wigner rotation matrix R{w) is: 



+ — iv,R^,Vk - v^R,,Vk) + JY^yJ^YT^'''^''''' 

and the factor 

i' = ii ( 1 + — 

Matrix R{w) is written in terms of the vector parameter to, which is a function of v /x' and 
V, given by 

_ F{v\Q^R{^.')v^Q) 
l + G{v',0,R{^x')v,Oy 
and the parameter 0, such that R{4>) = R{w)R{fx') is 

If any one of the two velocities v or v' vanishes, R{w)ij = 6ij. 

The composition law is obtained by the homomorphism between the Lorentz group £ and 
the group SL{2, C) of 2 x 2 complex matrices of determinant +1. The Lie algebra of this group 
has as generators J = —ia/2 and K = cr/2, where cxj are Pauli spin matrices. A rotation of 
angle a around a rotation axis given by the unit vector n is given by the 2x2 unitary matrix 
exp(Q; ■ J), 

R{a) = cos(a/2)ao - in ■ crsin(a/2). (2.194) 
In terms of the vector /Lt = tan(a/2)n, 

R{ti) = —L={ao - ■ a), (2.195) 
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where ctq is the 2x2 unit matrix. A pure Lorentz transformation of normal parameters /3j is 
represented by the hermitian matrix exp(/3 ■ K). This matrix is: 



L(/3) = cosh(/3/2)ao + ^-^ sinh(/3/2). 

P 



(2.196) 



In terms of the relative velocity parameters, taking into account the functions cosh/3 = 7(f), 
sinh/3 = 7f/c and the trigonometric relations cosh(/?/2) = a/ (cosh (3 + l)/2 and tanh(/3/2) = 
sinh/5/(l + cosh/3), the matrix can be written as 



Liv] 



1 + 7 



ao + . 2.197 

1 + 7 c y 

Then, every element of SL{2, C) is parametrized by the six real numbers {v, /x), and inter- 
preted as 

A{v,n) = L{v)R{n). (2.198) 

We thus see that every 2x2 matrix A G SL{2, C) can be written in terms of a complex 
four- vector a'* and the four Pauli matrices cr^. As A = a'^cr^, and detA = 1 leads to a^a^ = 1 
or (a°)^ — = 1. The general form of (j2.198jl is 



cr ■ V 



A{v,^i) 



1 + 7 
2(l + /i2 



(To 1 - i 



1+7 



+ cr 



u + rt X ^ 

1 + 7 



(2.199) 



7(f )v/c. 



here the dimensionless vector u 

Conversely, since Tr (cr^a^) = 25^y, we obtain = (l/2)Tr (Aa^j). If we express ()2.199p in 
the form A{v, /x) = a^a^ we can determine fi and v, and thus n, from the components of the 
complex four-vector a'^ as: 

Im (a) 



Re (a°) 



u = 2 [Re (a°)Re (a) + Im (a°)Im (a) + Re (a) x Im (a 



(2.200) 
(2.201) 



where Re (a'^) and Im (a^) are the real and imaginary parts of the corresponding components 
of the four- vector a^. When Re (a°) =0 expression (|2.2()()jl is defined and represents a rotation 
of value 71 along the axis in the direction of vector Im (a). 

If we represent every Lorentz transformation in terms of a rotation and a boost, i.e., in the 
reverse order, A{v,iJ,) = R{fi)L{v), then the general expression of A is the same as ()2.199p 
with a change of sign in the cross product term u x /j,. Therefore, the decomposition is also 
unique, the rotation R{fi) is the same as before but the Lorentz boost is given in terms of the 
variables a'^ by 

n = 2 [Re (a°)Re (a) + Im (a°)Im (a) + Im (a) x Re (a)] . 

Note the difference in the third term which is reversed when compared with ()2.20H] . 

In the four-dimensional representation of the Lorentz group on Minkowski space-time, a 
boost is expressed as L{l3) = exp(/3 ■ K) in terms of the dimensionless normal parameters j3i 
and the 4x4 boost generators Ki given by 
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If we call B = (3 ■ K = PiKi, we have 
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with = + /3| + Pi and 5^ = and so on for the remaining powers of i?, so that the 
final expression for L{f3) = exp(/3 ■ K) is 
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where S = sinh/3 and C = cosh/5. What is the physical interpretation of Let us assume 
that observers O and O' relate their space-time measurements x and x' by x'^ = L{f3Y^x^ . 
Observer O sends at time t and at a later time t + dt two light signals from a source placed at 
the origin of its Cartesian frame. These two signals when measured by O' take place at points 
r' and r' + dr' and at instants t' and t' + dt\ respectively. Then they are related by 

cdt' = L^Qcdt, dx'^ = Uocdt 



because dx^ = 0. The quotient dx'^dt' is just the velocity of the light source v\ i.e., of the 
origin of the O frame as measured by observer O', and then this velocity = cUq/L^q = 
c{Pi/P)S/C, such that the relation between the normal parameters and the relative velocity 
between observers is 

V (3 



P 



tanh P 



and therefore tanh/3 = v/c. Function cosh/3 = 7(f) = (1 — v"^ /c^)~^/'^ and when the trans- 
formation is expressed in terms of the relative velocity it takes the form of the symmetric 
matrix: 
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(2.202) 



= L{—v). The orthogonal 4x4 rotation matrix takes the 
Ri^^)=ii (2.203) 



where R^n) is the 3x3 orthogonal matrix (j2.156|l . When a Lorentz transformation is expressed 
in the form K{v, /x) = L{v)R{^), then by construction the first column of A(v, ^t) is just the 
first column of (I2.2()2|] where the velocity parameters v are defined. Therefore, given the 
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general Lorentz transformation A{v,fj,), from its first column we determine the parameters v 
and thus the complete L{v) can be worked out. The rotation involved can be easily calculated 
as L{—v)A{v,^) = Ri^fi). If expressed in the reverse order A{v,^) = R{n)L{v), then it is 
the first row of A that coincides with the first row of ()2.2()2jl . It turns out that, given any 
general Lorentz transformation A{v, then A{v, /j,) = L{v)R{^) = R{^)L{v') with the same 
rotation in both sides as derived in (|2.2nn|) and L{v') = R{—^)L{v)R{^) = L{R{—fj,)v), i.e, 
the velocity v' = R{—fj,)v. In any case, the decomposition of a general Lorentz transformation 
as a product of a rotation and a boost is a unique one, in terms of the same rotation -R(^t) and 
a boost to be determined, depending on the order in which we take these two operations. 

Matrix A can be considered as a tetrad (i.e., a set of four orthonormal four-vectors, one 
time-like and the other three space-like) attached by observer O' to the origin of observer O. 
In fact, if the matrix is considered in the form A{v, fj,) = L{v)R{n), then the first column of 
A is the four- velocity of the origin of the O Cartesian frame and the other three columns are 
just the three unit vectors of the O reference frame, rotated with rotation -R(/x) and afterwards 
boosted with L{v). 



Chapter 3 

Quantization of the models 



Quantization of generalized Lagrangian systems will suggest that wave functions for ele- 
mentary particles must be squared integrable functions defined on the kinematical space. 

We shall use Feynman's quantization method to show the structure of the wave function 
and the way it transforms under the kinematical or symmetry group of the theory. Once the 
Hilbert space structure of the state space is determined, this leads to a specific representation 
of the generators of the group as self-adjoint operators and the remaining analysis is done 
within the usual quantum mechanical context, i.e., by choosing the complete commuting set of 
operators to properly determine a set of orthogonal basis vectors of the Hilbert space. Special 
emphasis is devoted to the analysis of the different angular momentum operators the formalism 
supplies. They have a similar structure to the classical ones, and this will help us to properly 
obtain the identification of the spin observable. 

The structure of the spin operator depends on the kind of translation invariant kinematical 
variables we use to describe the particle, and the way these variables transform under the 
rotation group. Since in the Galilei and Poincare case, as we have seen previously, these 
variables are the velocity u and orientation ot and they transform in the same way under 
rotations in both approaches, then the structure of the spin operator is exactly the same in 
both relativistic and nonrelativistic formalisms. 

As we have seen in the classical description the position of the charge of the particle 
and its center of mass are different points, and spin is related to the rotation and internal 
motion (zitterbewegung) of the charge around the center of mass of the particle. The magnetic 
properties of the particle are connected only with the motion of the charge and therefore to 
the zitterbewegung part of spin. It is this double spin structure that gives rise to the concept 
of gyromagnetic ratio when expressing the magnetic moment in terms of the total spin. If the 
Lagrangian shows no dependence on the acceleration, the spin is only of rotational nature, 
and the position and center of mass position define the same point. Spin 1/2 particles arise 
if the corresponding classical model rotates but no half integer spins are obtained for systems 
with spin of orbital nature related only to the zitterbewegung. On the manifold spanned by 
non-compact variables u no half-integer spins can be found, because the spin operator has the 
form of an orbital angular momentum and eigenvectors are but spherical harmonics. 

Dirac's equation will be obtained when quantizing the classical relativistic spinning particles 
whose center of charge is circling around its center of mass at the speed c. In that case, the 
internal orientation of the electron completely characterizes its Dirac algebra. 
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3.1 Feynman's quantization of Lagrangian systems 

Let us consider a generalized Lagrangian system as described in previous chapters and 
whose evolution is considered on the kinematical space between points xi and X2- 

For quantizing these generalized Lagrangian systems we shall follow Feynman's path inte- 
gral method. ^ The Uncertainty Principle is introduced in Feynman's approach by the condition 
that if no measurement is performed to determine the trajectory followed by the system from 
Xi to X2, then all paths x{t) are allowed with the same probability. Therefore a probability 
definition P[x(r)], must be given for every path. 

But instead of defining the probability associated to each possible path P[x(r)], this is 
calculated in terms of a probability amplitude, 0[a;(r)] for that path such that P[a;(r)] = 
|(/)[a;(r)]p, where < P < 1. But in general does not need to be a positive real number; in 
fact it is a complex number. Thus, to every possible trajectory followed by the system, x(r) 
in X space, Feynman associates a complex number 0[x(r)] called the probability amplitude of 
this alternative, given by 

0[x(r)] =iVexp|^^ L(x(r),i;(r))dr| = iVexp|^v4[,](xi,X2)| , (3.1) 

where is a path independent normalization factor, and where the phase of this complex 
number in units of h is the classical action of the system A[^](xi,X2) along the path x(r). 
Once we perform the integration along the path, this probability amplitude becomes clearly a 
function of the initial and final points in X space, Xi and X2, respectively. 

In this Feynman statistical procedure, the probability amplitude of the occurrence of any 
alternative of a set of independent alternatives is the sum of the corresponding probability 
amplitudes of the different independent events. The probability of the whole process is the 
square of the absolute value of the total probability amplitude. This produces the effect that 
the probability of the whole process can be less than the probability of any single alternative 
of the set. This is what Feynman calls interfering statistics. 

Then, the total probability amplitude that the system arrives at point X2 coming from Xi, 
i.e., Feynman's kernel K{xi,X2), is obtained as the sum or integration over all paths, of terms 
of the form of Eq. (13.111 . Feynman's kernel K{xi,X2), will be in general a function, or more 
precisely a distribution, on the X x X manifold. If information concerning the initial point is 
lost, and the final point is left arbitrary, say x, the kernel reduces to the probability amplitude 
for finding the system at point x, i.e., the usual interpretation of the quantum mechanical wave 
function $(x). By the above discussion we see that wave functions must be complex functions 
of the kinematical variables. 

We thus see that Feynman's quantization method enhances the role of the kinematical 
variables to describe the quantum state of an arbitrary system, in spite of the independent 
degrees of freedom. We consider that this is one of the reasons why the kinematical variables 
have to play a leading role also in the classical approach. 

We are used to consider in quantum mechanics, instead of a single function $(x), mul- 
ticomponent wave functions, i.e, a set of linearly independent functions ipiit^r) defined on 
space-time and labeled with a discrete subindex that runs over a finite range, such that it can 
be considered as a vector valued function in a finite dimensional complex vector space. In 
general this finite space carries some irreducible representation of the rotation group and each 



^ R.P. Feynman and A.R. Hibbs, Quantum Mechanics and Path Integrals, MacGraw Hill, NY (1965), p. 36. 
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component ijji represents a definite spin state of the system. Nevertheless, our wave function 
$(a;) depends on more variables than space-time variables. Once we define later the complete 
commuting set of observables to obtain, in terms of their simultaneous eigenvectors, an or- 
thonormal basis for the Hilbert space of states, we shall find that can be separated in two 
parts. One part r) depending on space-time variables and another part x that depends 
on the remaining translation invariant kinematical variables, that in our case will reduce to 
the velocity u and orientation a. It is this possible separation of our wave function that will 
produce the emergence of the different components of the usual formalism. 



3.1.1 Transformation of the wave function 

To see how the wave function transforms between inertial observers, and therefore to obtain 
its transformation equations under the kinematical groups, let us consider that O and O' are 
two inertial observers related by means of a transformation g E G, such that the kinematical 
variables transform as: 

x" = r{x,g). (3.2) 

If observer O considers that the system follows the path x(r), then it follows for O' the 
path x'{t) = f{x{T),g) and because the action along classical paths transforms according to 
Eq. ()1.28|) . the probability amplitude for observer O' is just 

(j)'[x'{T)] = AT exp I ^ r L{x\T),i'{T))dT 



h 



Tl 



j.e. 



iVexp <; - I L{x{r),^{r))dT exp <j - ^ da{g-^x{T)) 

(f)'[x'{T)] = (1>[x{t)] exp 1^ ia{g;x2) - a{g;xi)) 

where the last phase factor is independent of the integration path. If we add all probability 
amplitudes of this form, it turns out that Feynman's kernel transforms as: 

K'{x[,X2) = K{xi,X2)exp |^ {a{g;x2) - a{g;xi))^ . (3.3) 

If information concerning the initial point Xi is lost, the wave function transforms as the part 
related to the variables X2, up to an arbitrary function on G, 

$'(x'(x)) = 'i>'{gx) = $(x) exp |^ {a{g; x) + 9{g))^ , (3.4) 

or in terms of unprimed x variables 

$'(x) = ^g-'x) exp 1^ {a{g; g-'x) + 9ig)) | , (3.5) 

where 9{g) is some function defined on G but independent of x. 

Since our system is somewhere in X space, the probability of finding the system anywhere is 

1. Then we have to define the way of adding probabilities at different points x E X. If we define 
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a measure on X, such that dfi{x) is the volume element in X space and \^{x)\'^dfj,{x) is 

interpreted as the probability of finding the system inside the volume element dfi{x) around 
point X, the probability of finding it anywhere in X must be unity, so that 

/ \<l>{x)\^dfi{x) = 1. (3.6) 
Jx 

Since from ()3.5|1 

\<^'{x')\' = \<l>{x)\', (3.7) 

it is sufficient for the conservation of probability to assume that the measure to be defined yu(x) 
is group invariant. In that case, equation (j3.7|) implies also that inertial observers measure 
locally the same probability. This will have strong consequences about the possibility of invari- 
ance of the formalism under arbitrary changes of phase of the wave function. But the phase 
can be changed in a different manner at different points x. We can use this fact to further 
impose the local gauge invariance of the theory. It must be remarked that this arbitrary change 
of phase P{x) is not only a phase on space-time, but rather on the whole kinematical space 
of the system and this enlarges the possibilities of analyzing different transformation groups 
that can be more general than the original kinematical groups, because they act on a larger 
manifold. 

Consequently, the Hilbert space H whose unit rays represent the pure states of the system 
is the space of squared-integrable functions L^(X, yu) defined on the kinematical space X, /i(x) 
being an invariant measure such that the scalar product on H is defined as 

< $1^ >= / $*(a;)^(a;)d/i(x), (3.8) 
Jx 

$*(x) being the complex conjugate function of ^{x). There is an arbitrariness in the election 
of the invariant measure /i(x) but this will be guided by physical arguments. Nevertheless, the 
invariance condition will restrict the possible measures to be used. 



3.1.2 Representation of Observables 

Wigner's theorem, ^'^ implies that to every symmetry g E G of a continuous group, there 
exists a one to one mapping of unit rays into unit rays that is induced on ?i by a unitary 
operator U{g) defined up to a phase that maps a wave function defined on x into an arbitrary 
wave function of the image unit ray in x'. The Relativity Principle is a strong symmetry 
of physical systems that defines the equivalence between the set of inertial observers whose 
space-time measurements are related by means of a transformation of a kinematical group G. 
Now, if we interpret $(x) as the wave function that describes the state of the system for the 
observer O and $'(a;) for O', then we have 

U{gMx) = $'(x) = ^g-'x)exp Ua{g- g-'x) + e{g)] . (3.9) 



^ E.P. Wigner, Group theory and its application to the quantum mechanics of atomic spectra, Acad. Press, 
NY (1959). 

3 V. Bargmann, J. Math. Phys. 5, 862 (1964). 
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Since the 9{g) function gives rise to a constant phase we can neglect it and then take as 
the definition of the unitary representation of the group G on Hilbert space H 



*'W^t/te)*W^*(,-'.)e.p|,.te,-.)|. ,3.10) 

Gauge functions satisfy ()1.3np . and therefore the phase term can be replaced by 

"(S-; 9~^x) = -a{g'^] x) + a(0; x) + ^(^f, ^f"^) = -a{g'^] x) + ((5-), (3.11) 

because gauge functions can always be chosen such that a(0;x) = and the group function 
Ci.9) = i{.9i9~^) giving rise also to a constant phase, can be suppressed. We thus define the 
transformation of the wave function by 

^,„*W^*(,-..)exp{-i„fe-.;.)}. ,3.12) 

If the unitary operator is represented in terms of the corresponding self-adjoint generators of 
the Lie algebra in the form 

f/((7) =exp|--^(7'^X^|, (3.13) 

then, for an infinitesimal transformation of parameters dg'^ its inverse transformation has 
infinitesimal parameters —Sg'^, we obtain at first order in 6g'^ 



while 



and 



U{SgMx) = [l--6g''X,\ <l>(x) = $(x) - -6g''X,^x), 



H5g-'x) ^ ^fix,5g-')) = ^x) - 5g''ui{xf^^''^ 



dx' 



exp!^-^a{5g ^] x)^ = 1 - ^a{5g ^]x). 
But because a{0; x) = 0, 

a(i9-';x) = ^ (-V). 

Og" 9=0 

and the substitution of the above terms in (j3.12j) and further identification of the first order 
terms in Sg'^ imply that the self-adjoint operators X^- when acting on the wave functions have 
the differential representation 

X^ = ^uiix)^-v^ix), (3.14) 

where 



dg- 



9=0 



If we restrict ourselves to transformations of the enlarged configuration space (t, qi) that 
can be extended to the whole kinematical space x = {t,qi, . . . , ^ ^■*), then, using the same 
notation as in p.l2|) - p.l5L if the infinitesimal transformation is of the form 

t' = t + MJg-, q[ = q, + M.Jg'^, . . . , q'f^ = q^'^ + M^^dg'', 
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these generators take the form 




+ ... + MI 



la 



(k-1) 



dql 



d 



) 



(3.16) 



When compared with the Noether constants of the motion ()1.22|) written in the form 



we see a certain kind of 'correspondence recipe'. When restricted to kinematical groups, 
the functions B„{x) of ()1.22p . are obtained from the Lagrangian gauge functions a{g;x), by 
(ll.29jl . which is exactly the same derivation as the functions Va{x) above in ()3.15jl . Now, by 
identifying the different classical observables and generalized momenta that appear here in 
(I3.17jl with the corresponding differential operators of ()3.16|1 that multiply the corresponding 
MIJ function, we get: the generalized Hamiltonian H = p^^^g- — L, which is multiplied in 
(|3.17jl by the function M^, is identified with the operator ihd/dt which is also in front of 
the function M„ in (|3.16L and similarly, the generalized momentum the factor that 

multiplies the function M-^\ with the differential operator —ihd/dqf \ for s = 0, . . . , — 1, 
because the functions Va^x) = B^^^x), are the same. 

(s) 

Remember that pIs+i) are canonical conjugate variables. Then, each generalized 

momentum is replaced by (h/i) times the differential operator that differentiates with respect 
to its conjugate generalized coordinate and the generalized Hamiltonian by ihd/dt. 

The Heisenberg representation is that representation in which the time dependence has 
been withdrawn from the wave function by means of a time dependent unitary transformation. 
Then the wave function in this representation depends on the kinematical variables with the 
time excluded, i.e., it depends only on the generalized coordinates gl^"*. Therefore, when acting 
on the wave function in the Heisenberg representation ?/'(gj, . . . , gf'^"^'*), the observables 
qf"^ and p^^^ satisfy the canonical commutation relations 



If functions Va{x) in (I3.14jl vanish, the X^j generators satisfy the commutation relations of 
the group G. But if some v^^x) ^ the X„ generators do not satisfy in general the commutation 
relations of the initial group G where they come from, but rather the commutation relations of 
a central extension of G. The group representation is not a true representation but a projective 
representation of G as shown by Bargmann. ^ 

In fact, from (I3.1()jl we get 



-N^ = -HM^ + p\,^,)Mt^ - B^{x) 



(3.17) 



UigiMx) 



$(^1 ^x) exp{-a(^i; ^/x)} 



acting now on the left with U{g2), 



U{g2)U{gi)<l>{x) = t/(^2)$(^i ^x)exp{-a{gi;g^ V)} 

1 1 

= ^{{92gi)~^x) exp{-a{g2; g2^x)} exp{-a{gi; (^2^i)~V)}, 



(3.18) 



^ V. Bargmann, Ann. Math. 59, 1 (1954). 
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while acting on ^{x) with U{g2gi) 



U{g2gi)'^{x) = <i>{{g2gi) ^x) exp{-a{g2gi] {g2gi) ^x)}. 



(3.19) 



If we define {g2gi) ^x = g^ ^g2 ^x = z, then giz = g2 ^x and because gauge functions satisfy 
(jl.3()|] . we write 

a{g2; 9iz) + a{gi; z) = a{g2gi] z) + ^(5(3, 5-1), (3.20) 



and by comparing (|3.18ll with ()3.19ll . taking into account (13.2011 . we obtain 



U{g2)U{g,)^{x) = U{g2giMx) e^Uid^. Qi)}. (3.21) 



Since $(a;) is arbitrary, we have a projective unitary representation of the group G characterized 
by the non-trivial exponent ^{g,g'). 

For both Galilei and Poincare particles the kinematical space is the ten-dimensional mani- 
fold spanned by the variables (t, r, u,a.),t being the time, r the charge position, u the velocity 
and a the orientation of the particle. Thus in the quantum formalism the wave function of 
the most general elementary particle is a squared-integrable function $(t, r, u, ex.) of these 
kinematical variables. For point particles, the kinematical space is just the four-dimensional 
space-time, so that wave functions are only functions of time and position, but spinning par- 
ticles will have to depend on the additional variables like velocity and orientation. The spin 
structure will thus be related to these additional variables. 

3.1.3 Nonrelativistic spinning particles. Bosons 

Now let us apply the formalism to the most interesting case of spinning particles. Let us 
consider next Galilei particles with (anti)orbital spin. This corresponds for example to systems 
for which X = Q/S0{3) and thus the kinematical variables are time, position and velocity. A 
particular classical example is given in Chapter 2, Section IT2I bv the free Lagrangian 



with u = dr/dt. For the free particle, the center of mass q = r — k has a straight motion while 
the relative position vector k follows an elliptic trajectory with frequency 00 around its center 
of mass, being the spin related to this internal motion. It is expressed as S = —mk x dk/dt. 
The kinematical variables transform under Q in the form 



The wave functions are functions on X and thus functions of the variables {t,r,u). On this 
kinematical space the gauge function is the same as in (I2.33jl , where m defines again the mass 
of the system. Taking into account as in the previous example the correspondence recipe for 
the Hamiltonian H — > ihd/dt, the first generalized momentum = P ^ —iKV and the other 




(3.22) 




t(r) + b, 

R{a)r{T) + vt{r) + a, 
R{ct)u{T) + V. 



(3.23) 
(3.24) 
(3.25) 
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generalized momentum P2 = U ^ —ihVu-, the generators of the projective representation are 
given by 

H^ih—, P=-V, K ^mr -t-W - -Vu, (3.26) 
at I II 

h h 

J X - V + UX -Vu^ L + Z, (3.27) 

where V is the gradient operator with respect to = r variables and the gradient operator 
with respect to the q2 = u variables. It is important to stress that this representation of the 
generators is independent of the particular Lagrangian that describes the system. It depends 
only on the kinematical variables (t, r, u) and the usual Galilei gauge function. 

If we define q = r — k = {K + Pt)/m, it satisfies the commutation relations with P, 

[qi,Pj] = ihSij, 

which are the canonical commutation relations between the linear momentum and position for 
a point particle and therefore these canonical commutation relations between the total linear 
momentum and the center of mass position for a spinning particle are already contained in the 
commutation relations of the extended Lie algebra of the kinematical group. Therefore the 
quantum mechanical operator 

q = r- — V„, (3.28) 
im 

can be interpreted as the center of mass position operator. Discussion of other possibilities for 
the center of mass position operator can be found in the book by the author. 

In this representation, one Casimir operator is the internal energy H — P"^ /2m. We see 
that the spin operator is defined as usual 

S ^ J - —K xP^uxU + kxP^ux -Vu + — V„ x -V; 

m i im i 

written in terms of two non-commuting terms, it satisfies 

[S, S] = ihS, [J, S] = ihS, [S, P] = [S, H] = [S, K] = 0, 

i.e., it is an angular momentum operator, transforms like a vector under rotations and is 
invariant under space and time translations and under Galilei boosts, respectively. The second 
part of the spin operator is of order so that it produces a very small correction to the first 
Z part. 

Operators Z satisfy the commutation relations 

[z, z] = ihz, [J, z] = ihz, [z, P] = [z, H] = 0, 

[Z, K] = -ihU = -ffVu, 

i.e., Z is an angular momentum operator, transforms like a vector under rotations and is in- 
variant under space and time translations but not under Galilei boosts. It is usually considered 
as the quantum mechanical spin operator. 

We see however, that the angular momentum operator J is split into two commuting terms 
r X P and Z. They both commute with H, but the first one is not invariant imder space 
translations. The Z operators are angular momentum operators that only difl'erentiate the 
wave function with respect to the velocity variables, and consequently commute with H and 
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P, and although it is not the true Galilei invariant spin operator, we can find simultaneous 
eigenstates of the three commuting operators H — P'^ /2m, and Z3. Because the Z operators 
only aff'ect the wave function in its dependence on u variables, we can choose functions with 
the variables separated in the form $(t, r, u) = J2i i^iit^ so that 

{H - P^/2m)iJi{t, r) = EiJi{t, r), (3.29) 

Z\,{u) = z{z + l)h\i{u), (3.30) 

ZsXiiu) = m^hxi{u). (3.31) 

The space-time dependent wave function 'ipi{t,r), satisfies Schroedinger's equation and is un- 
coupled with the spin part xi'^)- 

Due to the structure of Z^ in terms of the u variables, which is that of an orbital angular 
momentum, the spin part of the wave function is of the form 

x{u) = f{u)Y:^^{e,<j>), (3.32) 

f{u) being an arbitrary function of the modulus of u and F™^(^,0) the spherical harmonics 
on the direction of u. 

For the center of mass observer, S = Z and both angular momentum operators are the 
same. But for an arbitrary observer, Z operators do not commute with the boosts generators so 
that its absolute value is not Galilei invariant, while S is. But the splitting of the wave function 
into a multiple-component function that reflects its spin structure is an intrinsic property that 
can be done in any frame. 

It turns out that if for an arbitrary observer Z is not the spin of the system, r x P is 
not the conserved orbital angular momentum, because r does not represent the position of the 
center of mass of the particle. 

When there is an interaction with an external electromagnetic field, equation ()3.29|) is 
satisfied for the mechanical parts Hm = H — ecj) and Pm = P — eA and we thus obtain the 
usual equation 

(^H-e(P- ^ ^,(t, r) = Ei^iit, r). (3.33) 

This formalism, when the classical spin is of orbital nature, does not lead to half integer 
spin values, and therefore, from the quantum mechanical point of view these particles can be 
used only as models for representing bosons. 

3.1.4 Nonrelativistic spinning particles. Fermions 

Other examples of nonrelativistic spinning particles are those which have orientation and 
thus angular velocity. For instance, if X = Q/M.^, being the subgroup {R^, +} of pure 
Galilei transformations, then the kinematical space is spanned by the variables (t, r, a.). This 
corresponds for instance to the Lagrangian system described by 

^=2 (a) +2" - 

The particle travels freely at constant velocity while it rotates with constant angular velocity 
uj. The classical spin is just S = Joj, and the center of charge and center of mass represent 
the same point. 
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To describe orientation we can think of the three orthogonal unit vectors Cj, i = 1,2,3 
linked to the body, similarly as in a rigid rotator. If initially they are taken parallel to the 
spatial Cartesian axis of the laboratory inertial frame, then their nine components considered 
by columns define an orthogonal rotation matrix Rij{cx) that describes the triad evolution with 
the initial condition Rij{t = 0) = 6ij. 

Now, kinematical variables t, r and p transform under Q in the form 

t'ir) = t{T) + b, (3.35) 

r'(r) = R{a)r{r)+vt{r) + a, (3.36) 
P'ir) = /^ + Mr) + /xxp(r) _ 
1 - p ■ p(r) 

On the corresponding Hilbert space, the Galilei generators are given by: 

H = ih—, P=-V, K = mr -t-V, (3.38) 

ot I I 

J = ^ r X V + ^ {Vp + p X Vp + p(p ■ Vp)} = L + W, (3.39) 

Vp being the gradient operator with respect to the p variables and in the p parameterization 
of the rotation group. 

The W part comes from the general group analysis. The group generators in this parametriza- 
tion Xi will be obtained from (|3.37|) and according to (|1.38p and ()1.4fl|) . They are obtained 

as 



X,; 



dp''' 



d 

dp'' ' 



that can be written in vector notation as 

X = Vp + pxVp + p{p- Vp) 
They satisfy the commutation relations 

[Xi, Xk] = —2eikiXi 

h 

and therefore operators Wk = ^^k, or in vector notation 

W = ^^{V, + pxV, + pip-Vp)}, (3.40) 

will satisfy the angular momentum commutation relations 

[W,W]=ihW. (3.41) 

In this way since L and W commute among each other, we also get [J, J] = ihJ. 

In this example the center of mass and center of charge are the same point, L = r x P 
is the orbital angular momentum associated to the center of mass motion and W = S is the 
spin operator. The spin operator commutes with H, P and K and the wave function can be 
separated as $(t,r,p) = '^i4'i{t,r)xi{p) leading to the equations 

(H - P^/2m)iJi{t, r) = Eij,{t, r), (3.42) 
S\i{p) = s{s + l)h\i{p), (3.43) 
S^Xi{p)=mshxi{p). (3.44) 
Bopp and Haag ^ succeeded in finding s = 1/2 solutions for the system of equations ()3.43p 
^ F. Bopp and R. Haag, Z. Naturforschg. 5a, 644 (1950). 
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and ()3.44|1 . They are called Wigner's functions. ^ Solutions of ()3.43p for arbitrary spin s are 
but a linear combination of the matrix elements of a (2s + 1) x (2s + 1) irreducible matrix 
representation of the rotation group as can be derived from the Peter- Weyl theorem on finite 
representations of compact groups. ^'^'^ We shall deal with the s = 1/2 functions in the 
Appendix Section 13. 2| where explicit expressions and a short introduction to the Peter- Weyl 
theorem, will be given. 

To describe fermions, the classical particles must necessarily have compact orientation 
variables as kinematical variables, otherwise no spin 1/2 values can be obtained when the 
classical spin is related only to the zitterbewegung. 

3.2 Appendix: Spinors 

In this section of mathematical content we shall review the main properties of spinors, in 
particular those connected with the possible representation of the wave function to describe 
spin 1/2 particles. We shall describe the representations in terms of eigenfunctions of the 
different commuting spin operators. But it must be remarked that in addition to the spin 
operators in the laboratory frame we also have spin operators in the body frame, because our 
general spinning particle has orientation, and therefore, a local Cartesian frame linked to its 
motion. This produces the result that for a spin 1/2 particle the wave function necessarily is 
a four-component object. 

The general wave function is a function of the ten kinematical variables, ^{t,r,u, p), and 
the spin part of the system related to the translation invariant kinematical variables u and p 
is 

S = uxU + W = Y + W, (3.45) 

where Y and W are given by 

Y = ux^\/u, W = ^{Vp + px Vp + p(p-Vp)}, (3.46) 

in the tan(a/2) representation of the rotation group, as has been deduced in previous sections. 
V„ and Vp are respectively the gradient operators with respect to u and p variables. These 
operators always commute with the H = ihd/dt and P = —iKV operators, and therefore 
they are translation invariant. This feature allows the separation of the general wave function 
in terms of space-time variables and velocity-orientation variables to describe the translation 
invariant properties of the system. 

The above spin operators satisfy the commutation relations 

[Y, Y] = ihY, [W, W] = ihW, [Y, W] = 0, (3.47) 

and thus 

[5, S] = ihS. 

Because we are describing the orientation of the particle by attaching to it a system of three 
unit vectors e,, whose orientation in space is described by variables p or ck, then, if at initial 

® L.C. Biedenharn and J.D. Louck, Angular Momentum in Quantum Physics. Theory and AppUcation, 
Cambridge U. P., Cambridge, England (1989). 

^ A.R. Edmonds, Angular Momentum in Quantum Mechanics, Princeton U. P., Princeton NJ (1957). 
* N. Ja. Vilenkin, Fonctions speciales et Theorie de la representation des groups, Dunod, Paris (1969). 
^ A.O. Barut and R. Raczka, Theory of group representations and applications, PWN, Warszawa (1980). 
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instant t — we choose the body axes coincident with the laboratory axes, the components 
of the unit vectors Cj at any time are 

{Cijj = Rji{ot) = SjiCosa + njni{l — cos a) — ejikfiksma, (3.48) 

in the normal parametrization and also in the p parametrization by 

{ei)j = Rji{p) = - P^)^ji + 2piPi - 2ejikPk), (3.49) 

where the Cartesian components of the rotation axis unit vector n are: 

ni = sin 6* cos 0, n2 = sin 6* sin 0, 713 = cos ^, (3.50) 

where 9 is the polar angle and the usual azimuth angle. Explicitly: 

eii — cos q; + sin^ cos^ 0(1 — cos a), 

ei2 = cos sin a + sin^^^ sin (6 cos 0(1 — cos a), 

ei3 = — sin 6' sin sin q; + sin cos cos 0(1 — cos a), 

621 = — cos^sino; + sin^ ^sin0cos0(l — cosa), 

622 = cosQ; + sin^^sin^0(l — cosa), 

623 = sin ^ cos sin a + sin ^ cos ^ sin 0(1 — cos a), 

631 = sin ^ sin sin a + sin ^ cos ^ cos 0(1 — cos a), 

632 = — sin ^ cos sin q; + sin 9 cos 9 sin 0(1 — cos a) , 

633 = coso; + cos^ ^(1 — cosa), 

in the oc = an, or normal parametrization of the rotation group. In the p = tan(Q;/2)n 
parametrization the body frame is 

en = {1 + pI-pI-pI)/{1 + p'), 
ei2 = (2pip2 + 2p3)/(l + p'), 
ei3 = (2pip3-2p2)/(l + p'), 

e2i = (2p2Pi-2p3)/(l + p'), 
e22 = {1-pI + pI-pI)/{1+p'), 

623 = (2p2P3 + 2pi)/(l + p'), 

esi = (2pip3 + 2p2)/(l + p'), 
e32 = (2p3P2-2pi)/(l + p2), 

633 = (1-P?-P^+P^)/(1+P'), 

where p'^ = pf + pi + pi = tan^(a/2). 

In addition to the different components of the spin operators Si, Yi and Wi in the laboratory 
frame, we also have another set of spin operators. They are the spin projections on the body 
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axes Cj, i.e., the operators Ri = Bi - S, = ei-Y and Ti = ei - W , respectively. In particular, 
spin operators Tj, collecting terms from ()3.49l) and (13.46(1 . take the expression 



n 



Ti = ^{ei)kWk = ———^^{{l-p'^)5ik + 2piPk-2ek^jpj) 



k=l 

d d 

^klrPr 



pk{p ■ Vp)^ , 



dpk dpr 

and after some tedious manipulations we reach the final result, written in vector notation as 

T = ^ {Vp - p X Vp + p(p • Vp)} . (3.51) 

We see, by inspection, that this result can also be obtained from the expression of W in (j3.46|) . 
just by replacing p by — p, followed by a global change of sign. This is because we describe the 
orientation of the particle by vector p in the laboratory frame from the active viewpoint, i.e., 
with the laboratory reference frame fixed. However, its orientation with respect to the body 
frame is described by the motion of the laboratory frame, whose orientation for the body is 
— p, and the global change of sign comes from the change from the active point of view to the 
passive one. This is the difference in the spin description in one frame or another. 
It satisfies the following commutation relations 

[T, T] = -ihT, [T, W] = 0. 

and in general all spin projections on the body frame Ri, Mi and T,, commute with all the spin 
projections on the laboratory frame Si, Yi and Wi. This is in agreement with the quantum 
mechanical uncertainty principle, because spin components with respect to different frames are 
compatible observables. 

To find eigenstates of the spin operator we have to solve equations of the form: 

S^xiu, p) = s{s + l)h^xiu, p), S^xiu, p) = mhxiu, p). 

But we also have the orientation of the particle, and therefore the spin projections on the 
body axes. These projections commute with 5^ and 5*3, and it is possible to choose another 
commuting spin operator, like the T3 operator, and therefore our wave function can be taken 
also as an eigenvector of T3, 

T3X{u,p) = nhx{u,p), 

so that the complete commuting set of operators that describe the spin structure must also 
include spin projections on the body axes. 
The spin squared operator is 

= Y^ + W^ + 2Y ■ W, (3.52) 

and we see from ()3.47|1 that is expressed as the sum of three commuting terms and its eigen- 
vectors can be obtained as the simultaneous eigenvectors of the three commuting operators on 
the right-hand side of ()3.52|) . Operators Y and W produce derivatives of the wave function 
with respect to u and p variables, separately. Thus, each x(**, p) can again be separated as 

X(n,p) = ^f/,(n)\/,(p), (3.53) 
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where the sum runs over a finite range, and where Uj{u) will be eigenfunctions of Y"^ and 
Vj{p) of W"^ , respectively. 

Functions Uj{u) are multiples of spherical harmonics defined on the orientation of the 
velocity vector w, because the Y operator has the structure of an orbital angular momentum 
in terms of the u variables, and thus its eigenvalues are integer numbers. The global factor 
left out is an arbitrary function depending on the absolute value of the velocity u. 

It turns out that to find the most general spinor is necessary to seek also solutions of the 
Vj{p) part, depending on the orientation variables. This goal will be achieved in the next 
section, where we consider the action of the rotation group on itself as a transformation group. 



3.2.1 Spinor representation on SU(2) 

We shall describe now in detail the orientation part of the general wave function, V{p). 
If there is no contribution to spin from the zitterbewegung part Y , the spin operator (13.4511 
reduces to the W operator given in (j3.46j) . To solve the corresponding eigenvalue equations 
we shall first represent the spin operators in spherical coordinates. 

If we represent vector p = tan(a/2)n = rn in spherical coordinates (r, 9, 0), with r = \p\ = 
tan(a/2) and 9 and the usual polar and azimuth angles, respectively, then unit vector n has 
the Cartesian components given in ((Hinni)- If from now on we take h = 1, the spin operators 
()3.46jl are represented by the differential operators 



d 

[1 + r^) sin 9 cos d) — 
or 



1 



- cos 9 cos — sin — — 



Wo 



1 



'1 + r ) sin 9 sin i 



d_ 

Or 



d 



cos 6* sin + cos i 



89 

d_ 
89 



sm I 



+ 



r sin 9 
cos 9 sin 



COS^ COS( 

sin^^ 

cosri 



sin 9 



r sin 9 



'1 + r ) cos 6* — 

or 



d sm9 d 



39 



d 



The Casimir operator of the rotation group W"^ is: 



1 + 



'l + r' 



I — 

dr" 



+ 



2(1 



r^) d 



dr 



cos 9 d 



sin' 9 ^ 



The up and down spin operators defined as usual by W± = Wi ± iW2 are 



~2i 



2i 



d 



l + r^)sin6'^ + 
(1 + r^)sin6' 



dr 
d 



cos 9 



ir 



dr 



cos 9 — ir 



d_ 

89 ~ 

d_ 

09 



r cos 9 — 
r sin 9 

r cos 9 + i 
r sin 9 



d 



d 



They satisfy the commutation relations 



d 



d 



[W+,WJ\=2W^. 
(WJ)*, where * means to take the complex 



We can check that {Wi)* = -Wi and W+ = 
conjugate of the corresponding operator. 

If F^{r,9,(j)) is an eigenfunction of W^ and W3, it satisfies the differential equations: 



W,Fr(r,9, 



mF^{r,9,. 
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To find solutions of the above system we know that we can proceed in the following way. Let 
us compute first the eigenfunctions of the form F/. Then operator W+ annihilates this state 
W+F^ = and by acting on this function with operator W- we can obtain the remaining 
eigenstates of the same irreducible representation characterized by parameter s and for 
—s<m<s. Then our task will be to obtain first the functions. 

Now, let us consider eigenfunctions that can be written in separate variables as Ff (r, 6, 0) 
A{r)B{e)C{(j)). Then 

W3A{r)B{e)C{(j)) = sA{r)B{e)C{<l)) 

gives rise to 

Sin 

(1 + r^) cosOA'BC AB'C + ABC' = 2isABC 

r 

where A' is the derivative of A and so on, and by dividing both sides by ABC we have 

2, n^'ir) sine B' (9) C"(0) 
^ ' A(r) r B{e) C{(f)) 

Now, the third term on the left-hand side must be a constant, because the remaining terms 
are functions independent of (p. Therefore, this term is written as C '{(/))/ C {(p) = ik and thus 
C(0) = e*^*^ up to an arbitrary constant factor. Since C(0 + 27r) = C(0) this implies that the 
constant k must be an integer. The other two functions satisfy 

r(l + r^) cos 9A'B - sin 9AB' + ir{k - 2s) AB = 0. (3.54) 

If there exist solutions with real functions A and B, then necessarily k = 2s so that the 
eigenvalue s can be any integer or half integer, and equation (I3.54jl can be separated in the 
form: 

.,A'(r) sm9 B'(9) 

rl + r^— ff = = p = constant, 3.55 

^ ^ A{r) cos 6* B{9) ^ ' v ^ 

where, up to constant factors, the general solution is 

1 + r2 

By acting on this solution F/ = A{r)B{9)C{(j))^ with W+, since W+F^ = 0, it gives: 

r(l + r2)sin2 9A'B + (sin 9 cos 9 + ir sin 9)AB' — 2s{ir cos 9 + 1)AB = 0. 

By dividing all terms by AB, taking into account ()3.55|1 . we get the condition {p — 2s) (1 + 
ircos6') = 0. Then there exist real solutions in separate variables whenever p = 2s = k. They 
are given, up to a constant factor, by 

F:ir,0,<P) = (y^) {sin9re'''t (3.56) 
For s = 1/2 and after the action of W_ we obtain the two orthogonal spinors 



^W=(t^) , B{9) = {sm9^P 
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that produce a two-dimensional representation of the rotation group. We can similarly check 
that = 0. 

By inspection of the structure of PV± operators, if we take the complex conjugate of expres- 
sion W+F^ = we get — = and therefore ~ so that taking the complex 
conjugate spinors of the above representation we obtain another pair of orthogonal s = 1/2 
spinors, 



1/2 
1/2 



r cos 9 — i 



-1/2 
1/2 



sin 6 e 



-i<f> 



The remaining representations for higher spins can thus be obtained by the same method, 
or by taking tensor products of the above two-dimensional representations. For instance, for 
s = 1 we can obtain the following three orthogonal representations. From ()3.56l) with s = 1 
and acting with the W- operator we get 



1/2V 



^?=(*;/2)(</') 



sin 9 e 



j20 



o'i-4> 



1 + r2 
r 

sin 9 (i + r cos ^ ) e" , 

1 + 

(i + r cos 6*)^ 



1 + r2 



that can also be obtained as the tensor product \1/ (g) 

If we work in the normal or canonical representation of the rotation group, where the 
parameters are a = an, this amounts to replacing the variable r = tan(Q;/2) in terms of 
parameter a and expressing the differential operator d/dr in terms of d/da, and then the spin 
operators are given by 



Wo = ^ 



2 sin 9 cos d) -— 



d_ 

da 



2i 



2 sin 9 sin 6 —— 



d_ 

da 



cos 9 cos(f) 
tan(a;/2) 

/ cos 6* sin</) 
V tan(a/2) 

1 



d 
d 

+ cos (p I — 



sm ( 



2i 



d 



89 
sin 9 



tan(a/2) sin^ 
cos^ sind) 



+ 



COS^ COS( 



sin 9 



COS( 



sin^^ tan(a/2)sin0 



2 cos 9 

da tan(a/2) 39 



d d 



1 



d 



1 



tan(a/2)9a 4sin2(a/2) {89'^ sin 9 89 sin^ 9 



cos 9 d 



+ 



1 



w. 



~2i 
~2i 



d 

2sin9— + 
oa 

d 

2sin9— + 
oa 



cos 9 



tan(a/2) 
cos6' 



+ i 



— I 



d_ _ 
89 ~ 

d_ 

89 



cos^tan(a/2) — A d 



tan(a/2) sin^ 

cos^^tan(a/2) + 2^ d 
tan(a/2) sin 6* 



d 



d 



tan(a/2) 

and the orthogonal spinors of the two two-dimensional representations can be written as 

a 



and 



^1/2 



^1/2 



i sin — sin 6' e**^, 



a n 

cos — + ism — cos 9, 



^1/2 



a n 

COS I sm — cos 9 

2 2 



-1/2 
1/2 



-i sin — sin 6* e 
2 



(3.57) 
(3.58) 
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We have mentioned that the different spinors are orthogonal. To endow the group manifold 
with a Hilbert space structure it is necessary to define a hermitian, definite positive, scalar 
product. The Jacobian matrix of variables p' in terms of variables p given in (j;i;^7L has the 
determinant 

^ [dp') (i-^-pr 

and thus the transformation of the volume element 



We also get from (fT37|l that 



(1 - p- py 



and then the measure 



-a p 



l+p'2)2 ^(l+^2)(l + p2)y (l_^.p)4 ^ (l+p2)2 



is in fact an invariant measure. 

In spherical coordinates it is written as 



sin 6 



drdOdcj) 



(l + r2)2 

and in the normal representation is 

sin^(a/2) ?,m.9dad6d(f). 

Since the rotation group is a double-connected group, the above measure must be defined on 
a simply connected manifold, i.e., on the universal covering group of S'0(3), which is SU{2). 
The SU{2) group manifold in the normal representation is given by the three-dimensional 
sphere of radius 27i and where points on the surface of this sphere represent a unique SU{2) 
element, namely the 2x2 unitary matrix —I. The normalized invariant measure becomes 

diiN{a,9,(f)) = -— sin^(a;/2) sin 6*^0; ^6*^0. (3.59) 
Therefore, the hermitian scalar product will be defined as 

<f\g>= da del d(l) f*{a,e,(f))g{a,e,(f))sm^{a/2)sme, (3.60) 

471" Jo Jo Jo 

where /* is the complex conjugate function of /. 

All the previous computed spinors are orthogonal vectors with respect to the group invariant 
measure ()3.59|1 . In particular, the normalized s = 1/2 spinors are those given in ()3.57p - ()3.58p . 
multiplied by \/2. 

The spin projection operators on the body axis Cj linked to the particle, are given in (|3.5H1 
in the p parametrization, and we have seen that they differ from the spin operators W only 
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in the change of p ^ — p, and a global change of sign. In the normal parametrization this 
corresponds to the change a — > —a, followed by a global change of sign. 
It can be checked as mentioned before, that 

[Ti,Tk] = -ieikiTi, (3.61) 

m,n]=0. (3.62) 

Since W'^ = we can find simultaneous eigenvectors of the operators W'^, W3 and T3, 
which will be denoted by Dmn{a.) in such a way that 

Since W3{a)Dmn{a) = mDmnia), by producing the change a — *• —a we get W3{—a)Dmn{—a) 
mDmn{—a) and the subsequent global change of sign it reduces to 

-Wsi~a)Dl:li-a) = T,{a)D^:^^{-a) = -mD^) (-a), 



so that the above spinors (|3.57p - ()3.58j) are also eigenvectors of T3. 

With this notation, the four normalized spinors, denoted by the corresponding eigenvalues 

|s, m, >, are 

$1 = 1 1/2, 1/2, 1/2 > = v^(cos(a/2) + zcosesin(a/2)), (3.63) 
$2 = 11/2,-1/2,1/2 > = iV2sin(a/2) sin^e-*"^, (3.64) 
$3 = 11/2,1/2,-1/2 > = iV2sin(«/2) sin^e^'^. (3.65) 
$4 = 11/2,-1/2,-1/2 >= V2(cos(a/2) -z cos ^sin(a/2)), (3.66) 

(3.67) 

They form an orthonormal set with respect to the normalized invariant measure (j3.59jl and with 
the scalar product defined in ()3.6fl|l . We can check that the lowering operators l^_$i = $2, 
= 0, W_^3 = $4, W_^^ = 0, and simmilarly T_$i = 0, T_$3 = $1, T_$2 = 0, and 
T_$4 = $2, and the corresponding up relations when acting with the rising operators W+ and 
T+, respectively. Remark that because the opposite sign in the commutation relations of the 
Tj operators, here the T± operate in the reverse direction. 

The important feature is that if the system has spin 1/2, although the s = 1/2 irreducible 
representations of the rotation group are two-dimensional, to describe the spin part of the 
wave function we need a function defined in the above four-dimensional complex Hilbert space, 
because to describe orientation we attach some local frame to the particle, and therefore in 
addition to the spin values in the laboratory frame we also have as additional observables 
the spin projections in the body axes, which can be included within the set of commuting 
operators. 

3.2.2 Matrix representation of internal observables 

The matrix representation of any observable A that acts on the orientation variables or in 
this internal four-dimensional space spanned by these spin 1/2 wave functions is obtained 
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as Aij —< >, i,j = 1,2,3,4. Once these four normalized basis vectors are fixed, when 

acting on the subspace they span, the differential operators Wi and Tj have the 4x4 block 
matrix representation 

where <t are the three Pauli matrices and I represents the 2x2 unit matrix. We have included 
Planck's constant into the angular momentum operators. 

If we similarly compute the matrix elements of the nine components of the unit vectors 
(ej)j, i,j = 1,2,3 we obtain the nine traceless hermitian matrices 

We can check that the % = S ■ Si = Si ■ S . We see that the different components of the unit 
vectors e^, in general do not commute. The eigenvalues of every component ej^, in this matrix 
representation of definite spin, are ±1/3. However, the matrix representation of the square of 
any component is (e^)^ = 1/3, so that the magnitude squared of each vector e? = 'Yliji^ijY — ^ 
when acting on these wave functions. The eigenvalues of the squared operator (cj^)^ are not 
the squared eigenvalues of Cj^-. This is because the function eij(^k does not belong in general 
to the same space spanned by the $fc, k = 1, ... ,4 although this space is invariant space for 
operators Wi and Tj. In fact, each function ejj$fc is a linear combination of a spin 1/2 and a 
spin 3/2 wave function. 

We do not understand why any component of a classical unit vector e-ij of a Cartessian 
frame, can have as eigenvalues ±1/3 in the quantum case and its square {^j-Y — II/3 instead 
of 1/9. 



3.2.3 Peter- Weyl theorem for compact groups 

The above spinors can also be obtained by making use of an important theorem for repre- 
sentations of compact groups, known as the Peter- Weyl theorem, ^'^ which is stated without 
proof that can be read in any of the mentioned references. 

Theorem.- Let D'^^\g) be a complete system of non-equivalent, unitary, irre- 
ducible representations of a compact group G, labeled by the parameter s. Let dg 
be the dimension of each representation and D^^^{g), 1 < i, j < dg the correspond- 
ing matrix elements. Then, the functions 

^/IgDlfig), l<i,j<ds 
N. Ja. Vilenkin, Fonctions speciales et Theorie de la representation des groupes, Dunod, Paris (1969), p. 

39. 

A.O. Barut and R. Raczka, Theory of group representations and applications, PWN-Polish Scientific Publishers, 
Warszawa (1980), p. 174. 

F. Peter and H. Weyl, Math. Ann. 7, 735 (1927). 
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form a complete orthonormal system on G, with respect to some normalized in- 
variant measure fiNid) defined on this group, i.e., 

/ ^/d,Dl^>{g) ^/d;Dl^,\g)df^M = S'''5.,J,i. (3.71) 
Jg 



That the set is complete means that every square integrable function defined on G, f{g), admits 
a series expansion, convergent in norm, in terms of the above orthogonal functions Dlj\g), in 
the form 



(s) 

where the coefficients, in general complex numbers obtained by 

a^= [ VdsD^^*ig)fig)df^M- 



G 



In our case SU{2), as a group manifold, is the simply connected three-dimensional sphere of 
radius 27r, with the normalized measure as seen before ()3.59ll . 

d^iy(a,9,(f)) = — - sin 9 sm( a/ 2)"^ dad9d(f>. 

In the normal parametrization, the two-dimensional representation of SU (2) corresponds 
to the eigenvalue s = 1/2 of S*^ and the matrix representation is given by 

L)(i/2)(^) ^ cos(a/2)I - i sin(a/2)(u ■ <t), 

i.e., 

p.(i/2)( N _ /cos(a/2) — i cos^ sin(Q;/2) — i sin ^ sin(a/2) e"*''^ \ 
~ \ -i sin 9 sin(a/2) e'^ cos(a/2) + i cos ^ sin(a/2) ) ' 

If we compare these four matrix components with the four orthogonal spinors given in 
(ESl-eni we see that 



In the three-dimensional representation of 5*0(3), considered as a representation of SU(2) 



Df^ {a) = 6ij cos a + UiUj{l - cos a) + eikjUk sin a = ej. 

we get another set of nine orthogonal functions. Multiplied by a/3 they form another orthonor- 
mal set orthogonal to the previous four spinors. It is a good exercise to check this orthogonality 
among these functions. 



3.2.4 General spinors 

In the case that the zitterbewegung content of the spin is not vanishing we can also obtain 
spin 1/2 wave- functions as the irreducible representations contained in the tensor product of 
integer and half-integer spin states coming from the U{u) and V{p) part of the general wave 
function (I3.53|l . 



3.2. APPENDIX: SPINORS 
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The total spin operator of the system is of the form 

S = uxU + W = Y + W, 

where Y = —iKVu and W is given in ()3.46|1 . Spin projections on the body axes, i.e., operators 
Ti = Bi ■ W, are described in ()3.5H) . They satisfy the commutation relations 

[Y, Y] = iY, [W, W] = iW, [T, T] = iT, 

[Y,W] = Q, [Y,T]=Q, [W,T]=Q. 

These commutation relations are invariant under the change p by —p in the definition of the 
operators W and T, because they are changed into each other. The expression of the body 
frame unit vectors Cj is given in (I3.48jl and (j3.49j) . 

We can see that these unit vector components and spin operators Wj and Tj satisfy the 
following properties: 

1) eij{-a, 9, (j)) = -ej.{a, 9, (p). 

2) e,-W = j:^eijW,=T,. 

3) E,e,T, = W. 

4) For all the action WiCj. = 0, with no addition on index i. 

5) For all the action TjCjj = 0, with no addition on index i. 

6) For all fc, with i ^ j, we have that WiCkj + WjCki = 0, and in the case that i = j, 
it leads to property 4. 

7) For all k, with i ^ j, we have that TjCj^ + Tjei^ = 0, and similarly as before in the 
case i = j it leads to property 4. 

This implies that Si ■ W = W ■ = T^, because of property 4, since when acting on an 
arbitray function /, 

{W ■ e,)/ ^ Yl ^Mjf) = / E ^^■(^^^■) + E ^ij^j(f) = Tiif), 
j j j 

because WjCij = 0. 

In the same way J2j ^j'^j = '^j^j = 

Now we fix the value of spin. Particles of different values of spin can be described. Let us 
consider systems that take the lowest admissible spin values. For spin 1/2 particles, if we take 
first for simplicity eigenfunctions V{p) of W"^ with eigenvalue 1/2, and then since the total 
spin has to be 1/2, the orbital Y part can only contribute with spherical harmonics of value 
y = and y = I. 

If there is no zitterbewegung spin, y = 0, and Wigner's functions can be taken as simulta- 
neous eigenfunctions of the three commuting W'^, W3, and T3 operators, and the normalized 
eigenvectors \w,W3,t3 > are explicitly given by the functions ()3. 63113. 66|) . 

If we have a zitterbewegung spin of value y = I, then the U{u) part contributes with the 
spherical harmonics 

Y,\9,4>) = |l,l>=-sin(0~)e^^yA, (3.73) 

Y,'{9,4>) ^ \l,0>=cos{9)\^, (3.74) 

V 47r 

Y,-\9,4^) ^ |l,-l>=sin(^~)e-*^yj, (3.75) 



86 



CHAPTER 3. QUANTIZATION OF THE MODELS 



normalized with respect to the measure 

sin 



JO 

2 



which are the indicated eigenfunctions \y,y3> oiY and F3, and where the variables 9 and 
determine the orientation of the velocity u. 

The tensor product representation of the rotation group constructed from the two irre- 
ducible representations 1 associated to the spherical harmonics (I3.73jl - (l3.75jl and 1/2 given in 
(ISinnil-dSinEI) is split into the direct sum 1 ® 1/2 = 3/2 © 1/2. 

The following functions of five variables 9, 0, a, 9 and 0, where variables 9 and correspond 
to the ones of the spherical harmonics Yj^, and the remaining a, 9 and 0, to the previous spinors 
$j, are normalized spin 1/2 functions |s, 53,^3 > that are eigenvectors of total spin S^, and 5*3 
and T3 operators 

^ (f°<I>i - V2Y,'^2) , (3.76) 

(3.77) 

(r0$3 - v^n^-f 4) , (3.78) 

i=(-ro$, + V2Ffi<|.3), (3.79) 

such that \l/2 = S-'^i and similarly = 5'_\E'3, and also that ^^3 = T_\E'i, and = T_\&2- 
They are no longer eigenfunctions of the W3 operator, although they span an invariant vector 
space for S'^, S3 and T3 operators. In the above basis (I3.76|] - ()3.79jl formed by orthonormal 
vectors "^i, the matrix representation of the spin is 

^ = ^ + ^=?fn !V (3.80) 



^1 = 


|l/2,l/2,l/2> 


^2 = 


11/2,-1/2,1/2 > 


^3 = 


|l/2,l/2,-l/2> 


^4 = 


11/2,-1/2,-1/2 > 



(-F°$2 + v^rr^<f 1) , 



2 V cr 

while the matrix representation of the Y and W part is 



y-'4(: w^4(: !i. (3.81) 



3V0'^/' 6V0'^ 
which do not satisfy commutation relations of angular momentum operators because the vector 
space spanned by the above basis is not an invariant space for these operators Y and W . 

The spin projection of the W part on the body axis, i.e., the T operator, takes the same 
form as before (13.691) 



T.Ail n, T^^Ai ° T:.-'Jl ° I. (3.82) 



2vioy' " 2\-ii Q J ' 2vo-n 

because and ^^2 functions are eigenfunctions of T3 with eigenvalue 1/2, while \E'3 and \E'4 are 
of eigenvalue —1/2, and thus the spinors span an invariant space for Si and T, operators. 
In fact the basis is formed by simultaneous eigenfunctions of total spin S*^, 6*3 and T3, and the 
ket representation is the same as in the case of the given in (I3.63j) - (|3.66jl . 

The expression in this basis of the components of the unit vectors are represented by 



Chapter 4 



Dirac equation and analysis of Dirac 
algebra 



4.1 Quantization of the u = c model 



For Luxons we have the nine-dimensional homogeneous space of the Poincare group, spanned 
by the ten variables (t, r, u, a.) similarly as before, but now u is restricted to m = c. For this 
system, since u ■ u = and n 7^ 0, we are describing particles with a circular internal orbital 
motion at the constant speed c. 

In the center of mass frame, (see Fig l4.1|] the center of charge describes a circle of radius 
Rq = S/mc at the constant speed c, the spin being orthogonal to the charge trajectory plane 
and a constant of the motion in this frame. Let us consider the quantization of this u = c 
model whose dynamical equation is given by (I2.126I1 . 

If we analyse this particle in the centre of mass frame it becomes a system of three degrees 
of freedom. These are the x and y coordinates of the point charge on the plane and the phase 
a of the rotation of the body axis with angular velocity u. But this phase is the same as the 
phase of the orbital motion, as we shall see later, and because this motion is a circle of constant 
radius only one degree of freedom is left, for instance the x coordinate. In the centre of mass 
frame the system is equivalent to a one-dimensional harmonic oscillator of angular frequency 
u = mc^/S in its ground state. 

Identification of the ground energy of the one-dimensional harmonic oscillator huj/2 with 
the rest energy of the system in the center of mass frame +mc^, for if > particles, implies 
that the classical constant parameter S = h/2. All Lagrangian systems defined with this 
kinematical space, irrespective of the particular Lagrangian we choose, have this behaviour 
and represent spin 1/2 particles when quantized. 



4.2 Dirac' 



s equation 



The kinematical variables of this system transform under V according to 



7t(r) + 7(1; ■ R{fj,)r{r))/c^ + b, 



(4.1) 





{v ■ R{fj.)r{T))v + a, 



(4.2) 



87 



CHAPTER 4. DIRAC EQUATION AND ANALYSIS OF DIRAC ALGEBRA 



R{h)u{t) +-fv + {v- R{fi)u{T))v-f^/{l + 7) 



,2 



// X _ M + p(r) + /X X p(r) + /x; -a(r), p(r)) 

l-M-p(r) + G>,/x;n(r),p(r)) ' 




Figure 4.1: Motion of the charge in the CM. frame. 

where the functions and Gc are given in ()2.92|] and ()2.93|] . respectively. When quantized, 
the wave function of the system is a function r, u, p) of these kinematical variables. For 
the Poincare group all exponents and thus all gauge functions on homogeneous spaces are 
equivalent to zero, and the Lagrangians for free particles can thus be taken strictly invariant. 
Projective representations reduce to true representations so that the ten generators on the 
Hilbert space, taking into account ()4.H) - ()4.4p and ()3.15|) are given by: 

H = ^h^. P = -V, K = r^|-t-V-45xn, (4.5) 

h 

J = rx-V + S, (4.6) 

where as we shall see, the angular momentum operator S represents Dirac's spin operator and 
is given by the differential operator 

h h 

S = ux - Vu + j^{Vp + P^'^p + p{p- Vp)} = uxU + W = Su + S^, (4.7) 

and where the differential operators V„ and Vp are the corresponding gradient operators with 
respect to the u and p variables as in the Galilei case. The operator S is not a constant of the 
motion even for the free particle, and although it is not the angular momentum of the system 
with respect to its center of mass we keep this notation because it is the equivalent to Dirac's 
spin operator. Of course, it reduces to the true spin S in the center of mass frame. 
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To obtain the complete commuting set of observables we start with the Casimir invariant 
operator, or Klein-Gordon operator 

In the above representation, H and P only differentiate the wave function with respect to time 
t and position r, respectively. Since the spin operator S operates only on the velocity and 
orientation variables, it commutes with the Klein-Gordon operator ()4.8|l . Thus, we can find 
simultaneous eigenfunctions of ()4.8I1 . S^, and S3. This allows us to try solutions in separate 
variables so that the wave function can be written 

r, u,p) = J2 ^i(^> 'r)xiiu, p), (4.9) 

i 

where ipi{t,r) are the space-time components and the Xii'^^P) represent the internal spin 
structure. Consequently 

{H"^ - c^P"^ -m^c^)^i{t,r) = (4.10) 

i.e., space-time components satisfy the Klein-Gordon equation, while the internal structure 
part satisfies 

S\,{u, p) = sis + l)h''xi{u, p), (4.11) 

SsXiiu, p) = mshxiiu, p). (4.12) 

Eigenfunctions of the above type have been found in Section lX2| in particular we are interested 
in solutions that give rise to spin 1/2 particles. These solutions, which are also eigenvectors of 
the spin projection on the body axis T3, become a four-component wave function. 

For spin 1/2 particles, if we take first for simplicity eigenfunctions x(p) of S"^ with eigenvalue 
1/2, then since the total spin has to be 1/2, the orbital zitterbewegung part Y = u x U can 
only contribute with spherical harmonics of value 2/ = and y = 1. This means that we can 
find at least two different kinds of elementary particles of spin 1/2, one characterized by the 
singlet ?/ = (lepton?) and another hy y = 1 (quark?) in three possible states according 
to the component y^. If we call to the spin part Y the colour, we can have colourless and 
coloured systems of spin 1/2. The three different colours y^ are unobservable because the 
states ()3. 76113. 79|l are eigenstates of 5*3 and T3 but not eigenstates of Y^,. Nevertheless this 
interpretation of this spin part as representing the colour, as in the standard model, is still 
unclear and will be discussed in more detail in the forthcoming workshop Spin05 next week. 

For y = 0, the spin 1/2 functions Xi(p) are linear combinations of the four $j functions 
()3.63|) - ()3.66ll and in the case y = 1 they are linear combinations of the four of ()3.76l) - ()3.79|) . 
such that the factor function in front of the spherical harmonics is 1 because for this model 
M = c is a constant. It turns out that the Hilbert space that describes the internal structure 
of this particle is isomorphic to the four-dimensional Hilbert space C^. 

If we have two arbitrary directions in space characterized by the unit vectors u and v 
respectively, and Su and Sv are the corresponding angular momentum projections Su = u - S 
and Sv = V ■ S, then S^u = ~Su, and [Su, Sv] = ihSuxv- In the case of the opposite sign 
commutation relations of operators Tj, we have for instance for the spin projections [Ti,T2] = 
—ihT^, thus suggesting that ei x 62 = —63, and any cyclic permutation 1 — 2 — > 3, and thus Cj 
vectors linked to the body, not only have as eigenvalues ±1/3, but also behave in the quantum 
case as a left-handed system. In this case vectors are not arbitrary vectors in space, but 
rather vectors linked to the rotating body and thus they are not compatible observables, so 
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that any measurement to determine, say the components of e,, will produce some interaction 
with the body that will mask the measurement of the others. We shall use this interpretation 
of a left-handed system for particles later, when we analyse the chirality in section I1.2.2I For 
antiparticles it will behave as a right handed one. 

Operators Si and Tj have the matrix representation obtained before which is just 

where we represent by cr the three Pauli matrices and I is the 2x2 unit matrix. 

Similarly, the matrix elements of the nine components of the unit vectors (ej)j, i,j = 1, 2, 3 
give rise to the two alternative sets of representations depending on whether the zitterbewegung 
contribution is y = or ?/ = 1. In the first case we get 

while in the y = 1 case the representation is 

"^ = -9U oj'^^ = -9i-.^ 0j'^-^ = -9U -^)- ^'-''^ 

It must be remarked that the different components of the observables are not compatible in 
general, because they are represented by non-commuting operators. 

We finally write the wave function for spin 1/2 particles in the following form for y = 

$(o)(t, r, u,a) = Y, r)<!>i{a, 9, 0), (4.17) 

i=l 

independent of the u variables, and in the case y = 1 hy 

i=A 

$(1) (t, r, n, a) = ^ V^,(t, r)^,(^, 0; a, 9, <j)). (4.18) 

i=l 

where 9 and (p represent the direction of vector u. Then, once the $j or functions that 
describe the internal structure are identified with the four orthogonal unit vectors of the inter- 
nal Hilbert space C^, the wave function becomes a four-component space-time wave function, 
and the six spin components Si and Tj and the nine vector components (ej)j, together the 
4x4 unit matrix, completely exhaust the 16 linearly independent 4x4 hermitian matrices. 
They form a vector basis of Dirac's algebra, such that any other translation invariant internal 
observable that describes internal structure, for instance internal velocity and acceleration, 
angular velocity, etc., must necessarily be expressed as a real linear combination of the men- 
tioned 16 hermitian matrices. We shall see in Sec. 14.31 that the internal orientation completely 
characterizes its internal structure. 

The spin operator S = u x U + W which, as seen in ()3.8()ll and (|4.13jl . coincides with the 
usual matrix representation of Dirac's spin operator. 
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If we consider the expression of the kinematical momentum for u = c particles 



K = —r~tP 



c 



-S X u 



and we take the time derivative of this expression followed by the scalar product with u, it 
leads to the Poincare invariant operator (Dirac's operator): 



When Dirac's operator acts on a general wave function $(0) or we know that H and P 
have the differential representation given by (I3.38jl and the spin the differential representation 
()3.46L or the equivalent matrix representation (I4.13|] . but we do not know how to represent 
the action of the velocity u and the (du/dt) x u observable. However, we know that for this 
particle u and du/dt are orthogonal vectors and together with vector u x du/dt they form an 
orthogonal right-handed system, and in the center of mass frame the particle describes a circle 
of radius Rq = h/2mc in the plane spanned by u and du/dt. 



Figure 4.2: Representation of the local body frame and the different observables for the 
(a) H > solution and (b) H < solution. This orientation produces Dirac equation in 
the Pauli-Dirac representation 

Let us consider first the case y = 0. Since u and du/dt are translation invariant observables 
they will be elements of Dirac's algebra, and it turns out that we can relate these three vectors 
with the left-handed orthogonal system formed by vectors Ci, 62 and 63 with representation 
()4.15|] . Then, as shown in part (a) of Figure for the H > system, we have u = asi and 
du/dt X u = bes, where a and b are constant positive real numbers. Then the third term in 
Dirac's operator is {b/c^)e^ ■ S = {b/c^)T3, and ()4.19p operator becomes 




(4.19) 




(b) 



H-aP-ei--n = 0. 



(4.20) 
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If we make the identification with the H < solution of part (6) of Figure 14.21 the relation 
of the above observables is opposite to the previous one but now with the coefficients —a and 
—6, respectively, i.e., we get 

H + aP-ei + \n = 0, (4.21) 

which clearly corresponds to the change H —H in equation (|4.2np . 

Multiplying ()4.2ip by ()4.2n|l we obtain an expression whihch is satisfied by both particle 
and antiparticle 

ij2_^p2j_iiLi = o, (4.22) 

9 4c4 ' ^ ' 

and which is an algebraic relation between and P^. By identification of this expression with 
the Klein-Gordon operator ()4.8|) . which also contains both H > and H < solutions, leads 
to a = 3c and b = 2mc^/h = c^/Rq and by substitution in ()4.2np we obtain Dirac's equation: 

H -cP-a- pmc^ = 0, (4.23) 

where Dirac's matrices a and P are represented by 

and thus Dirac's gamma matrices are 

/./J=f' °), ^.o»a=f ° Z), (4.25) 



^0 -ly ' ' - ' " y-a ^ 

i.e., Pauli-Dirac representation, where 3ei plays the role of a unit vector in the direction of the 
velocity. Substitution into (|4.2H) corresponds to the equivalent representation with the change 

This representation is compatible with the acceleration du/dt lying along the vector 62. In 
fact, in the center of mass frame and in the Heisenberg representation, Dirac's Hamiltonian 
reduces to H = Pmc^, and the time derivative of any observable A is obtained as 

dA i dA 

such that for the velocity operator n = co:, 

du i r 2n n 2mc^ f icr\ ^ / . 

- = ^K^/3,col = — (_._^ ,}=j^3e,, (4.27) 



c^/Ro being the constant modulus of the acceleration in this frame, and where 3e2 plays the 
role of a unit vector along that direction. 

The time derivative of this Cartesian system is 

de-2 i ^ , c , , 

_ = _^3,„e^eJ = --e„ (4.29) 

^ = ^[8mc=,e3| = 0, (4.30) 
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since 63 is orthogonal to the trajectory plane and does not change, and where c/Rq = u is the 
angular velocity of the internal orbital motion. This time evolution of the observables e, is 
the correct one if assumed to be a rotating left-handed system of vectors as shown in Figure 
I4.2l -(a). It is for this reason that we considered at the beginning of this chapter that the body 
frame rotates with the same angular velocity as the orbital motion of the charge. 

To be consistent with the above consideration as 3ej as unit vectors, this means that the 
spin in the center of mass frame should be along 363. This is the case for the upper components 
while for the lower components (which in this representation correspond to < states) the 
orientation is the opposite. This means that for particles the corresponding set of axis forms 
a left handed system while for antiparticles they behave as a right handed system, showing a 
clear chirality difference between particles and antiparticles. 

In general 

— = -\H, S] = -\cP ■ a + (3mc^, S] = cP x a = P x u, 
dt n n 

is not a constant of the motion, but for the center of mass observer, this spin operator uxU+W 
reduces to the equivalent of the classical spin of the particle S and is constant in this frame: 

f 4[/3»c^S1^0. (4.31) 

Only the T3 spin component on the body axis remains constant while the other two Ti and 
T2 change because of the rotation of the corresponding axis. 

When analyzed from the point of view of an arbitrary observer, the classical motion is a 
helix and the acceleration is not of constant modulus c^/Rq, and the spin operator S is no 
longer a constant of the motion, because it is the total angular momentum J = r x P + S 
that is conserved. 

Identification of the internal variables with different real linear combinations of the 
matrices lead to different equivalent representations of Dirac's matrices, and thus to different 
expressions of Dirac's equation. 

For instance if we make the identification suggested by Figure 14.31 u = —ae^ and the 
observable du/dt x u = bei with positive constants a and b, we obtain by the same method 









f-,7 0\ 




or 




■. a 



/^=U n)' «=( n (4-35) 

and thus gamma matrices 



, _ a _ / I N 10 <T 



/-/3=l, . 7-7"«=l „ 1. (4.36) 



i.e., Weyl's representation. 
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u 




Figure 4.3: Orientation in the Weyl representation. 

When we compare both representations, we see that Weyl's representation is obtained 
from Pauli-Dirac representation if we rotate the body frame 7r/2 around 62 axis. Then the 
corresponding rotation operator 

R{n/2,e,) = exp(^^e2 • S) = ^M'^T,) = -L (^J "^^ . 

We can check that R'jpjjK^ = 7(^, where and '^y^r are gamma matrices in the Pauli- 
Dirac and Weyl representation, respectively. 

We can similarly obtain Dirac's equation in the case of zitterbewegung y = 1, by using the 
set of matrices (|4.16p instead of (|4.15p . because they are multiples of each other and only some 
intermediate constant factor will change. 

4.2.1 PCT Invar iance 

In Figure 14.41 we represent the parity reversal of the description of the electron as given 
by this model of luxon which is circling around the center of mass at the velocity c and that 
under P and in the center of mass frame it changes according to 

P : {r ^ -r,u^ -u, du/dt -du/dt, S ^ S,H ^ H}. 

In the Pauli-Dirac representation as we see in Figure 14. 2| this amounts to a rotation of value 
TT around axis 63 and thus 

P = R{n, 63) = exp(z7re3 ■ S/h) = exp^iTrT^/h) = i'jo, 

which is one of the possible representations of the parity operator ±70 or ±270. In Weyl's 
representation this is a rotation of value n around ei which gives again P = i'jQ. 
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Figure 4.4: Space reversal of the electron in the center of mass frame is equivalent to a 
rotation of value tt along S. 




Figure 4.5: Time reversal of the electron produces a particle of negative energy. 
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In Figure 14.51 we represent its time reversal also in the center of mass frame 

T : {r ^ r,u ^ -u, du/dt du/dt, S ^ S,H ^ -H}, 

but this corresponds to a particle of of if < such that the relative orientation of spin, velocity 
and position, given by equation (I2.127jl agrees with the motion depicted in this figure. 

A Dirac particle is a mechanical system whose intrinsic attributes are mass m > and 
spin h/2. We also see that the sign of H is also Poincare invariant and it is also an intrinsic 
property which establishes two different systems of the same value of m and S. The system 
with H > is called the particle and the other with H < the antiparticle. The value of 
the mass attribute is introduced by hand. To characterize its interaction with an external 
electromagnetic field, we also introduce by hand another intrinsic property the electric charge 
e, located at the point r. This implies that in addition to the mechanical properties m and S 
the system has as electromagnetic properties the electric charge e and because of the charge 
location separated from its center of mass and its motion at the speed of light, an electric 
dipole moment d and a magnetic moment /x, respectively. The electric charge can also have 
either a positive or negative sign. 



m 
H 
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The PCT transformation transforms particle into antiparticle and conversely, while keeping 
invariant the mechanical attributes m and S and the electromagnetic attributes n and d. The 
PCT invariance of the system establishes a relationship between the sign of H and the sign 
of e, although an indeterminacy exists in the election of the sign of the charge of the particle. 
The product eH is PCT invariant. 

This implies that particle and antiparticle have a magnetic moment and an oscillating 
electric dipole in a plane orthogonal to the spin. Once the spin direction is fixed, the magnetic 
moment of both have the same relative orientation with the spin, either parallel or antiparallel, 
according to the election of the sign of the electric charge. The electric dipole moment oscillates 
leftwards for particles and rightwards for antiparticles which shows a difference between them 
which is called chirality. If as usual we call the electron to the system of negative electric 
charge the particle, the above PCT transformation transforms the system (a) of figure Hl^ into 
the system (6). If what we call the particle is of positive electric charge, then the spin and 
magnetic moment are opposite to each other for both particle and antiparticle. 

However, to our knowledge no explicit direct measurement of the relative orientation be- 
tween spin and magnetic moment of the free electron, can be found in the literature although 
very high precision experiments are performed to obtain the absolute value of g, the gyromag- 
netic ratio. 

A plausible indirect experiment ^ has been proposed to measure the relative orientation 
between spin and magnetic moment for one outer electron atoms like Rb or Cs. 



^M.Rivas, Are the electron 
|ArXiv:physics/0112057, 



spin and magnetic moment parallel or antiparallel vectors?, LANL 
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a) 



b) 



Figure 4.6: Electromagnetic attributes fi and d for (a) a negatively charged particle and 
its PCT transformed (6), and their relative orientation with the spin, in the center of 
mass frame. The electric dipole of the particle oscillates leftwards and rightwards for the 
antiparticle. 

Rb*'' atoms have one electron at the level 5s. Its nucleus has spin 3/2 and the ground state 
of the atom has a total spin 1, and therefore the outer electron has its spin in the opposite 
direction to the spin of the nucleus. The magnetic moment of the atom is basically the magnetic 
moment of this outer electron because the inner shells are full and the magnetic moment of 
the nucleus is relatively smaller. 

Ultracold Rb®^ atoms in an external magnetic field will be oriented with their magnetic 
moments pointing along the field direction. If in this direction we send a beam of circularly 
polarized photons of sufficient energy ~ 6.8GHz to produce the corresponding hyperfine tran- 
sition to fiip the electron spin in the opposite direction and thus leaving the atom in a spin 
2 state, only those photons with the spin opposite to the spin of the outer electron will be 
absorbed. Measuring the spin orientation of the circularly polarized beam will give us the spin 
orientation of the electron thus showing its relationship with the magnetic moment orientation. 
Now the task is to check also the relative orientation for positrons. 



The classical model which satisfies Dirac's equation when quantized gives rise to two possi- 
ble physical systems of if > and H < 0. The if > is usually called the particle. According 
to the previous analysis the internal motion of the charge takes place on a plane orthogonal to 
the spin direction and in a leftward sense when we fix as positive the spin direction. For the 
antiparticle the motion is rightwards. For particles, the local orientable frame of unit vectors 
Si behaves as a left handed system rotating with an angular velocity in the opposite direction 
to the spin, while for antiparticles it can be considered as a right handed one. 

If we should have started the analysis by considering first the antiparticle, then in order 



4.2.2 Chirality 
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Figure 4.7: Relative orientation of the body axis for the antiparticle that leads to Pauli- 
Dirac representation. It behaves as a rotating right handed Cartessian frame around the 
spin direction. 



to get the same Pauli-Dirac representation as before we have to consider the body axis as the 
ones depicted in figure H771 i.e., in the opposite direction to the ones we chose before and this 
leads by the same arguments that the 7^ matrices have to replaced by the —7^, so that the 
Hamiltonian in the center of mass frame is —(3mc^. In this way the motion of the body frame, 
instead of ()4.28ll4.:'i();i is 



dei 

~dt 
de2 

"dt 

dt 



-[-(3mc^, 



-uje2, 



[-(3mc^, 62] = cjei, 



-[-/3mc^ 



0, 



(4.37) 
(4.38) 
(4.39) 



with 00 = c/Rq, which clearly corresponds to a rotating right handed system with an angular 
velocity around the spin direction. 

Matter is left and antimatter is right in this kind of models as far as the charge motion and 
the rotation of the local body frame are concerned, so that particles and antiparticles show a 
clear chirality. 

Although the local motion of the charge, which takes place in a region of order of Compton's 
wavelength, is probably physically unobservable, this motion corresponds nevertheless to the 
oscillation of the instantaneous electric dipole moment, which oscillates at very high frequency, 
but its sense of motion, once the spin direction is fixed, refiects this difference between particle 
and antiparticle. This electric dipole motion is independent of whether the particle is positively 
or negatively charged. 

Finally, when we compare the spin operator and the vector 63 we see 

2 V o- J ' ~ 3 I -0 



S=- 
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that the two upper components of the Dirac spinor correspond to positive energy solutions 
and therefore the upper components of these operators are related by S ~ 63, while the lower 
components correspond to negative energy solutions and for this components these operators 
behave as S ~ —63, a vector relationship which is clearly depicted in the figures and HTTI 
respectively. 



4.3 Dirac's algebra 

The three spatial spin components Si, the three spin projections on the body frame Tj and 
the nine components of the body frame (ej)j, i,j = 1,2,3, whose matrix representations are 
given in the y = case in (14 .1511 or in (|4.16jl in the y = I case, together with the 4x4 unit 
matrix I, form a set of 16 linearly independent hermitian matrices. They are a linear basis of 
Dirac's algebra, and satisfy the following commutation relations: 

[Si,Sj] = iheijkSk, [Ti,Tj] = -iheijkTk, [Si,Tj] = 0, (4.40) 

[Si, {ej)k] = iheikr{ej)r, [Ti, {ej)k] = -iheijr{er)k, (4.41) 
and the scaled 3ej vectors in the y = case 

4i 

[(3ei)fc, {Sej)i] = — {SijekirSr - 4«eiir^r) , (4.42) 

showing that the operators transform like vectors under rotations but they are not commut- 
ing observables. In the case y = I, the scaled — 9ej, satisfy the same relations. 

If we fix the pair of indexes i, and j, then the set of four operators 5*^, Si, Tj and (ej)j 
form a complete commuting set. In fact, the wave functions ^j, i = 1, ...,4, given before 
()3.63|) - ()3.66|] . are simultaneous eigenfunctions of S^, S3, T3 and (63)3 with eigenvalues s = 1/2 
and for S3, t3, and 633 are the following ones: 

$1 = |l/2, 1/2, 1/3 >, $2=1- 1/2, 1/2, -1/3 >, (4.43) 

$3 = |l/2, -1/2, -1/3 >, $4=1- 1/2, -1/2, 1/3 >, (4.44) 



and similarly for the spinors of ()3.76|) - ()3.79|l 

^1 = |l/2, 1/2, -1/9 >, ^2=1- 1/2, 1/2, 1/9 >, (4.45) 

^3 = |l/2, -1/2, 1/9 >, ^4=1- 1/2, -1/2, -1/9 > . (4.46) 
The basic observables satisfy the following anticommutation relations: 

{S„S,} = {T„T,} = ^6,,I, (4.47) 



{5„T,} = -(3e,)„ (4.48) 



{S„ (3e,),} = 2 6,kT„ {Ti, {Se,)^} = 2 5,,Sk, (4.49) 
2 2 

{{ei)j, {ek)i} = - Sik6jil + -tikrejis{er)s- (4.50) 
y o 
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If we define the dimensionless normalized matrices: 

2 2 

= 3{ei)j, {or aij = ~9{ei)j), Si = -Si, U = -Ti, (4.51) 

together with the 4x4 unit matrix I, they form a set of 16 matrices Fa, A = 1, . . . , 16 that are 
hermitian, unitary, linearly independent and of unit determinant. They are the orthonormal 
basis of the corresponding Dirac's Clifford algebra. 

The set of 64 unitary matrices of determinant +1, ±Tx, ±iTx, A = 1, . . . , 16 form a finite 
subgroup of SU{4:). Its composition law can be obtained from: 

Oij Clkl = ^ik^jl^ + i^ik^jlr Sr — iSjltikr tr + ^ikr^jls O-rs, (4.52) 

aijSk = iejkiau + 6jkti, (4.53) 

aijtk = -ieikiaij + 6ikSj, (4.54) 

Sittjk = ieikiaji + 6iktj, (4.55) 

SiSj = ieijkSk + SijI, (4.56) 

Si tj tj Si G'ji') (^4.5T^ 

tiCijk = -if^iji^ik + SijSk, (4.58) 
titj = -ieijktk + SijI, (4.59) 

and similarly we can use these expressions to derive the commutation and anticommutation 
relations ()4. 40114. 50|) . 

Dirac's algebra is generated by the four Dirac gamma matrices 7^, = 0, 1, 2, 3 that satisfy 
the anticommutation relations 

{7^7^} = 2r/^% (4.60) 

7]^^'^ being Minkowski's metric tensor. 

Similarly it can be generated by the following four observables, for instance: Si, S2, Ti and 
T2. In fact by ()4.56j) and ()4.59|1 we obtain S3 and T3 respectively and by ()4.57|1 . the remaining 
elements. 

Classically, the internal orientation of an electron is characterized by the knowledge of the 
components of the body frame {ei)j, i,j = 1,2,3 that altogether constitute an orthogonal 
matrix. To completely characterize in a unique way this orthogonal matrix we need at least 
four of these components. In the quantum version, the knowledge of four (ej)j matrices and by 
making use of ()4.52ll - ()4.59ll . allows us to recover the remaining elements of the complete Dirac 
algebra. It is in this sense that internal orientation of the electron completely characterizes its 
internal structure. Dirac's algebra of translation invariant observables of the electron can be 
generated by the orientation operators. 



Chapter 5 



Some spin features 



5.1 Gyromagnetic ratio 

The g = 2 gyromagnetic ratio of the electron was considered for years a success of Dirac's 
electron theory. ^ Later, Levy-Leblond ^ obtained similarly g — 2 but from a s = 1/2 non- 
relativistic wave equation. Proca ^ found = 1 for spin 1 particles and this led Belinfante ^ 
to conjecture that the gyromagnetic ratio for elementary systems is g = 1/s, irrespective of 
the value s of its spin. He showed this to be true for quantum systems of spin 3/2, and a 
few years later the conjecture was analyzed and checked by Moldauer and Case ^ to be right 
for any half-integer spin, and by Tumanov ^ for the value s = 2. In all these cases a minimal 
electromagnetic coupling was assumed. 

Weinberg ^ made the prediction g = 2 ioi the intermediate bosons of the weak interactions 
when analyzing the interaction of W bosons with the electromagnetic field by requiring a 
good high-energy behavior of the scattering amplitude. The discovery of the charged 
spin 1 bosons with g = 2, contradictory to Belinfante's conjecture, corroborated Weinberg's 
prediction and raised the question as to whether g — 2 for any elementary particle of arbitrary 
spin. 

Jackiw ^ has given another dynamical argument confirming that the gyromagnetic ratio of 
spin-1 fields is g = 2, provided a nonelectromagnetic gauge invariance is accepted. He also 
gives some ad hoc argument for s = 2 fields, consistent with the g = 2 prescription. 

Ferrara et al. ^ in a Lagrangian approach for massive bosonic and fermionic strings, by the 
requirement of a smooth fixed-charge M — > limit, get g — 2 as the most natural value for 
particles of arbitrary spin. However the only known particles which fulfill this condition are 
leptons and charged bosons, i.e., charged fermions and bosons of the lowest admissible 
values of spin. No other higher spin charged elementary particles have been found. 

The aim of this section, instead of using dynamical arguments as in the previous attempts, 

1 P.A.M. Dirac, Proc. Roy. Soc. London A117, 610 (1928). 

2 J.M. Levy-Leblond, Comm. Math. Pbys. 6, 286 (1967). 

3 A. Proca, Compt. Rend. 202, 1420 (1936); Journ. Phys. Radium, 49, 245 (1988). 

4 F.J. Belinfante, Phys. Rev. 92, 997 (1953). 

5 P. A. Moldauer and K.M. Case, Pbys. Rev. 102, 279 (1956). 
^ V.S. Tumanov, Sov. Phys. JETP, 19, 1182 (1964). 

S. Weinberg, in Lectures on Elementary Particles and Quantum Field Theory, edited by S. Deser, M. 
Grisaru and H. Pendleton, MIT press, Cambridge, MA (1970), p. 283. 
« R. Jackiw, Phys. Rev. D 57, 2635 (1998). 

^ S. Ferrara, M. Porrati and V.L. Telegdi, Phys. Rev. D 46, 3529 (1992). 
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is to give a kinematical description of the gyromagnetic ratio of elementary particles which 
is based upon the double content of their spin operator structure. 

The general structure of the quantum mechanical angular momentum operator in either 
relativistic or nonrelativistic approach is 

h 

j = rx-V + S = rxP + S, (5.1) 

where the spin operator is 

h 

S = ux -Vu + W, (5.2) 

and Vu is the gradient operator with respect to the velocity variables and W is a linear 
differential operator that operates only on the orientation variables a and therefore commutes 
with the other. For instance, in the p = ntan(tt/2) parametrization W is written as 

= I [Vp + p X V, + p(p ■ Vp)] . (5.3) 

The first part in ()5.2|) . related to the zitterbewegung spin, has integer eigenvalues because 
it has the form of an orbital angular momentum in terms of the u variables. Half-integer 
eigenvalues come only from the operator ()5.3|) . This operator W takes into account the change 
of orientation, i.e., the rotation of the particle. 

We have seen in either relativistic or non-relativistic examples that if the only spin content 
of the particle S is related to the zitterbewegung part Z = u x U, then the relationship 
between the magnetic moment and zitterbewegung spin is given by 

11= - kx — = Z, 5.4 

^ 2 dt 2m ' ^ ' 

i.e., with a normal up to a sign gyromagnetic ratio g = I. If the electron has a gyromagnetic 
ratio g = 2, this implies necessarily that another part of the spin is coming from the angular 
velocity of the body, but producing no contribution to the magnetic moment. 

Therefore for the electron, both parts W and Z contribute to the total spin. But the 
W part is related to the angular variables that describe orientation and does not contribute 
to the separation k between the center of charge and the center of mass. It turns out that 
the magnetic moment of a general particle is still related to the motion of the charge by the 
expression (15. 4|) . i.e., in terms of the Z part but not to the total spin S. It is precisely when we 
try to express the magnetic moment in terms of the total spin that the concept of gyromagnetic 
ratio arises. 

Now, let us assume that both Z and W terms contribute to the total spin S with their 
lowest admissible values. 

For Dirac's particles, the classical zitterbewegung is a circular motion at the speed of light 
of radius R = S/mc and angular frequency u = mc^ / S, in a plane orthogonal to the total spin. 
The total spin S and the Z part, are both orthogonal to this plane and can be either parallel 
or antiparallel. Let us define the gyromagnetic ratio hj Z = gS. For the lowest admissible 
values of the quantized spins z = 1 and w = 1/2 in the opposite direction this gives rise to a 
total s = 1/2 perpendicular to the zitterbewegung plane and then g = 2. 

For s = 1 particles the lowest possible values compatible with the above relative orientations 
are z = 2 and w = 1 in the opposite direction, thus obtaining again g = 2. The possibility 
z = 1 and w = is forbidden in the relativistic case because necessarily w 7^ to describe 
vector bosons with a multicomponent wave- function. 
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5.2 Instantaneous electric dipole 

The internal motion of the charge of the electron in the center of mass frame is a circle at 
the speed of light. The position of the charge in this frame is related to the total spin by eq. 
(Uni, i.e., 

k = — - Sxu, (5.5) 

where S is the total constant spin and u = dk/dt, with m = c is the velocity of the charge. In 
addition to this motion there is a rotation of a local frame linked to the particle that gives rise 
to some angular velocity, but this rotation has no effect on the electric dipole structure. (See 
Fig. 15.11 where the angular velocity and the local frame are not depicted). 




Figure 5.1: Electron charge motion in the CM. frame. 

Now, from the point of view of the center of mass observer, the particle behaves as though 
it has a magnetic moment related to the particle current by the usual classical expression 

1 /" I ■ j3 6 / dk 
u = - / fc X 7 d r = - k X —, 

2 dV 

where e is the charge and jir — k) = e dk/dt (5^(r — k) is the particle current density. The 
orbital term kx dk/dt is related to the zitterbewegung part of spin that quantizes with integer 
values and which for spin 1/2 and spin 1 charged particles is twice the total spin S, giving rise 
to a pure kinematical interpretation of the gyromagnetic ratio g = 2 for this model as seen in 
the previous section. 

But also in the center of mass frame the particle has an oscillating instantaneous electric 
dipole moment d = ek, that is thus related to the total spin by 

d=-^Sxu. (5.6) 

This instantaneous electric dipole, which fulfills the usual definition of the momentum of 
the point charge e with respect to the origin of the reference frame, is translation invariant 
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because it is expressed in terms of a relative position vector k. It can never be interpreted as 
some kind of fluctuation of a spherical symmetry of a charge distribution. Even in this kind 
of model, it is not necessary to talk about charge distributions, because all particle attributes 
are defined at single points. 

In his original 1928 article, Dirac obtains that the Hamiltonian for the electron has, in 
addition to the Hamiltonian of a free point particle of mass m, two new terms that in the 
presence of an external electromagnetic field are 

eh ieh 

1: B + a E = -fx B - d E, 5.7 

2m 2mc 



where 

a- 
cr 



and CK = 7o7, 



i.e., Xl is expressed in terms of cr Pauli-matrices and a is Dirac's velocity operator when written 
in terms of Dirac's gamma matrices. 

We shall show that the quantum counterpart of expression is in fact the electric 

dipole term of Dirac's Hamiltonian ()5.7|l . The remaining part of this section is to consider the 
representation of the 'cross' product in ()5.6p in terms of the matrix (or geometric) product of 
the elements of Dirac's algebra that represent the quantum version of the above observables, 
so that a short explanation to properly interpret these observables as elements of a Clifford 
algebra is given in what follows. 

Both, velocity operator u = ca and spin operator S are bivectors in Dirac's algebra, 
considered as elements of the Geometric or Clifford algebra of space-time in the sense of 
Hestenes. 

In fact, Dirac's alpha matrices are written as a product of two gamma matrices = 7o7i 
and also the spin components Sj = {ih/2) •jkli, j, k, I cyclic 1, 2, 3, and where the four gamma 
matrices, 7^, = 0, 1,2,3 are interpreted as the four basic vectors of Minkowski's space-time 
that generate Dirac's Clifford algebra. They satisfy 7^ ■ 7,/ = rj^^,, i.e., 7q = 1 and •yf = — 1, 
where the dot means the inner product in Dirac's Clifford algebra. We thus see that velocity 
and spin belong to the even subalgebra of Dirac's algebra and therefore they also belong to Pauli 
algebra or geometric algebra of three-dimensional space. Under spatial inversions 70 70 and 
7i — 7j, the velocity operator changes its sign and it is thus a spatial vector, while the spin 
is invariant under this transformation as it corresponds to a spatial bivector or pseudovector. 

The relationship between the cross product and the outer and inner product of two vectors 
a and b in Pauli algebra is, 

a X b = —ia Ab = b ■ (ia), (5.8) 

where A represents the symbol for the outer product in geometric algebra, the imaginary unit 
i represents the unit three-vector or pseudoscalar and ia is the dual bivector of vector a. 

The inner product of a vector b and a bivector A is expressed in terms of the geometric 
product in the form 

b-A=^{bA-Ab) (5.9) 

where in Dirac's or Pauli algebra the geometric product bA is just the ordinary multiplication 
of matrices. 



" P.A.M. Dirac, Proc. Roy. Soc. London, A117, 610 (1928). 

D. Hestenes, Space-Time algebra, Gordon and Breach, NY (1966); D. Hestenes and G. Sobczyk, Clifford 
Algebra to Geometric Calculus, D. Reidel Pub. Co. Dordrecht, (1984). 
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7z7s 



(a) 



Figure 5.2: A basis for vectors (o) and bivectors (pseudovectors) (6) of Pauli algebra. 

If we choose a basis of vectors and pseudovectors as in Fig. 15. 2[ where the double-lined 
objects of part (6) represent the dual vectors of the corresponding spatial bivectors, and express 
in these bases the observables of Fig. 15. H then the spatial velocity vector u = C7072 and the 
pseudovector S = {h/2)j2l3 and therefore, using (I5.8jl and (j5.9|) we get 



icfi ( 1 

S xu = u- (iS) = — - (70727273 - 72737o72) 



-ich 



-7o73- 



Now vector k = -R7073 with R = h/2mc, and substituting in (15. 6j) we get the desired result. 



5.3 Classical Tunneling 

As a consequence of the zitterbewegung and therefore of the separation between the center 
of mass and center of charge, we shall see that spinning particles can have a non-vanishing 
crossing of potential barriers. 

Let us consider a spinning particle with spin of (anti)orbital type, as described in Section 
12. 2L under the influence of a potential barrier. The Langrangian of this system is given by: 



^ mr 
^ 2"T 



m u 



2cj2 t 



eV{r)i. 



(5.10) 



Sharp walls correspond classically to infinite forces so that we shall consider potentials that 
give rise to finite forces like those of the shape depicted in Fig. 15.31 where Vq represents the 
top of the potential. 

Then the external force F{x), is constant and directed leftwards in the region x G (— a, 0) 
and rightwards for x G (0,6), vanishing outside these regions. 

Potentials of this kind can be found for instance in the simple experiment depicted in Figure 
15.41 in which an electron beam, accelerated with some acceleration potential K, is sent into the 
uniform field region of potential Vq contained between the grids or plates A, C and B. 

In Figure 15.41 from a strict classical viewpoint a spinless electron stops at the dotted line 
and is rejected backwards. But a classical spinning electron can cross the barrier provided its 
kinetic energy is above some minimum value, although below the top of the potential. This 
minimum value depends on the separation between plates. 
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Figure 5.3: Triangular potential barrier. 
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Figure 5.4: Electron beam into a potential barrier. A classical spinless electron never 
crosses the dotted line. A spinning particle of the same energy might cross the barrier. 
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Let us assume for simplicity that the spin is pointing up or down in the z direction such 
that the point charge motion takes place in the XOY plane. Let Qx, Qy and qz = 0, be the 
coordinates of the center of mass and x, y and z = 0, the position of the charge. 

The dynamical equations are 

''''-'f(.), ^1 = 0. (5.11) 



dt'^ m ' dt"^ 

dt- ' " (^-^^)=0' ^ 



^ =0, ^ + = 0, (5.12) 



where 




for X G (—a, 0), 
F{x) = { eVo/b, for x E (0,6), 
otherwise. 

Equations ()5.1HI are nonlinear and we have not been able to obtain an analytical solution in 
closed form. We shall try to find a numerical solution. To make the corresponding numerical 
analysis we shall define different dimensionless variables. Let R be the average separation 
between the center of charge and center of mass. In the case of circular internal motion, it 
is just the radius Rq of the zitterbewegung. Then we define the new dimensionless position 
variables: 

qx = qx/R, qy = qy/R, x = x/R, y = y/R, d = a/R, b = b/R. 

The new dimensionless time variable a = ut is just the phase of the internal motion, such 
that the dynamical equations become 



d'^x ^ ^ d'^y 



where A{x) is given by 



-eVo/amu;^i?^, for x G (— a,0), 
A{x) = { eVo/bmcu^R^, for x G (0,6), 
0, otherwise. 

In the case of the relativistic electron, the internal velocity of the charge is uR = c, so that 
the parameter e/mc^ = 1.9569 x 10~^V~^, and for potentials of order of 1 volt we can take the 
dimensionless parameter eVo/mu^R^ = 1.9569 x 10~^. 

If we choose as initial conditions for the center of mass motion 

g^(0) = 0, dqy{0)/da = 0, 

then the center of mass is moving along the OX axis. The above system reduces to the analysis 
of the one-dimensional motion where the only variables are (jx and x. Let us call from now 
on these variables q and x respectively and remove all hats from the dimensionless variables. 
Then the dynamical equations to be solved numerically are just 

d^Q d'^x 

5^+.-« = 0. (5.13) 
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where A{x) is given by 

C -1.9569 X lO"'' a' Vo, forxG(-a,0), 

A{x) = < 1.9569 X IQ-H-^Vo, for x G (0,6), (5.14) 

[ 0, otherwise. 

Numerical integration has been performed by means of the computer package Dynamics 
Solver. The quahty of the numerical results is tested by using the different integration 
schemes this program allows, ranging from the very stable embedded Runge-Kutta code of 
eight order due to Dormand and Prince to very fast extrapolation routines. All codes have 
adaptive step size control and we check that smaller tolerances do not change the results. 
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Figure 5.5: Kinetic Energy during the crossing for the values a = b = 1. 

With a = b = 1, and in energy units such that the top of the barrier is 1, if we take an 
initial kinetic energy K below this threshold, K = mg(0)^/2eVo = 0.41 we obtain for the center 
of mass motion the graphic depicted in Fig. 15. 5| where is shown the variation of the kinetic 
energy of the particle K{q), with the center of mass position during the crossing of the barrier. 
There is always crossing with a kinetic energy above this value. In Fig. l5.6[ the same graphical 
evolution with a = 1 and 6 = 10 and K = 0.9055 for a potential of 10'^ Volts in which the 
different stages in the evolution are evident. Below the initial values for the kinetic energy of 
0.4 and 0.9 respectively, the particle does not cross these potential barriers and it is rejected 
backwards. 

If in both examples the parameter a is ranged from 1 to 0.05, thus making the left slope 
sharper, there is no appreciable change in the crossing energy, so that with a = 1 held fixed 
we can compute the minimum crossing kinetic energies for different b values, Kc{b). 

To compare this model with the quantum tunnel effect, let us quantize the system. In the 
quantization of generalized Lagrangians developed in the Chapter O the wave function for this 
system is a squared-integrable function ip{t,r,u), of the seven kinematical variables and the 
generators of the Galilei group have the form: 

d 

H = th—, P = -ihV, K = mr-tP + ihWu, J = r x P + Z, (5.15) 

where is the gradient operator with respect to the u variables. These generators satisfy 

J.M. Aguirregabiria, Dynamics Solver, computer program for solving different kinds of dynamical systems, 
which is available from his author through the web site < http : 77tp . Ic . ehu . es/ j ma . html'> at the server of 
the Theoretical Physics dept. of The University of the Basque Country, Bilbao (Spain). 
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Figure 5.6: Kinetic Energy during the crossing for the values a = 1, b = W. 



the commutation relations of the extended Gahlei group, ^'^ and the spin operator is given by 
Z = —ihu X Vn. 

One Casimir operator of this extended Gahlei group is the Galilei invariant internal energy 
of the system which in the presence of an external electromagnetic field and with the minimal 
coupling prescription is written as, 



£ = H-eV 



1 

2m 



(P - eAf 



(5.16) 



where V and A are the external scalar and vector potentials, respectively. 

In our system A = 0, and V is only a function of the x variable. It turns out that because 
of the structure of the above operators we can find simultaneous eigenfunctions of the following 
observables: the Casimir operator (|5.16p . H, Py, P^, Z"^ and Z^. The particle moves along 
the OX axis, with the spin pointing in the OZ direction, and we look for solutions which are 
eigenfunctions of the above operators in the form: 



H - eV(x) - — p2 



^ = Hip = Ftp, Pytp = 0, P,ip = 0, 



(5.17) 



Z^iP = s{s + l)h^i), Z,il) = ±shil), (5.18) 

so that ip is independent of y and z, and its time dependence is of the form exp{—iEt/h). Since 
the spin operators produce derivatives only with respect to the velocity variables, we can look 
for solutions with the variables separated in the form: 



x)x\u) 



and thus 



+ E - eV{x) -E\(p{x 



0, 



(5.19) 



_ 2m dx^ 

Z\{u) = s{s + l)h\{u), Z^xH = ±shx{u), (5.20) 

where the spatial part 0(x), is uncoupled with the spin part x(**)i ^^nd E — eV{x)—S represents 
the kinetic energy of the system. The spatial part satisfies the one-dimensional Schroedinger 



J.M. Levy-Leblond, Galilei Group and Galilean Invariance, in E.M. Loebl, Group Theory and its applica- 
tions, Acad. Press, NY (1971), vol. 2, p. 221. 
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equation, and the spin part is independent of the interaction, so that the probabihty of quantum 
tunneling is contained in the spatial part and does not depend on the particular value of the 
spin. If the particle is initially on the left-hand side of the barrier, with an initial kinetic 
energy Eq = E — £, then we can determine the quantum probability for crossing for a = 1 and 
different values of the potential width b. 

The one-dimensional quantum mechanical problem of the spatial part for the same one- 
dimensional potential depicted in Fig. 15.31 is: 



+ Re 
CiAi{D{l 
C3Ai(L(l - 



-ikx 



G- 



X\ 

a) ~ 



G2Qi{D{l 
- C4Bi(L(l 



G-^ 



)), 



X < —a, 
—a < X < 0, 
0<x<b, 
X > b. 



(5.21) 



where x is the same dimensionless position variable as before, and the constants 



k 



E 



D 



eV^a? 
2mc^ ' 



eVob^ 



2mc^ 



G 



E 



(5.22) 



Functions Ai(x) and Bi(x) are the Airy functions of x. The six integration constants R, T, and 
Ci^i = 1,2,3,4, can be obtained by assuming continuity of the functions and their first order 
derivatives at the separation points of the different regions. The coefficient represents the 
probability of the particle to be refiected by the potential and |Tp its probability of crossing. 
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Figure 5.7: Classical and Quantum Probability of crossing for different potentials. 

Computing the T amplitude for a = 1 and different values of the potential width b, and 
for energies below the top of the barrier eVo, we show in Fig. 15. 7[ the average probability for 
quantum tunneling for four different potentials Vq of 10^, 10^, 10^ and 10^ Volts. This average 
probability has been computed by assuming that on the left of the barrier there is a uniform 
distribution of particles of energies below eVb. 

L. Landau and E. Lifchitz, Mecanique quantique, Mir Moscow (1988), 3rd. edition. 
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If we consider for the classical spinning particle the same uniform distribution of particles, 
then, the function P{b) = 1 — Kc{b), where Kc{b) is the minimum dimensionless kinetic energy 
for crossing computed before, represents the ratio of the particles that with kinetic energy 
below the top of the potential cross the barrier because of the spin contribution. 

This function P{b), is also depicted in Fig. 15.71 We see that for the different potentials 
shown in that figure the classical average probability of crossing is smaller than the quantum 
one, but for stronger potentials this classical probability, coming from the spin contribution, 
becomes relatively important. 

Because the tunnel effect is a function of h and the spin of elementary particles is also 
of order of h it is very difficult to separate from the outcome of a real experiment involving 
elementary particles, which part is due to a pure quantum effect and which is the contribution 
to crossing coming from the spin structure. From ()5.19|1 and ()5.2np it is clear that the quantum 
probability of tunneling is independent of the spin. 

To test experimentally this contribution, it will be necessary to perform separate experi- 
ments with particles of the same mass and charge but with different values of the spin. Thus, 
the difference in the outcome will be related to the spin contribution. This can be accomplished 
for instance, by using ions of the type A'^^ that could be either in a singlet, (s = 0) state or 
in a triplet (s = 1) state. 

But if there exists a contribution to crossing not included in the usual quantum mechanical 
analysis we have to modify the quantum mechanical equations. To be consistent with the 
above analysis the Schroedinger-Pauli equation should be modified to include the additional 
electric dipole term. A term of the form —eERcosut, where E is the external electric field and 
R the radius of the zitterbewegung, should be considered to solve the corresponding quantum 
wave function. This term is of the order of the separation R between the center of mass and 
center of charge, which is responsible for the classical crossing. This additional electric dipole 
term is already included in Dirac's equation but is suppressed when taking the low velocity 
limit, as it corresponds to this low energy example. Nevertheless, although this is a low energy 
process and the time average value of the electric dipole vanishes, there are very high field 
gradients. 

We see that the separation between the center of mass and center of charge that gives 
rise to the spin structure of this particle model justifies that this system can cross a potential 
barrier even if its kinetic energy is below the top of the potential. 

5.3.1 Spin polarized tunneling 

I like to point out the following ideas to discuss whether they can be useful in connection 
with the interpretation of the magnetoresistance of polycrystaline films. This is known in the 
literature as the spin polarized tunneling. 

The main feature of the "classical" spin polarized tunneling we have seen in the previous 
section is not a matter of whether tunneling is classical or not, because this is a nonsense 
question. Matter at this scale is interpreted under quantum mechanical rules. But if we use a 
model of a classical spinning particle that, when polarized orthogonal to the direction of motion, 
produces a crossing that is not predicted by the Schroedinger-Pauli equation, it means that 
this quantum mechanical equation is lacking some term. The coupling term — /x ■ B, between 
the magnetic moment and magnetic field that gives rise to the Pauli equation, is inherited 

I'' V.N. Dobrovolsky, D.I. Sheka and B.V. Chernyachuk, Surface Science 397, 333 (1998); P. Raychaudhuri, 
T.K. Nath, A.K. Nigam and R. Pinto, cond-mat/9805258, preprint. 
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from Dirac's electron theory. But Dirac's equation also predicts another term —d ■ E, of the 
coupling of an instantaneous electric dipole with the electric field. It is this oscillating electric 
dipole term that we believe is lacking in quantum mechanical wave equations. In general, the 
average value of this term in an electric field of smooth variation is zero. But in high intensity 
fields or in intergranular areas in which the effective potentials are low, but their gradients 
could be very high, this average value should not be negligible. 

The conduction of electrons in synterized materials is completely different than the con- 
duction on normal conductors. The material is not a continuous crystal. It is formed by 
small grains that are bound together by the action of some external pressure. If we can depict 
roughly the electric current fiow, this is done by the jumping of electrons from grain to grain, 
through a tunneling process in which there is some estimated effective potential barrier con- 
fined in the gap between grains. Therefore these materials show in general a huge resistivity 
when compared with true conductors. 

The form of this potential is unknown. The simplest one is to assume a wall of thickness d, 
the average separation between grains, and height h. But it can also be estimated as one of the 
potentials of the former example. What we have shown previously is that for every potential 
barrier, there is always a minimum energy, below the top of the potential, that electrons above 
that energy cross with probability 1 when polarized orthogonal to the motion, even within 
a classical interpretation. But this effect is not predicted by "normal" quantum mechanics 
because tunneling is spin independent. 

Now, let us assume that we are able to estimate some average effective potential barrier in 
the intergranular zone of this polycristaline material. If the corresponding minimum crossing 
energy of this barrier for polarized electrons is below the Fermi level, then, when we introduce 
a magnetic field in the direction of the film and the magnetic domains in the grains become 
polarized, all electrons above that minimum energy of crossing will fiow from grain to grain as 
in a good conductor, with a classical probability 1. That's all. Here the difficulty is to estimate 
properly this potential barrier and therefore the corresponding classical crossing energy. 

It can be argued that the presence of the magnetic field to polarize electrons produces a 
change in the energy of particles. Nevertheless, even for a magnetic field of the order of 1 Tesla 
and in a potential barrier of 1 Volt, the magnetic term — ^t ■ B contributes with an energy of 
order of ±5.7 x 10~^eV, which does not modify the quantum probability of crossing. 

5.4 Formation of bound pairs of electrons 

If we have the relativistic and nonrelativistic differential equations satisfied by the charge of 
the spinning electrons we can analyse as an example, the interaction among them by assuming 
a Coulomb interaction between their charges. In this way we have a system of differential 
equations of the form ()2.1481 12.149|) in the relativistic case or in the form ()2.7H I2.72p in the 
nonrelativistic, for each particle. For instance, the external field acting on charge ei is replaced 
by the Coulomb field created by the other charge 62 at the position of ei, and simmilarly for the 
other particle. The integration is performed numerically by means of the numerical integration 
program Dynamics Solver 

In Figure Em we represent the scattering of two spinning electrons analysed in their center 

^^J.M. Aguirregabiria, Dynamics Solver, computer program for solving different kinds of dynamical systems, 
which is available from his author through the web site < http : //tp . Ic . ehu . es/ j ma . htmr> at the server of 
the Theoretical Physics Dept. of The University of the Basque Country, Bilbao, Spain. 
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Figure 5.8: Scattering of two spinning electrons with the spins parallel, in their center 
of mass frame. It is also depicted the scattering of two spinless electrons with the same 
energy and linear momentum. 

of mass frame. We send the particles with their spins parallel and with a nonvanishing impact 
parameter. In addition to the curly motion of their charges we can also depict the trajectories of 
their center of masses. If we compare this motion with the Coulomb interaction of two spinless 
electrons coming from the same initial position and with the same velocity as the center of mass 
of the spinning electrons we obtain the solid trajectory marked with an arrow. Basically this 
corresponds to the trajectory of the center of mass of each spinning particle provided the two 
particles do not approach each other below Compton's wave length. This can be understood 
because the average position of the center of charge of each particle aproximately coincides 
with its center of mass and as far as they do not approach each other too much the average 
Coulomb force is the same. The difference comes out when we consider a very deep interaction 
or very close initial positions. 

In Figure Em we represent the initial positions for a pair of particles with the spins parallel. 
The initial separation a of their center of masses is a distance below Compton's wave length. 
We also consider that initially the center of mass of each particle is moving with a velocity v 
as depicted. That spins are parallel is reflected by the fact that the internal motions of the 
charges, represented by the oriented circles that surround the corresponding center of mass, 
have the same orientation. It must be remarked that the charge motion around its center of 
mass can be characterised by a phase. The phases of each particle are chosen opposite to each 
other. We also represent the repulsive Coulomb force F computed in terms of the separation 
of the charges. This interacting force F has also been attached to the corresponding center of 
mass, so that the net force acting on point m2 is directed towards point mi, and conversely. We 
thus see there that a repulsive force between the charges represents an atractive force between 
their center of masses when located at such a short distance. 

In FigureEHniwe depict the evolution of the charges and masses of this two-electron system 
for a = 0.4Ac and v = 0.004c during a short time interval. Figure EHU represents only the 
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Figure 5.9: Initial position and velocity of the center of mass and charges for a bound 
motion of a two-electron system with parallel spins. The circles would correspond to 
the trajectories of the charges if considered free. The interacting Coulomb force F is 
computed in terms of the separation distance between the charges. 

motions of the center of masses of both particles for a longer time. It shows that the center of 
mass of each particle remains in a bound region. 




Figure 5.10: Bound motion of two electrons with parallel spins during a short period of 
time 



The evolution of the charges is not shown in this figure because it blurs the picture but it 
can be inferred from the previous figure. We have found bound motions at least for the range 
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< a < 0.8Ac and velocity < f < 0.01c. We can also obtain similar bound motions if the 
initial velocity v has a component along the OX axis. The bound motion is also obtained for 
different initial charge positions as the ones depicted in Figure EIHl This range for the relative 
phase depends on a and v but in general the bound motion is more likely if the initial phases 
of the charges are opposite to each other. 

We thus see that if the separation between the center of mass and center of charge of a 
particle (zitterbewegung) is responsible for its spin structure then this atractive effect and also 
a spin polarised tunneling effect can be easily interpreted. 



m,(0) 






\ mi(0) 



Figure 5.11: Evolution of the center of mass of both particles for a larger time 

A bound motion for classical spinless electrons is not possible. We can conclude that one 
of the salient features of this example is the existence from the classical viewpoint of bound 
states for spinning electron-electron interaction. It is the spin structure which contributes 
to the prediction of new physical phenomena. If two electrons have their center of masses 
separated by a distance greater than Compton's wave length they always repel each other as 
in the spinless case. But if two electrons have their center of masses separated by a distance less 
than Compton's wave length they can form from the classical viewpoint bound states provided 
some initial conditions on their relative initial spin orientation, position of the charges and 
center of mass velocity are fulfilled. 

The example analysed gives just a classical prediction, not a quantum one, associated 
to a model that satisfies Dirac's equation when quantised. The possible quantum mechanical 
bound states if they exist, must be obtained from the corresponding analysis of two interacting 
quantum Dirac particles, bearing in mind that the classical bound states are not forbiden from 
the classical viewpoint. Bound states for a hydrogen atom can exist from the classical viewpoint 
for any negative energy and arbitrary angular momentum. It is the quantum analysis of the 
atom that gives the correct answer to the allowed bound states. 
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